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A Spark Shadowgraphic Study of Body Waves in Water* 


J. Howarp McMILLEN** AND E. NEwTron HARVEY 
Department of Biology, Princeton University, Princeton, New Jersey 


(Received April 17, 1946) 


Shadows of shock waves in water were photographed by 
means of a point-spark of less than one-tenth micro- 
second duration; the waves were formed when high ve- 
locity spheres struck the water surface. The shock wave 
front had the form of an elliptical arc of small eccentricity; 
the ellipticity was caused by the greater strength and 
velocity of the wave at its center. The absence-of-light 
band had a varying width which depended on the strength 
of the wave. A large number of secondary waves behind 
the impact wave were observed; these were produced by 
the vibrating spherical missile as it traveled through the 
water. The period of vibration was determined from the 
spacing of the waves for 4/32’, 6/32”, and 8/32” spheres; 
these agreed to within 19 percent of the periods calculated 
on the basis of Lamb's theory. Interference of the vibration 
waves was also observed when two spheres were shot into 
the water simultaneously. The cavity behind the missile 
was nearly conical in shape with the forward tip being 
distorted by pressure changed immediately ahead of the 
missile. A sphere produced a cusp at the forward tip of 
the cavity; cylinders and cones each produced their own 


HE spark shadowgraph is one of the most 

useful methods of studying transient phe- 
nomena. In this method a point-source spark of 
microsecond duration is employed to record 
photographically the shadows of moving ob- 
jects. Not only is a sharp silhouette of the 
moving objects reproduced, but any changes of 


* This work has been carried out under a contract 
recommended ,by the Committee on Medical Research 
between the Office of Scientific Research and Development 
and Princeton University. 

** Now in the Research Department, Naval Ordnance 
Laboratory, Washington, D. C, 


characteristic shadow pattern. When a grid was placed in 
front of the tank, distortions in its shadowgram were 
observed. These were caused by the high pressures sur- 
rounding the missile and, from the shift in the grid pattern, 
the direction and magnitude of the pressure gradient 
could be inferred. When the shock waves were reflected 
with the wave front nearly normal to the surface there was 
an indication that the stem, which is characteristic of 
Mach reflections, was present. All other reflections were 
regular. Focusing of the wave by a brass mirror was 
demonstrated in the spark shadowgram. Waves reflected 
near normal incidence from substances having a lower 
acoustic impedance than water were reflected with a 
presence-of-light band ahead of the usual absence-of-light 
band. At glancing incidence this reversal occurred even 
when the reflecting substance had a higher acoustic 
impedance. Striations frequently appeared in the shadow- 
grams when waves were reflected. These striae were from 
0.02 to 0.08 cm apart and represented a quasi-stationary 
state. They were believed to be formed when the reflected 
vibration waves passed through the incident waves. 


density that might occur in the medium are also 
recorded. Density changes may be produced by 
waves of displacement, by regional alterations 
in temperature, and by the mixing ‘of different 
liquids. The first spark shadows were produced 
by Dvorak' in 1880; heated gases from a flame 
and striations in optical glass were some of the 
shadows studied. C. V. Boys,? in 1893, photo- 
graphed the shadows of fast moving bullets and 
the acoustic waves which accompanied the mis- 





!V. Dvorak, Ann. d. Physik 9, 502 (1880). 
2C. V. Boys, Nature 47, 415 (1893). 
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siles. R. W. Wood,’ in 1900, and Foley,‘ in 1912, 
obtained shadowgrams showing the reflection 
and refraction of the acoustic waves in air, the 
waves being produced by sparks. Since then 
many investigators have used the spark shadow- 
graph to study both the permanent and transient 
changes in fluid density. 

A report recently released for publication by 
the U. S. Navy Department reprints several 
spark shadowgrams showing waves produced in 
water by explosions. These photographs were 
made by P. Libessart and appear in an article 
by H. Polacheck and R. J. Seeger.® Schlieren 
photographs of waves in water have also been 
made; these were spark excited and were photo- 
graphed by O. v Schmidt.* Other photographs 
of body waves in water, using processes similar 
to spark photography, have been made by: D. P. 
MacDougall, C. H. Messerly, and E. M. Boggs; 
E. C. Campbell and C. W. Wyckoff; R. W. 
Spitzer and R. S. Price; and A. M. Shanes.5 

In the present investigation the spark shadow- 
graph was employed to study body waves in 
water when these waves were produced by the 
impact of an object on the water surface. Some 
of the results of this investigation have already 
been published.’ In all experiments the wave 
was produced by the impact of a fast moving 
sphere on the water surface. Many of these 
impact waves were strong enough to be shock 
waves and were observed to possess the usual 
characteristics of such waves. Others of these 
waves were either weak shock waves or simply 
acoustic waves, it being very difficult to distin- 





Fic. 1. Sketch of apparatus used to produce spark 
shadowgraphs in water. 


*R. W. Wood, Phil. Mag. 48, 218 (1899). 
(94s L. Foley and W. H. Souder, Phys. Rev. 35, 373 
). 
5H. Polachek and R. J. Seeger, Navy Bur. of Ord. Exp. 
Res. Rep. No. 14, Aug. (1944). 
® QO. v Schmidt, Zeits. f. tech. Physik 19, 554 (1938). 
7 J. Howard McMillen, Phys. Rev. 68, 198 (1945). 
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guish between these two wave types when they 
occur in water. 


1. APPARATUS AND METHOD 


The illuminating spark was patterned after a 
spark which was developed jointly by the Gen- 
eral Electric Company and Aberdeen Proving 
Ground. The spark was produced between two 
magnesium electrodes by discharging a condenser 
at 4400 volts having a 0.2 uf capacitance. The 
light from the spark traversed a tank of water, 
T, and fell on the photographic plate, PP, as 
shown in Fig. 1. Almost all the waves which 
were shadowgraphed were produced by small 
spheres that were shot vertically downward into 
the water, following the path, G, in the diagram. 
The spheres struck the water with velocities 
between 1000 and 5000 ft./sec. and their sizes 
varied from. 1/8” to 1/4” in diameter. The 
largest kinetic energy possessed by these missiles 
was about 1500 ft. lb. or 2.0310" ergs. 

The spark was triggered after the wire screen, 
O, was broken by the missile. The time interval 
between breaking the screen and ignition of the 
spark was controlled by an electronic delay sys- 
tem. Shadowgraphs of spheres moving through 
air with velocities of around 5000 ft.//sec. showed 
very little blurring, indicating that the effective 
duration of the spark did not exceed a tenth of 
a microsecond. 

The water tank had two Plexiglas walls which 
were either 1/4” or 1/2” thick, depending upon 
the size sphere being shot. The distance between 
the Plexiglas walls was 53} and between the 
other two walls it was 143”. The water was 
103” deep. The shadowgrams were recorded on 
11X14-inch lantern slide plates of medium 
contrast. 

The shadowgrams gave a magnification of 1.11 
diameters to a central object in the tank. The 
spark was placed in either one of two positions. 
In the first, the axial ray from the spark passed 
through the water surface and was perpendicular 
to the Plexiglas walls and missile path. In the 
other, the axial ray was lowered so as to be 10 
cm below its previous position. Because of the 
refraction of the rays as they passed from the 
water to the air, two rays originating in the tank 
and striking the plate would be about 18.2 
percent farther apart than they would be if 
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there were no refraction. Further details con- 
cerning the apparatus are given elsewhere.” 


2. WATER IMPACT WAVES 


When a fast moving object strikes a water 
surface a high pressure pulse is formed. This 
pulse or wave travels out from the point of 
impact along radial lines and its front bounds a 
region of high pressure, high density, and high 
refractive index. When the impacting object 
transfers a large amount of momentum to the 
water, the wave front marks a discontinuous 
change in pressure and the wave is called a shock 
wave. The shadowgram of Fig. 2 shows such a 
wave; the abrupt density change across the 
wave front produces the light and dark bands as 
it is recorded in the spark shadowgram. The 
wave was formed by a 4/32” steel sphere striking 
the water with a velocity of 4840 ft./sec. and 
at the time the wave was shadowgraphed the 
peak pressure was about 141 atmospheres. At 
the front of the wave the refractive index was 
increased by 0.14 percent while the density was 
raised by 0.64 percent. 

The width of the leading absence-of-light or 
dark band indicates the strength of the wave 
and can be used to measure the wave’s peak 
pressure. This peak pressure, as was shown by 
Hilton,’ is proportional to the three-halves power 
of the band width and inversely proportional to 
the square root of the wave’s radius of curvature. 
Inspection of the band width along the spherical 
wave front reveals’ that the pressure is greatest 
directly ahead of the sphere and that it tapers 
off along the arc as the surface of the water is 
approached. This variation in pressure is con- 
sistent with the angular distribution of the linear 
impulse delivered by the sphere to the water at 
entrance. 

Shock waves in air are distinguished from 
acoustic or sound waves by the fact that their 
velocity increases with their strength or peak 
pressure. The same is true in water but to a 
lesser degree. In the spark shadowgram this 
dependence of wave velocity on strength is 
revealed by the shape of the wave front. The 
center portion, being stronger, travels faster than 
the portion near the surface so that the wave is 





8 W. F. Hilton, Proc. Roy. Soc. A169, 174 (1938). 
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Fic. 2. Spark shadowgram of shock wave produced by 
4/32” steel sphere. Wave had advanced 10.7 cm into the 
water. Debris from trigger causes weak waves in front of 
main wave. 


not hemispherical but rather ellipsoidal. The 
long axis of the elliptical wave front shown in 
Fig. 2 is 3.4 percent larger than the short axis, 
showing that the average velocity of the central 
portion was supersonic by the same percent; 
the weak portion of the wave along the surface 
travels with the normal sound velocity. The 
velocity of sound in water at 15°C is 4789 ft./sec. 

The shock wave in water eventually weakens 
and travels with a constant velocity, the acoustic 
velocity. The sphere on the other hand is 
retarded by forces proportional to the square of 
its velocity. Just after entering the water the 
sphere in Fig. 2 is subjected to a negative 
acceleration equal to two million g. As a result, 
the wave moves farther and farther away from 
the sphere. The sphere in Fig. 2 entered the 
water with a velocity of 4760 ft./sec., but at the 
time the shadowgram was made the velocity 
had dropped to 2570 ft./sec. The sphere produces 
a cavity which is opaque to the light from the 
spark and which extends all the way to the 
surface. Above the surface a splash is shown 
moving vertically upward with a velocity which, 
in Fig. 2, is about 1600 ft./sec. This velocity is 
supersonic in air and the shock wave produced 
by the splash is clearly visible, as are the shock 
waves formed in air by the sphere in its down- 
ward flight. 

By measuring the band widths of shock waves 
formed by spheres having different diameters 
and impact velocities, it was found’ that the 
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Fic. 3. Spark shadowgram showing secondary waves 
»roduced by vibrations of a 4/32” steel sphere. The wave 
as advanced 4.0 cm into the water. The bottom of the 

horizontal black band is the water surface, the black band 
being a shadow of the meniscus on the Plexiglas wall. 


pressure in the wave varies as the 2.17 power of 
the impact velocity and is a linear function of 
the projection area of the sphere. The pressure 
which is developed at impact depends on how 
much water is compressed between the surface of 
the entering sphere and the outgoing shock wave. 


3. SECONDARY WAVES FORMED BY 
VIBRATING MISSILE 


When a missile strikes the water it is generally 
set into vibration and produces a series of 
acoustic waves in the water as it proceeds along 
its trajectory. These waves are found in the 
spark shadowgram just behind the main shock 
wave. Those formed by spherical missiles have 
been studied more extensively because of the 
sphere’s symmetry. 

A good example of secondary waves formed 
by a 4/32” steel sphere is shown in Fig. 3 
where the waves completely fill the region be- 
tween the vibrating sphere and the leading wave. 
These waves or bands are circular arcs with 
centers which lie between the surface of the 
water and the sphere, an observation which 
shows that they originate with the sphere as it 
-moves forward in the water. 

The bands generally form a continuous line 
throughout the arc covering an angular sweep 
of from 160° to 170°. However, it is not uncom- 
mon to find many bands disappearing or be- 
coming indistinct along their arc. Some bands 


544 








change along the arc from sharp to diffuse or 
from diffuse to sharp. This change seems to take 
place most frequently near the angle of about 
45°, as measured from the water surface. 

It is interesting to note that those secondary 
waves or bands which lie close to the missile are 
not perfect arcs of a circle but are distorted. 
This distortion is caused by the change in the 
refractive index near the missile which is pro- 
duced by the large pressure gradient there. 

In spite of the complex nature of these bands 
a definite band pattern can be discerned. This 
pattern is more distinct directly in front of the 
missile; it is also easier to identify after the wave 
has advanced into the tank an appreciable dis- 
tance. This pattern consists of a series of sharp 
absence-of-light bands with a broad absence-of- 
light band in between them. These dark bands 
are flanked by a narrower presence-of-light band. 
Figure 4 shows the band pattern produced by a 
4/32” and 6/32” steel sphere. One observes that 
the bands of the smaller spheres are closer to- 
gether. When the broad dark band is examined 
in the shadowgrams, as one moves away from 
the center along the arc, one observes that it is 
really two bands which only become resolved 
away from the center. Here it is observed to 
resolve itself into two narrow bands and thus 
the repetitive pattern is really composed of three 
bands, one sharp band and two bands which 
remain unresolved except at the side. 

The separation between the sharp bands was 


Fic. 4. Spark shadow- 
grams of the secondary 
waves produced by a vi- 
brating spherical missile. 
The one on the left was 
produced by a 4/32” steel 
sphere; theoneon theright 
by a 6/32” steel sphere. 
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measured for those waves which are in a line 
with the flight of the sphere. These are recorded 
in Table | for the 4/32”, 6/32”, and 8/32” steel 
spheres. It is observed that the band separation 
increases with the size of the sphere. 

Since these secondary waves have centers 
which advance in the direction of the sphere, as 
consecutive bands are considered, it is clear that 
these waves originate with the moving sphere. 
This is substantiated by the observation that 
the waves become further apart as the sphere is 
slowed down. These secondary waves are pro- 
duced by vibrations of the sphere; with each 
expansion of the sphere, an impulse is imparted 
to the water producing a region of increased 
density which travels away from the sphere as a 
wave. In order to calculate the time interval at 
which these waves are produced it is only 
necessary to know the band separation, S, the 
impact velocity of the sphere, Vo, and the wave 
velocity, C. The time interval, 7, is then given 
by T=5S/(C— Vo), where all measurements refer 
to those made in the line of flight of the sphere. 
The time interval between the production of 
successive sharp bands is recorded in the Table I. 
These intervals are proportional to the diameter 
of the spheres, showing that the band pattern is 
produced by oscillations of the sphere. By 
dividing the time interval between the sharp 
bands by the total number of bands enclosed by 
them, namely, three, one obtains the period of 
the vibrating sphere. The period varies from 
0.97 to 1.72 microseconds. 

The period of a vibrating sphere is generally 
expressed in terms of the time required for an 
elastic wave to traverse the diameter of the 
sphere twice. The period is given by T=2kD/W 
where D is the diameter of the sphere, W the 
velocity of the elastic wave, and k is a constant. 


TaBLeE I. The band separation, S, of the vibration 
waves are given for steel spheres having three different 
diameters, D. The period of vibration for these spheres 
is also recorded and compared with Lamb’s theoretical 
values for distortional vibrations. 











S Period Period 
Band of Lamb's 
separation vibration theory 
D cm microsecond microsecond 
1/8” 0.21 0.97 +0.08 1.18 
3/16” 0.32 1.55+0.05 1.76 
1/4” 0.52 





1.72+0.03 2.37 
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Fic. 5. Spark shadowgram of two 4/32” steel spheres 
shot into water simultaneously. The vibration waves 
produce an interference pattern; the arrows indicate the 
paths of reenforcement. The wave has advanced into the 
water 12.2 cm. 


The vibrations of a sphere which give periods 
closest to those measured in this investigation 
are those in which the sphere changes shape, 
altering from an oblate to a prolate spheroid. 
For such vibrations Lamb® has shown that 
k=0.595. Using this constant and a distortion 
wave velocity of 3220 meters/sec., the periods of 
vibration for the three different sized spheres 
were calculated and recorded in the table. They 
agree to within 19 percent of the measured 
values. 

The principal source of error here is not that 
of measurement but of the identification of the 
bands. The bands are not always clear in the 
photographs and further resolution may possibly 
show that the bands are more complex than the 
present experiments indicate. It is clear, however, 
that their origin is the vibrating sphere and that 
the spacing of the bands is proportional to their 
diameter. 

A confirmatory piece of evidence regarding 
the nature of the vibration waves appears when 
two spheres are shot into the water simultane- 
ously. Figure 5 shows two 1/8” spheres moving 
side by side in the water. Since both spheres are 


9A. E. H. Love, The Mathematical Theory of Elasticity 
(Cambridge University Press, New York, 1934). 
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producing waves which overlap, an interference 
pattern is set up. A region of reinforced waves 
and weakened waves is to be seen along radial 
paths spreading out from the spheres. These 
paths are indicated by arrows in the figure. The 
interference pattern is similar to the one observed 
with ripples when two oscillating prongs are 
moved slowly across the surface of the water. 
Measurement of the wave-length which produces 
this interference pattern, in which the distance 








Fic. 6. Comparison of cavity shadow as produced with 
spark (left) and with x-rays (right). Spheres have advanced 
about 3.5 cm into the water, having had an impact velocity 
of 2700 ft./sec. 


apart of the two spheres and the direction of 
reinforcement is used, shows that it is com- 
parable with the spacing of the vibration waves. 
For example, from the interference bands in 
Fig. 5 one calculates that the vibration wave 
separation is 0.14 cm. This is to be compared 
with the average spacing of 0.21 cm, found with 
a 4/32” sphere. More precise determinat:on is 
not feasible with the present shadowgrams. 


4. SHADOW OF THE CAVITY 


As the spherical missile moves through the 
water it tosses the water to one side, and a large 
cavity is produced. The pressure in the cavity 
is considerably less than atmospheric. For a 1/8” 
sphere traveling with a velocity of 3000 ft./sec. 
the cavity can have a volume which is 25,000 
times that of the sphere. The cavity is nearly 
conical in shape and is opaque to the rays from 
the spark. The pressure is very great immediately 
ahead of the sphere, usually being measured in 
thousands of atmospheres, and the refractive 
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index of the water in this vicinity is increased. 
Because of these changes in refractive index the 
tip of the cavity produces an image on the 
shadowgram which is quite different from the 
true geometrical shadow. The tip in the shadow- 
gram has a cusp, and it also flares out at the 
sides, making the shadow much wider than the 
cavity itself. 

The true shape of the cavity was disclosed by 
taking x-ray shadowgrams also of microsecond 
duration. Since x-rays are not refracted, the true 
geometrical shadow was observed. Comparison 
with the spark shadowgrams revealed that their 
shadows of the cavity were as much as 250 
percent wider at the tip than those produced by 
x-rays. For purposes of comparison, an x-ray 
and a spark shadowgram of a 6/32” sphere 
taken under almost identical conditions are 
reproduced in Fig. 6. It is clear that the center 
of the sphere is somewhere near the top of the 
cusp and the most advanced point of the sphere 
is level with the most advanced point of the 
cavity. Hence the light forming the cusp has 
been bent around and behind the sphere. The 








Fic. 7. A shadow of an aluminum slug (with spike), 
taken with the slug at-rest, is shown superimposed on the 
spark shadowgram of the slug as it is moving through the 
water at a high speed; the picture shows the true relation- 
ship between the slug and the cavity in the shadowgram. 
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Fic. 8. Shadowgrams of the tip of the cavity as produced by various missiles. (The magnification can 
be determined from (G) where the cavity width is 1.0 cm.) (A) A slow 4/32” steel sphere at 915 ft./sec. 
(B) A 6/32” steel sphere at 1950 ft./sec. (C) A 8/32” steel sphere at 1710 ft./sec. (D) A 1/16” cylindrical 
dart 19/32” long and pointed at the ends; velocity between 1000 and 3000 ft./sec. Long axis is horizontal 
and perpendicular to the plane of the paper. (E) A cylindrical dart like the one in (D) with the long 
axis inclined about 30° with horizontal and perpendicular to the plane of the paper, velocity between 
500 and 2000 ft./sec. (F) A 1/16” cylinder, 17/32’’ long with square ends. Velocity between 500 and 
2000 ft./sec.; long axis in plane of paper. (G) A 1/16” cylinder, 17/32” long with square ends. Velocity 
between 500 and 2000 ft./sec. The cylinder is toppling. (H) A spin stabilized slug of aluminum with a 
nose having a total cone angle equal to 90°. The velocity is between 1300 and 1800 ft./sec. Total length 


of slug is 1.5 cm. (I) A flat-nosed slug. 


position of the sphere relative to the shadowgram 
of the cavity was also verified by another method 
in which aluminum slugs were used. The alumi- 
num slugs were machined with hemispherical 
ogives to which a needle point was attached. 
The tip of the needle showed little or no refrac- 
tive effects and the hemispherical nose of the 
slug was located with respect to the tip of the 
needle by direct measurement. A spark shadow- 
graph of the slug at rest was superimposed on 
one of the slugs traveling through the water 
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with a velocity of 1500 ft./sec. and the composite 
photograph is shown in Fig. 7. The shadow of 
the cavity is about 38 percent wider than the 
cavity itself at the equatorial plane of the 
spherical ogive. The light forming the cusp bends 
in behind the spherical ogive at all points along 
a 50 or 60 degree arc. The length of this arc 
increases with the pressure in front of the sphere. 
The remainder of the semi-circular arc is near a 
low pressure region. In the low pressure region 
the rays are bent away from the sphere while in 
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Fic. 9. The distortion in the grid shows the magnitude 
and direction of the pressure gardient in the vicinity of a 
fast moving sphere. The 6/32” steel sphere was moving 
with a velocity of about 1900 ft./sec. The lines in the grid 
are 1 mm apart. 


the high pressure region they are bent toward 
the sphere. 

An interesting confirmation of this explanation 
for the large shadow of the cavity tip occurs 
when two spheres are shot into the water simul- 
taneously. The flow velocity between the spheres 
is increased by the proximity of the two spheres, 
resulting in an even greater reduction in pressure 
on one side of each sphere. This is shown in 
Fig. 5 where the shadowgram shows a broadening 
of that lobe of the cavity tip which is closest to 
the neighboring sphere. 

The pressures acting on the spheres are quite 
high. The force in dynes which acts to retard 
the sphere equals 0.15AV? where V is 
velocity in cm/sec., and A the projection area of 
the sphere in cm?*. If one assumes that this force 
is spread out over the entire projection area of 
the sphere, an assumption which gives the 
largest possible area, then the average pressure 
in dynes/cm? is equal to 0.15V?. For a sphere 
moving with a velocity of 4800 ft./sec. the 
average pressure is 3290 atmospheres, and a 
‘sphere having a velocity of 1200 ft./sec. has an 
average pressure of 206 atmospheres. The pres- 
sures are large enough to produce considerable 
deviation of the rays; the pressure of 3290 
atmospheres will produce a change in the index 
of refraction equal to 3.3 percent. 
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The cusp at the tip of the cavity is more 
pronounced when fast spheres are used. In Fig. 8 
are shown several of these cusps. The one for 
the 8/32” sphere, Fig. 8C, shows that the light 
has been so refracted as to focus into a straight 
vertical line behind the sphere. A slow 4/32” 
sphere traveling at 915 ft./sec. produces a pin- 
point of light instead of a cusp, as shown in 
Fig. 8A. The interesting cavities shown in Fig. 
8D and Fig. 8E are of steel darts viewed endwise. 





Fic. 10. A double spark shadowgraph, taken before and 
during the flight of a sphere through water, is shown when 
a one-millimeter point-grid is used. The 6/32” steel sphere 
was moving with a velocity of about 2000 ft./sec. The 
white streak to the right of the sphere is a blemish. 


In one of these the refraction effects cause the 
tip to be wider than any other part of the cavity. 
A cylinder traveling sidewards produces a shadow 
with bulbs at the end of the cylinder, as shown 
in Fig. 8F. The cylinder of Fig. 8G is toppling 
in a plane which makes an angle with the rays 
from the spark. An aluminum slug with a 90° 
conical nose produces the symmetrical shadow 
of Fig. 8H, while a cylindrical slug with a flat 
end gives the shadow shown in Fig. 81. By 
inspecting these shadowgrams one can determine 
in a general way where the points of greatest 
pressure are located. 
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5. PRESSURES BETWEEN THE CAVITY 
AND IMPACT WAVE 


The spark shadowgraph can be used to meas- 
ure pressures in other regions as well as those 
appearing in the shock wave. For this purpose a 
grid was placed between the tank and the spark. 
The grid was made by scratching lines one 
millimeter apart on a sheet of Plexiglas; it was 
placed on the outside wall of the tank. The 
pattern of the grid around the missile is distorted 
by the changes in pressure as is shown in the 
shadowgram of Fig. 9. The rays above a line 
running perpendicular to the axis of the cavity 
and tangent to the two lobes of the cusp are 
bent downward, while those below are bent 
upward. The distorted pattern is also symmetri- 
cal about the axis of the cavity. A disturbance 
in the pattern is observed to extend to at least 
ten sphere diameters in front of and to the side 
of the missile. 

_ When a grid made up of a series of opaque 
dots is used, a better idea of the distribution of 
pressure is obtained. This grid was made by 


Fic. 11. Wave 
originating at the 
surface is shown 
being reflected from 
the steel plates on 
both sides. The ar- 
row points to what 
appears to be a 
Mach reflection 
The wave is 13.5 
cm below the sur- 
face. 
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Fic. 12. A concave surface of brass mirror (dotted line 
shows surface in profile) is shown focusing the wave along 
R. Wave P originates at O, center of mirror at C. Concave 


surface of radius 8 cm was cut out of 1’°*5’’X*5” block. 
Transmitted wave at T. 


indenting a Plexiglas plate with a needle, the 
indentations being aligned in rows and columns 
one millimeter apart. The shadow cast by the 
indentation consists of an opaque part and a 
bright speck of concentrated light; the latter is 
black in the print. The result of recording a 
before- and during-the-shot shadowgraph on the 
same plate is shown in Fig. 10. The double 
image of the black dots very clearly shows the 
direction of the displaced ray and the orientation 
of the lines of the pressure gradient. Because the 
displacement of the rays occurs in planes which 
are perpendicular to the rays and which are in 
front of and behind that of the bullet, the 
pressure gradient map may be regarded as a sort 
of average map for all these planes. The gradient 
lines in front of the sphere appear to be nearly 
straight and to spread out radially from the 
nose of the spherical missile. Those behind the 
missile are curved, the curvature being greatest 
near the sphere. From measurements on the 
displacement of the rays it is possible to obtain 
values for the large pressures near the sphere. 
As yet these calculations have not been made 
because of their length; the distorted grid pat- 
terns have only been used in a qualitative way 
to show the direction of the pressure gradient 
and to indicate its magnitude. 
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Fic. 13. Concave brass surface with center at C and 
radius of 8 cm is shown focusing wave at Y. Vibration 
waves are shown converging and diverging at the focus. 
Block is 155 inches. 


6. REFLECTION OF WAVES 
Mach Reflection 


Shock waves traveling through a compressible 
medium such as air are known to be reflected in 
a manner quite different from any other waves. 
The incident and reflected waves do not always 
meet at the reflecting surface but at a point 
just above it. In a spark shadowgram this 
configuration looks like a Y in which the stem is 
the vertical section and the incident and reflected 
waves are the arms of the Y. These reflections 
are called Mach reflections.'° 

Water is less compressible than air and the 
Mach reflections are less easily produced. Mach 
intersections which are similar to the Mach 
reflections have been observed with waves that 
were generated by underwater explosions.® Mach 
reflections occur when the wave is either very 
strong or when the wave is reflected near 
glancing incidence. Since impact waves are much 
weaker than explosive waves, Mach reflections, 
if they are to appear at all, would appear at 
very small angles of incidence. In Fig. 11 is 
shown a wave being reflected from a steel plate 
at a glancing angle of 9.5°. The pressure in the 
wave is about 40 atmospheres. A stem, perpen- 
dicular to the plate, is visible although it is not 
nearly as well defined as those observed in air. 
The arrow in the picture points to the stem. 
For a good example of the Mach Y’s in air, one 


10 E. Mach, Sitz. Akad. Wiss. Wien 77, 819 (1878). 
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is referred to Fig. 2 and Fig. 4, where the wave 
from the splash and the head wave of the missile 
intersect. 


Reflection from a Spherical Surface 


When the acoustic waves in water strike a 
concave surface, the curvature of the wave is 
inverted in the same way that it is for light 
waves and sound waves in air, thus converting 
a diverging wave into a converging one. In Fig. 
12 is shown a diverging wave front being reflected 
by a brass mirror having-a radius of curvature 
equal to 8 cm. The dotted lines in the picture 
show the profile of the mirror; the origin of the 
wave is at O and the center of curvature of the 
mirror is at C. The envelope of the reflected 
wave now converging is labeled R and the 
diverging incident wave P. A shadowgram of 
another shot triggered at a later time is shown 
in Fig. 13 when the reflected wave has passed 
through the focal point and again become di- 
verging. The vibration waves are shown ap- 
proaching and receding from the focal point, Y; 
they are converging below the point and diverg- 
ing above. The waves are reflected by the brass 
block in the same way that light waves are 
reflected by mirrors in a reflecting telescope and 
other optical equipment. 

The reflection from a convex surface is demon- 
strated in Fig. 14 where a lenticular-shaped 
object of cork composition reflects the wave. 
The convex surface has a radius of 11 cm and 
the wave originates_at the water surface 11.5 cm 
away. The reflected wave has its center about 
3.8 cm below the cork reflector; it is between the 
principle focus of the mirror and the mirror 
surface, just like optical images in convex 
mirrors. , 


7. REVERSAL 
Reversal at Normal Incidence 


During a study of reflected waves it was 
observed that some reflected waves appeared 
with an additional presence-of-light band ahead 
of the wave. Reflections from air, paraffin, and 
cork, for example, showed the reflected wave 
with the presence-of-light band foremost, where- 
as the incident wave has an absence-of-light band 
in the lead. This addition to the band pattern is 
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shown in Fig. 15 where a wave reflected from a 
block of paraffin is shadowgraphed. The wave 
was formed at the water surface, 10 cm above 
the block, and has traveled about 4 cm after 
reflection. The incident wave can be seen at the 
side. The missile has not yet hit the paraffin 
although it has perforated the reflected wave. 
This new light band appears to be about as wide 
as the light band which regularly follows the 
black band. However, in most of the reflections 
from other media which also produce a leading 
light band, the bands are broad and indistinct. 
This is shown in Fig. 16 where a wave is reflected 
from a block of cork. (The cork was large grained 
but had been smoothed by a sander.) Soft wood 
also produces this type of reflected wave. Re- 
versal by cork-composition is seen in Fig. 14 
where the wave is reflected from a convex 
surface. Similarly when a wave is reflected from 
an air surface, light is thrown out ahead of the 
wave. In Fig. 17 the wave is reflected from the 
end of an air-filled cylinder having a thin rubber 
membrane for a cover. The membrane is some- 
what obscured in the picture by the cavitation 
which follows the reflected wave. 

Other materials such as steel, brass, Plexiglas, 





Fic. 14. Reflection from convex surface of cork showing 
reversal. Virtual origin of reflected wave is at O’ which 
lies between the mirror’s principle focus and its surface. 
Cork surface is 11.5 cm from the water surface. Note 
waves originating from internal reflections within cork lens. 
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Fic. 15. Reflection from paraffin block showing reversal. 
Wave is 4.1 cm above block at furthest point. 





Fic. 16. Reversal upon reflection from block of large- 
grained cork. Wave is 1.4 cm above block at furthest 
point. 





Fic. 17. Reversal of reflected wave from container of air 
covered with thin rubber diaphragm. Cavitation in water 
above rubber diaphragm causes surface to look fuzzy. 
Diameter of container is 12 cm, rubber is 0.028 cm thick. 
Wave is 3.5 cm above diaphragm at furthest point. 


and lead reflect waves without any change in the 
arrangement of the dark and light bands. Figure 
18 shows a wave being reflected from a 1/4” 
steel plate and Fig. 19 shows a wave reflected 
from a 1/2” sheet of Plexiglas. Both waves 
resemble the incident wave in regard to the 
order of bands. In the last shadowgraph there 
is also present. a wave reflected at the lower 
Plexiglas surface. Since the wave in this case is 
being reflected by water which is less dense than 
Plexiglas, the wave is reversed and possesses a 
light band in front. Figure 20 shows a reflection 
from a brass plate one-half of which has been 
roughened, demonstrating that the smoothness 
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Fic. 18. Wave is shown being reflected from 1/4” steel 
plate. Wave has advanced 3.1 cm after reflection. 


of the reflection surface does not affect the order 
of bands upon reflection. 

Reversal of the band order takes place when 
the wave is reflected from a medium of lower 
acoustical impedance, where acoustical imped- 
ance is defined as the product of density and 
wave velocity. The acoustic impedance for the 
materials encountered in this report is recorded 
in Table II. All substances below water in the 
table have low impedances and reflect water 
waves with a light band in front. Likewise, 
waves reflected from inside those substances 
which are listed above water in the table also 
are reflected at the substance-water interface 
with a light band foremost. This is in accord 
with the accepted theory of wave reflection, 
namely, that when a reflection occurs in a wave 
going from a medium of high to low acoustical 
impedance there is a change in phase and the 
high pressure wave is reflected as a low pressure 
wave. Just what this change of phase does to 
the steep pressure front of the shock wave is not 
known. It is possible that only a change in the 
distribution of the pressure is accomplished; on 
the other hand a definite negative pressure region 
may develop. In either case the effect on the 
rays from the spark is to throw the beam forward 
producing the reversal in band order. A wave 
possessing a front having a pressure lower than 
that of the medium through which it is traveling 


would refract the rays from the spark so as to 


throw them in advance of the wave; but it is 
also known that there are certain distributions 
6f positive pressure behind the front which will 
also refract rays so that they produce light in 
advance of the dark band. It is interesting to 
note that paraffin which possesses an acoustic 
impedance that is only slightly less than that of 
water produces a very narrow light band, while 
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Fic. 19. Wave R is reflected from 1/2” Plexiglas plate 
without reversal. However, second wave SR is reflected 
from lower surface of Plexiglas-and is reversed. Wave 
had advanced 4.8 cm after reflection. 





Fic. 20. Wave reflected from partially roughened brass 
plate showing no reversal. Left half was rough, right half 
smooth. Reflected wave is 5.9 cm from plate at furthest 
point. 


other substances having a very small acoustic 
impedance produce wide bands. 


Reversal at Glancing Incidence 


When the wave front is nearly normal to the 
reflecting surface (the rays are then at glancing 
incidence), it is found that reversal occurs even 
when the reflecting surface has a high acoustic 
impedance. This is illustrated in Fig. 21 where 
a wave is shown being reflected from a 1/2” 
plate of steel. The wave front formed by rays 
having glancing angles of about 45° or less 
consists of a band of light at the front, ahead of 
the wave followed by a dark band. This reversal 
does not depend critically on the angle and 


TABLE ITI. Acoustic impedance X 10~* g cm sec. 





Steel 4.67 
Brass 3.60 
Plexiglas 0.331 
Water 0.148 
Paraffin 0.126 
Cork 0.099 
Rubber 0.003 
Air 0.000043 
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Fic. 21. Wave showing reversal upon reflection at 
glancing incidence from 1/2” steel plate. S is the surface 
of water, R the reflected wave, and P the incident wave 
of radius 13.4 cm. 





Fic. 22. Wave is shown reversed when reflected at 
air-water interface. S is the water surface, P the incident 
wave of radius of curvature 25.6 cm, and R the reflected 
wave. Angles of reflected rays measured with respect_to 
the normal are also shown. 


appears under various orientations of the re- 
flecting plate. A reflection from an air surface 
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Fic. 23. Striations are shown when wave is reflected 
from steel plate. P is the incident wave and R the reflected. 
This is a detailed view of Fig. 18. 


also shows reversal as can be seen in Fig. 22. It 
was once believed that the reversal could be 
attributed to the fact that the reflected wave 
was close behind the leading wave. When, 
however, an identical arrangement of wave fronts 
was produced by shooting two spheres into the 
water simultaneously, no reversal was observed; 
see Fig. 5. 


8. STRIATIONS 


A very interesting set of striations put in their 
appearance when the waves from vibrating 


spheres are reflected. The striations do not 
appear to move, but form a stationary pattern 
which persists for a very short interval. They 
are from 0.02 to 0.08 cm apart. Figure 23 shows 
these striations when the wave is reflected from 
1/4” steel plate. The incident wave P with its 
vibration bands is also shown. The striae are 
parallel to the reflecting plate and extend 
throughout the region bounded by the plate and 
the reflected wave R. The incident and reflected 
rays are nearly normal to the reflecting surface 
in Fig. 23; when the rays are near glancing 
incidence, the striations produce a different 
pattern, as shown in Fig. 24. Here the striae 
follow curved paths behind the reflected wave 
front and are perpendicular to both the wave 
front and the reflecting steel plate. There are 
three separate waves present in this picture, for, 
in addition to the primary P and reflected wave 
R there is also the head wave H produced by the 
Rayleigh wave on the reflecting steel plate. In 
the region bounded by the head wave, reflected 
wave, and primary wave, the striations form a 
pattern of parallel straight lines. 
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» REFLECTING PLATE. 





Fic. 24. Striations are shown when wave is reflected at glancing incidence from steel plate. R is re- 
flected wave of radius 9.8 cm, P the incident, and H the head wave from the steel plate. 


Another example of these striations is shown 
in Fig. 25 where a wave RR is shown being re- 
flected from two lead bars, L and LL, 1/8” 1”. 
The bars are separated by about an inch and the 
diffracted waves D are also visible. The region 
traversed by the reflected wave is filled with 
striae. (The other objects which cast white 
shadows are supports that are not in the plane 
of the shot.) 

Since the striations always appear when a 
reflected wave passes through a region just 
traversed by the incident wave, it is believed 
that they are produced: by the secondary or 
vibration waves from the vibrating missile. 
These waves always occupy the region immedi- 
ately behind the leading shock wave and inter- 
sect each other when the wave is reflected. In 
support of this idea one finds that the distance 
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between striae increases with the spacing be- 
tween the vibration waves; the striae are farther 





Fic. 25. Striations are shown when wave is reflected 
from lead bars LZ and LL. D is the diffracted wave and 
R the reflected. Note the vertical bands where the reflected 
waves intersect. Lead bars are 1/8” thick. 
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apart when produced by large rather than small 
spheres. Since these striae are oriented in a 
direction which is normal to a line which bisects 
the incident and reflected rays, it is quite 
possible that these striae are shadows of the 
wake produced by the intersection of these 
vibration waves. This wake may be either a 
region of abrupt temperature change or simply 
one of great turbulance. It is interesting to note 
that the striae occupy the same relative position 
as the “‘slip stream’’ does when shock waves in 
air intersect. 

In addition to the striations in Fig. 25 a set 
of bands appears, where the diffracted wave 
from the two lead plates, L and LL, intersect. 
These bands resemble the optical diffraction 


pattern produced by a narrow slit. While the 
pattern resembles an optical diffraction pattern, 
the pattern is more likely to be the shadow of 
an interference between the acoustic waves 
themselves. It is quite plausible that the vibra- 
tion waves from the sphere, after being diffracted 
from the lead bars, intersect and produce these 
bands. The details of this wave interaction and 
the method by which the rays from the spark 
produce this effect are still obscure. 
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Lightning Recorded by Aeroplane Propeller 


J. M. Bryant AND M. NEWMAN 
Department of Electrical Engineering, University of Minnesota, Minneapolis, Minnesota 


(Received March 11, 1946) 


A lightning stroke to the propeller blades of an aeroplane in flight left a remarkably complete 
record of its own characteristics. Sputtering of the discharge on the rotating blades left a clear 
time scale of duration of the stroke. The lightning current direction was also recorded by mag- 
netization of steel rings inherent in the propeller mechanism. 


RECENT lightning stroke to an aeroplane 

in which the rotating propeller intercepted 

the discharge gives a measurement of the time 
duration of lightning at the cloud level that may 
be met in aircraft travel. The aeroplane was an 
all-metal twin-motored transport traveling at a 
normal cruising speed of 190 miles per hour at 
an altitude of 11,000 ft. A retouched photo- 
graphic representation of a probable discharge 
path involving the aeroplane is given in Fig. 1. 
The propeller tip protrudes furthest below the 
aeroplane fuselage, thus forming the highest 
gradient point and hence the most probable dis- 
charge terminus of the plane when in a strong 
vertical electrostatic field in or between clouds. 
The rotation of the propeller at 1200 r.p.m. 
spread out the pitting effects of the lightning 
discharge to a time scale, as indicated in Fig. 2, 
determined by the rate of rotation of the pro- 


VOLUME 17, JULY, 1946 


peller, and the distance of the track from the 
axis of rotation. 

Upon examination of the propeller, shown in 
Fig. 1, it was found that all three blades showed 
burning marks owing to relatively low current 
values. Instead of a single stroke of short dura- 
tion (40 to 100 usec.) usually associated with 
lightning discharges, the stroke lasted for at 
least 0.15 second. This is proven by the fact that 
each blade cut across the ionized path which 
shifted sufficiently to leave separate sweep traces 
so that blade No. 1 as shown in Fig. 2 was in 
contact with the arc at least four successive times. 

The polarity of the propeller was positive to a 
negative cloud region as determined by the mag- 
netic polarity of steel retaining rings about the 
base of the blades as shown at R in Fig. 3. The 
propeller blades and the housing are of non-mag- 
netic material, and the steel rings were insulated 
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Fic. 1. Photograph showing pit marks left by a multiple 
lightning discharge as the rotating propeller blades swept 
through the stroke channel. The pit marks were retouched 
for clarity. The lightning channel was drawn in a probable 
position at the completion of the discharge process, the 
channel being drawn out circumferentially by the last 
receding propeller blade contact point. 


from the blades by plastic material bushings 
inherent in the mechanism designed to change 
the propeller pitch while in flight. Currents in 
each blade were thus confined to flow entirely 
inside each ring. Each ring is split for assembly 
purposes and each of the three rings was found 
to be magnetized in the same circular direction, 
showing from the polarity that the current flowed 
through the engine shaft to each blade and thence 
to the cloud. 

Each complete revolution of the propeller repre- 
‘sents 0.05 second and the time scale between 
centers of blades represent 1/60 second. Rows 
of pit marks show in the forward face of each pro- 
peller blade and a more serious burning occurred 
. on the trailing edges of the blades. One of the 
pit marks is shown much enlarged in Fig. 4, 
the depth of the pit being about one-tenth mm, 
the melted material being thrown outward by 
centrifugal force toward the tip of the blade. 
The whole series of events points to the prob- 
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Fic. 2. Photograph showing the pitting which could be 
caused by successive current contacting the propeller 
blade as it swept through the discharge channel. 


ability that the plane was hit by a “multiple 


‘stroke” of lightning of relatively long duration. 


The nature of the discharge from the cloud en- 
tering the blades was such as not to burn the 
blade between pits, which could be accounted 
for by sudden peaks of discharge current occur- 
ring at short intervals, and thus damaging the 
blade in “‘pits’’ along the line of travel of the 





Fic. 3. Photograph showing location of the steel rings, 
R, and the insulating sleeves, S, in the propeller pitch 
control mechanism. 


JOURNAL OF APPLIED PHYSICS 














arc channel. The contact time with the trailing 
edges was so much longer, in proportion to the 
ratio of the space between blades to that between 
pits, as to permit considerable burning. 

Multiple strokes of lightning to the earth's 
surface have been previously recorded by photo- 
graphic methods,! by magnetization effects in 
the ‘‘fulchronograph,’”? and by special magnetic 
oscillographs. An oscillographic record is re- 
produced in Fig. 5 of a long duration discharge 
with repeated peaks, as recorded by McEachron.* 
However, the time scale recording method fur- 
nished by the rotating blades gives a spacing of 
the successive peaks in the case of this airplane 
stroke as of much shorter intervals, as seen in 
Fig. 2, these intervals being from 10 to 100 
usec. between the pits across the face of each 
blade. 

It is of interest that the aeroplane propeller 
left a record of effects of the lightning channel 
near the point of origin of the stroke. From its 
character of shorter intervals and absence of very 
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Fic. 4. Photographic enlargement of a typical pit mark, 
showing how the molten metal was thrown out radially 
by centrifugal force because of propeller rotation. 


1B. F. J. Schonland and H. Collens, Proc. Roy. Soc. 
A152, 595 (1935). 


2C. F. Wagner and G. D. McCann, Trans. A.I.E.E. 
60, 497 (October, 1941). 


3K. B. McEachron, Elec. World 113, No. 6, 430 (Febru- 
ary 10, 1940). 
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Fic. 5. Reproduction of a crater lamp oscillogram, from 
a paper by K. B. McEachron, showing successive peaks 
and continuing current in a stroke to the Empire State 
Building. 


severe currents, it is quite likely that the stroke 
was from cloud to cloud rather than from cloud 
to ground, in which case it would be the first 
recorded data on such a stroke. A short cloud to 
cloud discharge might account for the shorter 
time between the consequent peaks than in 
those cases of strokes to ground recorded in 
Fig. 5. Another more likely possibility based on 
some preliminary calculations is that the aero- 
plane intercepted a discharge from a charge 
center, progressing upward toward the positive 
charges in the higher regions of the same cloud 
rather than to another charged cloud or to earth. 
Such a theory would fit the observed polarities 
and most probable path of the discharge, and 
in such a case the current peaks could correspond 
to sudden jumps of current in the progressive 
development of the lightning channel such as 
those observed by Schonland,‘ for the formative 
stages of stroke to earth. The observed time in- 
tervals of the pitting on the propeller correspond 
approximately to the time intervals between the 
light-intensity jumps observed by Schonland. 
The fact that the over-all time corresponds to 
about ten times greater than the duration of the 
formative stage would fit with a case of a very 
long developing discharge directed through the 
aeroplane upwards from the cloud charge con- 
centration rather than down to earth. 

In conclusion the writers wish to express their 
appreciation of the cooperation of the Mid- 
Continent Airlines and the Northwest Airlines 
in the investigation of lightning strokes to 
aircraft. 


4B. F.J. Schonland, The Lightning Discharge (Clarendon 
Press, Oxford, England, 1938), p. 6. 
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A Microwave Dielectric Loss Measuring Technique 


WILLIAM R. MAcLEANn 
Electrical Engineering Department, Polytechnic Institute of Brooklyn, Brooklyn, New York 


(Received April 8, 1946) 


Dielectric loss is measured by Q determinations of a resonator partially filled with the 
sample. Double sample technique is used to eliminate dominant spurious losses. Metal loss is 
reduced by confining field mostly to sample. Detuning for Q determination is accomplished by 
large movement of small rod. Detuning is calculated by action theorem. 


INTRODUCTION 
M EASUREMENT of the dielectric constant, 


e, and the loss factor, d, of dielectrics at 
ultra-high frequencies has been described by 
several authors. For instance, Works, Dakin, 
and Boggs' have devised a method using a 
tunable reentrant resonator, and Englund? has 
used a transmission line method involving a 
plug made of the dielectric to be measured. All 
these methods and also the one described in this 
paper are based upon measurements of the Q of 
a resonator or other tuned circuit. 

When the loss factor is very small certain 
difficulties arise which might be tabulated ap- 
proximately as follows: 


(a) Resistance losses at joints, particularly at sliding 
joints, can cause trouble. 

(b) It often happens that the loss due to the dielectric is 
small compared with the loss due to the metal and, hence, 
one is in the unfavorable position of measuring a small 
quantity in the presence of a larger parasitic one. 

(c) The loss due to the metal cannot be calculated due 
to the variability of the “roughness factor.”’ 

(d) A measurement of the loss due to the metal is often 
placed in some doubt if it is measured at a different fre- 
quency or if it is measured at a different internal con- 
figuration (this problem arises since the resonator without 
the sample does not resonate at the same frequency as with 
the sample). 

(e) There is a certain amount of loss due to the coupling 
into and out of the resonator for which a correction might 
be necessary. ‘ 

(f) Since the Q of the resonator is high in the case of 
small loss factors, the detuning to half-power points 
involves an extremely small shift in either the resonator or 
the oscillator; for instance, in the transmission line method 
it is necessary to measure a movement of only a few 
thousandths of an inch of the sliding plug. 

' Works, Dakin, and Boggs, A.I.E.E. Tech. Pap. 44-161, 
May, 1944. 


*C. R. Englund, “Dielectric constants and power factors 
at centimeter wave-lengths,”’ Bell Sys. Tech. J. 23, 114- 
129 (1944). 
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(g) The half-power point of. the crystal itself is at times 
in doubt, particularly if an assumption of a square law is 
made or if the crystal is calibrated at a much lower fre- 
quency. 

The present method is restricted to the diffi- 
cult measurement of small loss factors. If the 
loss factor is large, it will not work. It requires 
also a preliminary approximate determination of 
the dielectric constant, e. The final measurement 
by the present method checks the value of ¢ 
and determines a small loss factor, d, by a method 
which minimizes most of the difficulties tabulated 
above. The means by which these are overcome 
are briefly the following : 


(a) By the use of a cylindrical resonator in a transverse 
electric mode no current flows across the cracks. 

(b) The sample is placed inside the resonator as a 
dielectric core which confines the field almost entirely to 
the dielectric. This practically eliminates the losses on the 
side walls. The losses on the end pieces can in principle be 
made relatively as small as desired by using a resonator 
and sample which is long enough. 

(c) A calculated value for the loss on the side walls 
although inaccurate is sufficient to prove it negligible. 

(d) The actual losses on the end pieces, whatever they 
nay be, are eliminated in the measurement by observations 
taken with first one and then two identical samples. 





Fic. 1. Resonator disassembled showing samples. 
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(e) Any loss due to the input and output coupling is also 
eliminated at the same time. 

(f) The resonator is detuned for the Q measurements by 
pushing a small metal rod through the center. The detuning 
created by the rod is so small that a large movement may 
be obtained no matter how high the Q. 

(g) The method of measurement includes a determina- 
tion of the half-power point of the crystals at the frequency 
of operation. 


DESCRIPTION OF APPARATUS 


It is perhaps best to describe first the physical 
apparatus used in an experimental test of this 
method involving the measurement of poly- 
styrene at 10-cm wave-length. Figures 1—3 are 
photographs of the equipment. In Figs. 1and 2can 
be seen a hollow metal cylinder and two identical 
dielectric samples and also a cover plate. Into 
the metal cylinder protrude two small rotatable 
and retractable coupling loops, one for the driv- 
ing oscillator and one for the detecting crystal. 
To assemble the equipment one sample, for in- 
stance, is placed inside and the cover inserted. 
The cover is a loose sliding fit in the metal cylin- 
der. The cover is then clamped in place with a 





Fic. 2. Resonator disassembled showing samples. 


clamping spider which is also visible in Figs. 1 
and 2. As can be seen, the samples are bored 
through the center with a small hole through 
which moves the small metal rod also visible in 
the figures. Attached to the cover plate is a 
millimeter scale for measuring the penetration 
of the tuning rod into the sample. The equipment 
is assembled with either one or both samples 
inside. In Fig. 3 the equipment is shown fully 
assembled and ready for measurement with two 
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Fic. 3. Resonator assembled. 


samples inside. It is observed, of course, that the 
dielectric sample is considerably smaller than 
the diameter of the metal cylinder. It is an essen- 
tial part of the scheme of measurement that this 
dielectric core confine the field to keep it away 
from the side walls. If the two samples are alike, 
the resonator will have a resonant frequency 
with two samples identical to that with one 
sample. 


METHOD OF MEASUREMENT 


The method of measurement is basically a Q 
measurement as is common to all these methods. 

Assume the resonator assembled with one 
sample and assume a signal transmitted through 
the resonator by means of the two loops. To 
achieve this condition it will, of course, be neces- 
sary to tune the oscillator very close to the 
resonant frequency. It will then be found that 
by means of the tuning rod the resonator can 
be tuned to resonance with high precision. The 
loops are, of course, retracted until the coupling 
is the smallest value for satisfactory response. 
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Since in the mode of oscillation used the mag- 
netic lines are vertical, the response of the crystal 
will vary as the loop is rotated in an approxi- 
mately sinusoidal manner. Actually, the open 
circuit voltage generated in the crystal loop will 
vary exactly in a sinusoidal manner but the 
rotation will have a very slight effect on the 
impedance seen by the crystal. The weaker the 
coupling the less is this effect. If it were not for 
this impedance variation the crystal could be 
calibrated to its half-power point with certainty 
by rotating the loop 45 degrees from the position 
of maximum response. In actual operation with 
a small value of loss factor and a fairly sensitive 
crystal galvanometer and a driving power of the 
order of a watt available in the oscillator this 
variation in impedance will be negligible and 
can be so proved by repeating the calibration 
with two samples instead of one in the resonator ; 
for in this case the impedance fluctuation differs 
by a factor of two and hence if both measure- 
ments give the same calibration this phenomenon 
is negligible. The crystal is now calibrated for 
its half-power points. With one sample in the 
resonator it can be tuned to resonance with the 
push rod, and thereafter it can be tuned to the 
upper and lower half-power points by advancing 
and retracting the push rod. Readings of the 
push rod penetration at the upper and lower 
half-power points will indicate the Q of the 
resonator if the characteristics of the push rod 
as a tuning element can be calculated. It will be 
shown later how this is done. If it is impossible 
to reach the half-power points with the push 
rod, then it is too small: on the other hand, if 
the half-power points lie too close together the 
push rod is too large. It will be shown later how 
to estimate the proper size for the push rod al- 
though this is not critical. 

By this method it is therefore possible to de- 
termine the Q of the resonator when it contains 
first one and then two samples. 

It is more convenient to calculate in terms of 
_ damping factors, d, which are the reciprocals of 
the Q’s. Let us suppose, for instance, that the 
stored energy of the resonator is given by W and 
the power loss is given by D. Then, as is well 
known: 


d=D/wW. (1) 
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If there are two sources of loss it is possible to 
consider a damping factor, d, due to each, after 
which the total damping factor is their sum. This 
is clear from (1). 

In actual operation we speak of a d; observed 
with one sample and a dz observed with two 
samples. These are determined experimentally. 
Now suppose that d; consists of several parts 
one of which, do, is due to the dielectric loss. 
We first seek do, i.e., the damping factor due 
to the dielectric when two pieces are used. Sup- 
pose that the same maximum field strength 
existed in both cases, then W and D of (1) are 
both doubled, hence do is the same in both cases. 
Now let us consider dz, i.e., the loss due to the 
end plates when one sample is present, and 
seek dg when two samples are present. In this 
case, W is doubled but D is the same, hence: 


dg?) =4dr. (2) 


Consider now the loss factor d. due to the 
coupling loops and suppose for the moment that 
in introducing a second dielectric sample we also 
introduced another identical pair of coupling 
loops with associated loads. It is obvious that 
then the double resonator is merely two copies 
of the original one and hence d, would not be 
changed. But then imagine one of the pairs of 
coupling loops removed. This would halve the 
loss. Therefore the loss factors d. obey the same 
law as dz. 

The losses on the tuning rod and the losses 
on the side walls do not behave like (2) but it 
turns out that they can be made negligible as 
will be shown later. 

Neglecting then the tuning rod and side wall 
loss we can express d,; and dz in terms of dy and 
d, which latter represents the sum of end plate 
and coupling losses. These would be related by 


d,=dy)+d,, 
daady+-Ade, 3) 


which can immediately be solved as: 
do = 2d2—d. (4) 


We now have by a very simple formula the 
loss due to the dielectric in terms of the two 
measurable loss factors. However dp is the loss 
due to the dielectric but it is not the loss factor of 
the dielectric. These differ because the field is 
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not entirely confined to the dielectric. However 
a simple correction factor can be calculated as 
will be shown later. 


DESIGN 


The problem of designing the equipment will 
be posed as follows. The dielectric constant, e, 
and the vacuum wave-length, A, at which the 
measurement is to be made are given. It is first 
of all necessary to calculate the diameter, 2a, of 
the sample, and its length, /. It is then necessary 
to calculate the diameter, 26, of the resonator. 
This can be so large that it has no appreciable 
effect on the frequency of resonance. That means 
the frequency should be practically the same as 
for a resonator of infinite diameter. This is 
readily achieved since the field is confined almost 
entirely to the dielectric. In addition the di- 
ameter of the resonator must be so large that 
the loss factor due to the side wall is negligible. 
It is then necessary to calculate the approximate 
size of the tuning rod and to calculate its calibra- 
tion curve as a tuning element. Then the factor 
relating the dielectric loss, d, with the dy of (4) 
is calculated. One should then calculate approxi- 
mately the loss due to the tuning rod. 

It is not necessary to construct a new resonator 
for every measurement since the diameter is not 
critical as long as it is big enough and since the 
lid slides inside and adjusts itself to the sample 
length. As soon as these quantities mentioned 
above are known one can immediately proceed 
to the measurements just described. The rest 
of this article will be devoted to the calculation 
in convenient form of the necessary design quan- 
tities. They will be given as functions of ¢ and \ 
so that they can be used directly for this purpose. 
It turns out that there is one arbitrary degree of 
freedom in the evaluation of these quantities 
and the curves given at the end of this article 
are for a particular choice of this parameter, 
which merely determines whether the resonator 
should be short and fat or long and thin. 


SIZE OF SAMPLE 


Since it is planned to construct a resonator so 
that the outer wall has a negligible effect, calcu- 
lations are first made assuming this radius, }, to 
be infinite. That is, the system is taken to be a 
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dielectric wave guide bounded on two ends by 
infinite conducting planes. A discussion of such 
a configuration is found for instance in Schel- 
kunoff.* For the present application the TEo: 
mode is desired. 

In cylindrical spaces it can be shown that the 
complete transverse electric field can be derived 
from a Hertz vector, Q, having only an axial 
component, Q. For this component the desired 
solution of the homogeneous wave equation is 
given by: 


Q=1 sin(yZ)Jo(pr) 
. . « inside dielectric core, (5) 
Q=A sin(yZ)Ko(ar) 


. outside dielectric core. 


J is the Bessel function and Ko is that modified 
Bessel function which behaves as a negative 
exponential at large arguments. In (5) the in- 
tensity factor has been made unity inside the 
dielectric and A outside dielectric. y is the axial 
propagation factor and p and a are the radial 
propagation factors related to the phase factor, 
B, by: 

vy’? +p? =e,B", 

7y?—a?=8?, 


(6) 


where ¢, is the relative dielectric constant of the 
material.‘ 

In this case, the formulas for E and H in 
cylindrical components become : 


Inside 
Es = —joup sin(yZ)Ji(or) 
H,= —~p cos(yZ)Ji(pr) 
H,=p? sin(yZ)Jo(pr) 


Outside 
E,= —Ajwya sin(yZ)K (ar) 
H,= —Avya cos(yZ)K;(ar) 
H,= —Aa? sin(yZ)Ko(ar). 


A must then be chosen to satisfy the boundary 
condition on the dielectric interface. This leads 
to equations for A which in matrix form would 


be: 
pJ (pa) — aK (aa) 1 e 
foo ee aie (8) 


3S. A. Schelkunoff, Electromagnetic Waves (D. Van 
Nostrand Company, Inc., New York, 1943), p. 425 ff. 
4 Rationalized Giorgi units are used. u is permeability. 
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Fic. 4. Functions for symmetrical modes of dielectric wave guides. 


Roots of the determinant yield the charac- 
teristic or eigenvalues of the system. 
To systematize the solution somewhat we de- 
fine the auxiliary quantities, x, R, by: 
x =pa (9) 
(The equations in boxes are those pertinent to 
the actual design.) 
R=p/a (10) 
The determinantal equation from (8) must be 


solved by numerical methods, and 
this we define the functions: 


to facilitate 


J \(x) 
g(x) =— | 
xJo(x) 
(11) 
K(x) 
(x) =———_ 
xKo(x) 


which are plotted to a log-log scale as Fig. 4. 
In these terms the determinental equation be- 
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comes: 


o(x) =y(x/R) (12) 


The design problem is usually posed in the 
following manner: ¢e, and 6 are given and it is 
desired to find a, the radius of the sample and J, 
the length of the sample. The degree of freedom 
left to the experimenter which was mentioned 
earlier can best be used to fix the value of R. 
Having done this, one solves (12) for the lowest 
root, x=pa. Then by means of (6) and (10) the 
desired quantities may be expressed as: 
2a x (1+R?)! 1 


X R 


21 (—) 
e-+ R? 


(e-—1)! 
(13) 


In these equations \ is the wave-length in air 
and / is such that the sample will be 4 wave- 
length long. 
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The value of A which is needed for later cal- 
culations now follows from (8) as: 


RJ; (x) R°Jo(x) | 
iets ES ie 
K,(x/R) Ko(x/R) | 


To set up convenient explicit directions for a quick 
application of the present method, a value for R has been 
chosen. I take: 


R=1.200 
which gives 

x =2.965, 

A =8.66. 


For this arbitrary choice of the degree of freedom con- 
tained in R, Eqs. (13) are plotted as a function of dielectric 
constant in Fig. 5. This finishes the calculation for the size 
of sample, since from Fig. 5 for any value of e, the sample 
dimensions can be immediately read off in terms of the 
air wave-length. 


CORRECTION FACTOR 


It was mentioned before that the value of do 
obtained from the measurements must be cor- 
rected to yield the true value of dielectric loss 
factor, d. To find this, it is necessary first to 
calculate the value of the stored energy, W. This 
involves a calculation of the integral of E? 
throughout the volume of the resonator. 

In terms of evaluated Bessel integrals: 


9 


By=—(JiX(#) ~ Jol) Jae) (15) 























Fic. 5. Curves giving size of sample as a function of e,. 
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(16) 
X (Ko(x/R)Ko(x/R) —K1?(x/R)) 
such a calculation gives: 
W =B*urh(e,B,;+ABz), (17) 


where h is the length of the resonator, i.e., / or 21. 

Similarly the loss factor due to the dielectric, 
dy, can be calculated and related to the loss 
factor of the dielectric, d. The latter is: 


d=vy¥,/we 


where 7, is the conductivity of the dielectric. 
Then the ratio d/do, i.e., the correction factor 
mentioned above can be expressed as: 





| ‘| 
| d A*B, 
| do eB, | 





(18) 


For our numerical choice of R=1.2000 we find: 


B,=1.075, 
B,=0.00284. 


For these numerical values the ratio (18) is plotted as 
the lower curve of Fig. 6. By means of this curve the ob- 
served value do can be immediately corrected to the desired 
value, d. 
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Fic. 6. Upper curve: ratio of maximum detuning to 
volume ratio of rod to sample. Lower curve: factor to 
correct for field outside sample. 
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SIZE OF RESONATOR 


Since the lid of the resonator slides freely in- 
side, the height of the metal cylinder is not 
relevant so long as it is large enough. Its di- 
ameter, 2b, must, however, be larger than two 
separate numbers: namely, it must be large 
enough so that the change in frequency due to 
its presence is small, and it must also be large 
enough so that the loss on the side walls is 
negligible. 

We first calculate a minimum diameter to 
assure only a negligible change in frequency. 
This calculation is carried out by imagining first 
an envelope placed at infinity which is subse- 
quently brought in to the finite radius, b. At 
very large values of the variable radius, r, the 
effect of the metal envelope is to cancel E and 
the perpendicular component of H, and also to 
double the parallel component of H. This can be 
proved by a direct calculation, but it is also 
obvious from the nature of electromagnetic 
reflection. Such an electromagnetic field being 
parallel to the wall exerts thereon a pure pres- 
sure, P, whose magnitude is given by: 


p=3u| H1|*. 


The change in frequency resulting from a 
movement inward of the surrounding envelope 
can be calculated from an action theorem! ac- 
cording to which the change in energy AW is 
related to the change in frequency Aw by: 


Aw/w=AW/W. (19) 


(For the purposes of this calculation the damping 
is ignored and it is assumed that the resonator 
is oscillating continuously without driving forces.) 

To find the change in energy of the field 
within the resonator when the envelope is brought 
from r= « to r=b, we calculate the work done 
‘against the radiation pressure which can be seen 
to be: 


D 


AW =2xhya'A *f K,*(ar)rdr. 


b 


This integral can be evaluated with sufficient 
accuracy by using for the modified Bessel func- 


Quar. App. Math. 2, 329-335 (1945). 
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tion Ko the asymptotic approximation : 


Ko(x) oe 


—, 0 
(wx /2)! (20) 


This gives for the fractional detuning, S, the 
value: 
Aw *a*A? c™ 


—__——_.. (21) 
€«-B,+A°B, 

However, the problem is not to determine the 
detuning but to determine the minimum value 
of 6 for an arbitrarily giyen value of detuning. 
Equation (21) can be re-arranged for this pur- 
pose to yield: 


(-) R : ( 2 
—-} =—In 
a/s 2x aB2S(1+ R?) 


€-—1 
nati 
1+(B,/A*B2)e, 


| 








(22) 


where the necessary minimum value of } as a 
ratio to a to satisfy the detuning condition is 
explicitly given. 

But 6 must also be large enough so that the 
loss on the side wall is negligible. This loss, D,, 
can be expressed in terms of an integral over the 
surface of the cylinder, and evaluated as: 

D,=wpb2rhbatA*Ky*(ab), 


where 6 is the depth of penetration into the 
metal. This value of loss yields for Q due to the 
shield, Q,=1/d,: 


B? €«-B,\+A*B, 
* 2bba4 A*Ko?(ab) 





(23) 


The problem posed however, as was the case 
before, is to find a minimum value of 6 for a 
given large value of Q,. By the use of asymptotic 
expression for Ko this can be evaluated as: 


-————-— anstiindnnmaemeriatiateniatinanemmmmnigpmaaets 


b R e-—1\3 
ne 
asq@ 2x 1+ R? 
8A7*5/rX 
x ——¢,| 
A*B.+e,B, 
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Fic. 7. Dotted curve: minimum ratio of resonator to 
sample diameter to assure frequency change of less than 
1 percent. Solid curves: minimum ratio of resonator to 
sample diameter to assure side wall Q in excess of half- 
million plotted for various vacuum \ in cm. 


For our choice of R=1.200, (22) and (24) have been 
plotted as a function of dielectric constant in Fig. 7 for 
the following arbitrary choice of S and Q,: 


S=1 percent, 
Q. =500,000, 


and where in (24) the value of 6 for copper has been used. 
(24) depends explicitly on \ which is involved in 6 and 
hence a family of the curves is necessary. 

A value for 6 must be chosen large enough to 
meet the requirements of both (22) and (24). 


TUNING ROD 


At the center of the resonator H has only a z 
component and E is zero. Hence a very small 
axial wire all or partly immersed in the resonator 
does not disturb this field appreciably. Consider 
the circumstance of a small rod of radius, ¢, 
projecting into the resonator a distance z. If this 
rod were shrunk to zero radius and then expanded 
to its finite but small radius it would produce 
practically no change in the magnetic field H,. 


However, in the process of expanding this rod, . 


work would be done against a radiation pressure : 


p=3|H,\?. (25) 
The total work done would be: 
aw-f prods. (26) 
0 


Therefore the total work done is found by 
evaluating (26), and if we use the expression (17) 
for W and also make use of the action theorem, 
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Fic. 8. Relative fractional detuning as function of fractional 
penetration of tuning rod for one and two samples. 


we can get for the fractional detuning, ¢, at 
penetration z the value: 


Aw 1\3 d, 2yz—sin 2yz 
t=—= (_-) (p°As)*— , (27) 
® 4n h ¢-B,+A*Bz 





wherein \, is the wave-length along the axis of 
the resonator. In our application the resonator 
mode used is either 7 Eo1; or TEo.2, or generically 
TEoin, where n is the number of samples used, 
either one or two. Hence: 


\,/h=2/n. (28) 


Let us call ¢o the detuning brought about by 
full immersion of the rod, i.e., z=. In these 
terms: 


eo —_—_——} 


| 


t ifs sin 2x(z/l) 
= ---—] I, (29) 
to ntl 2a | 

| 


where 
(p*ds)? 
¢,Bi:+A2By 


1 2 
“te 
4nr 

(z/l) is the fractional immersion of the rod 
and runs from 0 to 1 or from 0 to 2 depending 
on the number of samples. This function is 
plotted in Fig. 8. 

It is now convenient to introduce the ratio, 7, 
of the volume of the fully immersed rod, to the 
volume of the sample. This ratio is independent 
of the value of n. In terms of this quantity the 
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total detuning becomes: 


to x? R* e-—1 


icatie Si ianincaatiaiay ‘stein - (30) 
T 4A*B,1+R? 1+6¢,B;/A*B, 


For R =1.200 this latter function is plotted as the upper 
curve of Fig. 6. 

It is now also useful to check the Q due to the 
tuning rod. This is of course a variable quantity 
depending upon the immersion. However at its 
worst, i.e., lowest, which occurs for full immer- 
sion it can be calculated by similar processes to 


be: 
2a\* 247 /a’ 
0-(—) : (— )(cB:+4°B.) 
N x' \6¢ 


In the experimental application described, this 
Q, came out to be several hundred thousand i.e., 
negligible. 


CONCLUSION 


By accepting the arbitrary choice of R=1.200 
the curves of this article make it possible to 
proceed immediately to measurements by this 
method. Other curves could be calculated by the 
same process for other values of R which would 
yield samples of a different shape. 

In the practical application of this method the 
writer would like to call attention to one spurious 
phenomenon that may occur. Because of the 
configuration of the field of the TE ;, mode the 
tuning rod detunes the resonator but a very 
small amount: of the order of one part in a 
thousand. However it detunes other modes 
radically. It is therefore conceivably possible 

that in pushing the rod through the resonator 
some other mode becomes tuned to the operating 
-frequency. This mode will moreover probably 
not be confined to the dielectric and hence will 
couple the input and output loops very strongly. 
The practical result of such an occurrence would 
be that the galvanometer deflects sharply at one 
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highly critical position of the tuning rod. Unless 
this position coincides with a point at which a 
reading must be made it can merely be ignored. 
If it does so coincide the half-power points can 
readily be shifted to another region by a small 
change in the oscillator frequency. 

The equipment described in this paper was 
tested by measuring the loss factor of a sample 
of polystyrene. To determine the reproduc- 
ability of the measurements, a series of thirteen 
pairs of readings of d; and dz were taken and 
thirteen values for d calculated. These readings 
yielded for the loss factor a mean value of 
2.7 10-4 and an r.m.s. fluctuation of 0.1 10-4. 


As an example, the numerical values for the test sample 
and resonator and for one pair of measurements will be 
given: 
bedahe 


Desired \ = 10 cm; assumed e, = 
From Fig. 5: 


2a=9.9cem; /=3.92 cm. 


Two samples cut to size. [Actual measurement of samples: 
2a = 10.0, 1=3.92.] 

Radius of resonator chosen: 25=15.9 cm. Hence 
b/a =1.59. From Fig. 7 S<1 percent and Q,>}3 million. 
Sample resonated at \=10.03 checking closely assumed e. 


Tuning rod diameter 0.240 cm chosen, yielding 
T =1/1740. 


From Fig. 6: t9=3.91-1074. 
Crystal calibrated to half-power points and following 
readings taken: 


n=1 n=2 
‘ } 1.4 0.94 
4 power penetrations (cm) 2.74 15.83 
. : 0.28 0.24 
fractional penetrations 0.70 1.49 
1 
t/te from Fig. 8 Pony ae 
difference 0.723 0.70 
t=d, de 2.83X10™4 2.74X10"4 
calculation of do 2.65 X 104 
d from Fig. 6 2.86 10~* 


Since the stability of the oscillator used was poor, 
and since the means for measuring the penetra- 
tion of the tuning rod in the present equipment 
is not very accurate, it seems that better results 
could be obtained by better equipment. 
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Magnetically-Controlled Gas Discharge Tubes 


R. E. B. MAktnson, J. M. SOMERVILLE, K. RACHEL MAKINSON, AND P. THONEMANN 
University of Sydney, Australia 
(Received February 4, 1946) 


Discharge tubes are described in which a low pressure glow discharge is initiated by means 
of a magnetic field pulse. A tube has been developed in which a glow discharge so triggered 
leads to the formation of a cathode spot on a mercury pool; this tube will hold off 20 kv until 
fired and will then pass a current pulse of several hundred amperes. Its use in radar modulator 


circuits is envisaged. 


1. INTRODUCTION 


T is known '~* that the application of a mag- 
netic field of sufficient strength can initiate a 
discharge between charged electrodes in certain 
configurations, such as coaxial cylinders, when 
the gas or vapor pressure is too low to permit a 
discharge otherwise. The magnetic field lengthens 
the electron trajectories so that the average 
number of ionizing collisions suffered by each 
electron is sufficient for a glow discharge to be 
sustained. A detailed study of such glow dis- 
charges carrying continuous currents up to sev- 
eral amperes between coaxial cylinders has been 
made by Penning,* who has also used the same 
principle in the construction of a low pressure 
manometer and a low pressure ion gun.*® 
The work described below was mostly con- 
ducted with current pulses, the peak current 
being of the order of a hundred amperes and the 
discharge normally an arc, with consequent low 
voltage drop between the electrodes. The main 
object was the development of magnetically- 
controlled discharge tubes which would be suit- 
able for use as modulators in radar equipment. 
Both permanent gases and vapors were used, 
tube life in the former case being limited by 
clean-up of the gas. A satisfactory type of tube 
was developed having two coaxial cylindrical 
electrodes, together with a mercury pool as an 
auxiliary cathode* (Fig. (5a)); a glow discharge 
in the mercury vapor between the cylinders was 


1 Strutt, Proc. Roy. Soc. 89A, 68 (1913). 

2Engel und Steenbeck, Elektrische Gasentladungen 
(Verlagsbuchhandlung Julius Springer, Berlin, 1934), Vol. 
II, p. 236. 

3 Penning, Physica 3, 873 (1936). 

4 Penning,*Physica 4,871 (1937). 

5 Penning and Moubis, Physica 4, 1190 (1937). 

* The name “‘Hodectron” is proposed for tubes of this 
type. 
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initiated by a magnetic field pulse, and bombard- 
ment of the pool by ions from this glow discharge 
produced an arc spot on the pool through which 
the main current passed. These tubes passed 
currents up to 200 amperes in 2-megawatt pulses 
of duration 1 to 10 microseconds at repetition 
frequencies up to 1200 c/sec. with very small 
power loss in the tube itself. 

Although simple two-electrode concentric cy- 
lindrical tubes were not satisfactory as radar 
modulators it was found that in such tubes, if 
the glow discharge was not allowed to pass into 
an arc, currents up to 50 amperes could be 
switched off against supply voltages of a few 
thousand volts by removal of the magnetic field. 
It seems therefore that such tubes may be useful 
where a rapid high current or high voltage 
switch is needed. 
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Fic. 1. (a) Showing pulse-forming line and load, R; 
(b) Arrangement of biased diode to hold off H.T. from 
oscillograph. 
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Fic. 2. Examples of two-electrode tubes. 


2. CIRCUIT ARRANGEMENTS 


The magnetic field was produced by a short 
coil coaxial with the electrodes and was applied 
in pulses of about 50 ysec. duration at repetition 
frequencies which were usually 50 or 100 per 
second; the peak value of the field was of the 
order of 500 oersted. The current in the coil was 
produced by charging a condenser from a d.c. 
source and discharging it through the coil by 
means of a vibrator or mercury thyratron. 

The performance of the tubes was studied 
mainly in the circuit shown in Fig. 1 (a), which 
approximated that in which they would be re- 
quired to operate as radar modulators. In this 
circuit the resistive load R corresponded to that 
_ placed on the modulator by the radar oscillator. 
The condensers C, forming with the inductances 
L, and Lz an artificial line, were charged through 
the hold-off resistance S to between 1 and 20 kv, 
and were discharged through the tube and the 
load R on the initiation of a discharge by the 
magnetic field pulse. The variation of potential 
difference between the electrodes was studied by 
means of a triggered oscillograph, a diode being 
used as shown in Fig. 1(b) to hold off potentials 
greater than a few hundred volts. The oscillo- 
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graph was also used to study the current through 
the tube, represented by the voltage drop across 
the small resistance r. 


3. GENERAL NATURE OF THE DISCHARGE 
IN TWO-ELECTRODE TUBES 


Examples of two-electrode tubes studied are 
shown in Fig. 2. Essentially the electrodes were 
concentric cylinders; with the outer cylinder as 
cathode a smaller magnetic field was required to 
initiate the discharge, so this arrangement alone 
was used. Copper, aluminum, steel, and nickel 
were used as cathode material, copper, nickel, 
and tungsten for the anode. Tubes were filled 
with helium, argon, air, and carbon dioxide and 
with iodine, mercury, and water vapors. 

The tubes were initially studied while still 
sealed to the vacuum system, so that the pressure 
could be controlled. The pressure at which a dis- 
charge could be initiated by a magnetic field 
depended on the configuration of the electrodes 
and the nature of the gas or vapor, but in general 
was of the order of 0.1 mm for permanent gases. 
The range of pressures for a controlled discharge 
in one tube was about 10 to 1. 

When the applied magnetic field reached its 
critical value the potential drop across the tube 
in general fell rapidly and the condensers C dis- 
charged completely in a few microseconds. The 
line circuit constants determined that the cur- 
rent through the tube should flow as an approxi- 

















/¢ 
an 7 
oe 
all 
— 
= 4 
ra 
S 
z ‘ 
ye 
. / ed 
» # 
§ \ 
: 
ey 
he 
ize 
x 
a 
€ 
g 
= ¢ 
q 7) ir) s34V 

















ANODE VOL TAGE 


Fic. 3. Variation with anode voltage of the peak current in 
the field coil necessary to fire a two-electrode tube. 
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Fic. 4. Typical oscillograph traces of the voltage across the tube during discharge. 


mately square pulse of a few microseconds dura- 
tion; with the circuit used a small reverse current 
sometimes flowed after the main pulse. The cur- 
rents passed were of the order of 100 amperes. 
The hold-off resistance S was large enough to 
ensure that the condensers did not charge up to 
breakdown potential twice during the same mag- 
netic field pulse. The variation with anode volt- 
age of the minimum field-coil current necessary 
to fire the tube is shown in Fig. 3. 


Voltage Drop in Tube 


When visual inspection was possible the dis- 
charge was seen to fill the whole inter-electrode 
space. It was usually associated with a bright 
white spot which moved about on the cathode 
and particularly favored a metal-insulator junc- 
tion. Most of the current passed through this 
spot. When the spot was present the voltage 
drop across the tube was between 30 and 100 
volts, in its absence about 400 volts. It was con- 
cluded that the presence of the spot meant that 
the discharge was an arc. Rapid alternation be- 
tween arc and glow in successive discharges 
sometimes occurred. 

Typical oscillograph traces of the voltage 
across the tube during the discharge are shown 
in Fig. 4. The initial deep valley ABC shown in 
Fig. 4(a) is spurious, and is caused by the rapidly 
falling anode potential acting through the inter- 
electrode capacity of the protecting diode. The 
portion CD shows the potential across the tube 
while it is conducting (the potential drop across 
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the diode being negligible). When the tube be- 
comes non-conducting, the characteristics of the 
line produce the portion DEF of the trace. The 
valley ABC and the portion DEF have been 
omitted from Figs. 4 (b), (c), (d), and (e). 
Figures 4(a) and (b) are typical of arc discharges, 
the voltage dropping rapidly to a value of about 
50 volts, which does not vary greatly either with 
current, over a range of 100 to 1, or with the 
nature and pressure of the gas in the tube. Fig- 
ure 4(c) is typical of a glow discharge, the voltage 
usually falling to between 300 and 600 volts, 
although much ‘higher values have sometimes 
been observed; these traces vary markedly with 
tube conditions and current. Figure 4(d) shows 
a delayed transition from glow to arc. 

The arc discharge in a helium-filled tube with 
a copper cathode was a characteristic pink color 
and photographs of the spectrum showed both 
He and Cu lines. In tubes of this type another 
kind of discharge was also observed, bright 
green, with a voltage drop of the order of 700 
volts and a spectrum showing more intense cop- 
per lines than those of the arc discharge. 

The transition from the glow to the arc was 
favored by irregularities or discontinuities in the 
electrodes or the magnetic field and by an in- 
crease in the supply voltage (and consequently 
in the tube current). The glow-arc transition 
showed no marked dependence on the nature 
and pressure of the gas. 

Owing to the low power dissipation in the 
tube with the arc discharge, this mode is to be 
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preferred. Consistent striking of the arc during 
the whole life of the tube was not obtained using 
a solid cathode. Tubes are described in Section 4 
in which this difficulty was overcome by using a 
mercury pool as cathode. 

Figure 4(e) shows an unstable type of voltage 
trace, frequently observed, in which the voltage 
was subject to rapid, apparently random fluctua- 
tions, no two consecutive traces being the same. 
Because of its appearance on the oscillograph 
such a trace has been called “grass.’’ The height 
of the peaks varied from a few volts to about 200 
volts. In tubes of the mercury pool type shown 
in Fig. 5(b) (see Section 4), grass occurred only 
when the tube temperature was below about 
40°C. With permanent gases the appearance of 
grass was unaffected by changes in pressure and 
was also insensitive to changes in the magnetic 
field strength. No satisfactory theory of its oc- 
currence has been evolved, but it is suggested 
that it may be associated with the motion of the 
cathode spot; alternatively its cause may be in 
the mechanism of formation of the arc spot. 


Clean-Up of Gas 


Considerable clean-up of gas took place during 
the discharge, severely limiting the life of any 
sealed-off tube which depended for its satisfactory 
operation on the presence of permanent gas. 
Some reduction of the clean-up was obtained by 
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Fic. 5. Three-electrode tubes with anode, starting 
cathode, and mercury pool cathode; (a) with external 
resistance, (b) with internal resistance. 
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preliminary saturation of the walls with gas and 
by the addition of a buffer volume, but the long- 
est life so achieved was only 50 hours. The 
effects of clean-up were avoided in tubes de- 
scribed in Section 4, which would start running 
on mercury vapor at room temperature. At- 
tempts to run tubes on iodine and water vapor 
were unsuccessful. 


4. ARC DISCHARGES IN MERCURY POOL TUBES 


By addition of a mercury pool to the coaxial 
cylinder electrodes described in the preceding 
section a tube was obtained which operated on 
mercury vapor at room temperature and in 
which the discharge always took the form of an 
arc. Such tubes are shown in Fig. 5. The mercury 
pool acted as the cathode taking the main cur- 
rent and the outer cylindrical electrode merely 
as a starting cathode. Between the starting 
cathode and the anode a glow discharge was 
initiated by the magnetic field; bombardment 
of the mercury by ions from this discharge 
caused an arc to strike on the mercury surface. 

These tubes operated satisfactorily at ambient 
temperatures between a few degrees C and 80°C, 
but below about 30°C unsteadiness occurred in 
the current trace on first switching on, disappear- 
ing as the tube warmed up. To ensure that the 
arc regularly formed on the mercury and not on 
the starting cathode, a resistance was placed in 
the starting cathode circuit. This was either an 
external resistance as in Fig. 5(a), or an internal 
resistance as in Fig. 5(b) where the starting 
cathode was a thin coat of Aquadag or sputtered 
tungsten. The latter type of tube had a life of 
only a few hours under the test conditions, i.e., 
as in a radar modulator, failure resulting from 
one of two causes; sometimes the action of the 
arc spot as it moved around the line of contact 
of the mercury and the glass separated the coat- 
ing from the mercury pool; sometimes sputtering 
from the anode during current reversal after 
passage of the main pulse so lowered the re- 
sistance of the starting cathode that the arc 
spot formed on it instead of the pool, the voltage 
drop in the tube being then much higher and 
arcs sometimes forming again between the mag- 
netic field pulses. 

Of all types tested, the tube of Fig. 5(a), with 
the external resistance between starting cathode 


JOURNAL OF APPLIED PHYSICS 























and pool, gave the most satisfactory operation. 
Its life was limited by the formation of a con- 
ducting path of sputtered material and mercury 
between the pool and the lead to the starting 
cathode, with consequent arcing on the starting 
cathode as described above. Such tubes have 
handled powers up to two megawatts from a 
power supply of 10-20 kv at repetition fre- 
quencies up to 1200 c/sec. with pulse lengths from 
1 to 10 usec. One had a life of 150 hours; the con- 
ducting path mentioned above could then still 
be broken up by shaking. The interval between 
the beginning of the magnetic field pulse and the 
beginning of the current pulse through the tube, 
which was from 20 to 50 usec., never varied from 
pulse to pulse by more than 2 usec.; the varia- 
tion decreased as the supply voltage was in- 
creased, its average value being less than 1 
usec. with a supply voltage of 10 kv. A mag- 
netic field pulse of peak intensity 350 oersted 
was sufficient to fire the tube. 

The value of the resistance R; was not critical, 
values from 0.01 to 0.5 megohm being satis- 
factory. Operation was not affected by tilt of 
the tube so long as the mercury did not touch 
the starting cathode. Alignment of the electrodes 
was not critical, nor was that of the field coil. 

Graphite, which does not amalgamate with 
mercury and sputters at a slow rate in mercury 
vapor, was found to be a suitable material for 
starting cathode and anode. The clearance be- 
tween these two electrodes and the pool was 
made greater than }”’, since “‘stalactites”’ of mer- 
cury sometimes grew downwards from them 
towards the pool. For a similar reason, the 
clearance between the starting cathode and the 
glass was kept greater than 35”. , 

In operation of the tube, gases are cleaned up 
and it was found that rigorous outgassing or 
very low pressures were not necessary before 
sealing off. 

Figure 6 shows the oscillograph traces of tube 
current, anode potential, and starting cathode 
potential; these may be interpreted as follows. 
When the magnetic field reaches the critical 
value a glow discharge is initiated between the 
anode and the starting cathode, and the poten- 
tial of the latter rises. About 0.7 usec. later the 
arc strikes on the mercury, and the starting 
cathode potential drops nearly to the cathode 
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Fic. 6. Oscillograph traces of tube current, anode potential, 
and starting cathode potential. 


potential (zero), remaining at this value while 
the discharge continues, and afterwards so long 
as the anode is negative. When the anode be- 
comes positive the starting cathode potential 
also rises until deionization in the inter-electrode 
space is complete, when it falls again to zero. 

It was observed that on first operating new 
tubes of. this type the arc spot at first formed on 
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Fic. 7. Tubes with electrodes other than coaxial cylinders, 


the free surface of the mercury pool, then formed 
only with difficulty until (after about 10 minutes) 
the mercury had begun to wet the glass; the 
spot then formed regularly at the meniscus. 


Other Mercury Pool Tubes 


During work leading up to the design of the 
tube just described, it was found that discharges 
could be initiated between electrodes whose 
geometry differed widely from the concentric 
cylindrical types and the configurations used by 
Penning.*® Some examples are shown in Fig. 7, 
all of which required a permanent gas in addi- 
tion to mercury vapor for operation at room 
temperature and were thus subject to clean-up 
troubles. 


5. INTERRUPTION OF CURRENTS BY REMOVAL 
OF THE MAGNETIC FIELD 


The general features of the glow discharge 
initiated by a magnetic field have been discussed 
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in Section 3. It was found that, so long as the 
discharge did not pass into an arc, reduction of 
the magnetic field below the critical value in- 
terrupted the current through the tube, causing 
it to fall to a negligibly small value in less than 
one microsecond. Once the arc had been formed, 
however, removal of the magnetic field caused no 
change in the tube current, which continued to 
flow until the condensers were discharged. 

Some attention was given to the production 
of a tube which would switch moderate or heavy 
currents off as well as on. It will be seen that the 
chief problem involved in this was to prevent 
the transition of the discharge from a glow to an 
arc. Concentric cylinder two-electrode tubes 
were used, and to prevent arc formation it was 
found necessary to keep the supply below about 
3000 volts; this limit did not depend much on 
the current. Currents up to 50 amperes were 
interrupted in this way. j 

It is possible that in a tube with a large cathode 
area, containing a gas such as nitrogen, which 
has a large normal cathode current density, large 
currents could be switched on and off, passing 
with a few hundred volts drop across the tube. 
A cold cathode tube might thus be obtained 
which would do what is now only done with 
difficulty by high vacuum hot-cathode tubes. 

We wish to acknowledge the extensive co- 
operation and help of Professor V. A. Bailey, 
the C.S.I.R. Radiophysics Laboratory, Sydney, 
and of Professor L. H. Martin, L. G. Parry, and 
W. R. Watson. While this work was being car- 
ried out, one of us (J.M.S.) was a holder of a 
Commonwealth Research Fellowship and an- 
other (K.R.M.) was employed by the Radio- 
physics Laboratory. 
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Depreciation of Fluorescent Lamps without Passage of Current 


H. C, FROELICH 
Lamp Development Laboratory, General Electric Company, Cleveland, Ohio. 
(Received April 9, 1946) 


Fluorescent lamps with thin and with partial coatings of phosphor on bulbs of ultraviolet 
transmitting glass were irradiated with Hg 2537 radiation. Comparison with similar lamps 
which were burned showed a decline in lumen output which was of the same general character 
in all experiments, except for different time factors. This indicates that the process of lamp 
depreciation is fundamentally of photochemical nature, involving excited mercury, but not 
requiring either Hg ions or neutral atoms. It is assumed that superficial oxygen ions of the 
phosphor crystal react with Hg* to form links of the type-SiO-O-Hg-Hg..., thus tying 
mercury to the crystal with a resulting loss of phosphor brightness. 


INTRODUCTION 


T is common knowledge that the light output 

of fluorescent lamps declines on burning. It 
has also been shown that this depreciation is 
caused by a screen of condensed mercury and 
mercury compounds in an extremely thin though 
discontinuous surface spread; it is not the result 
of a breakdown of the phosphor crystals in 
their bulk.! 

It seemed to be of fundamental interest to 
find out whether this depreciation is necessarily 
associated with the passage of current through 
a lamp, or whether it could also be brought about 
simply by a photochemical reaction involving 
excited mercury or perhaps excited phosphor. 
In the absence of mercury the phosphors do not 
depreciate. If a photochemical reaction should 
occur between mercury and the solids of glass 
wall or phosphor, then a new approach to the 
study and understanding of phosphor quality 
and “stability” (lamp maintenance) would have 
been found. Such reactions would have to be 
confined to the surface of the crystals because 
the bulk (or interior) of the crystals remains un- 
changed in brightness, and also because of the 
ease with which phosphors can be demercurized. 

The problem, then, was to expose phosphors 
to the action of excited mercury and follow their 
brightness changes with time. Two types of de- 
vices were used for the purpose, the limitations 
of which with respect to rigorous treatment will 
be pointed out later. However, the accuracy 
desired for this exploratory work is believed to be 
sufficient to establish either the presence or ab- 


1H. C. Froelich, Trans. Electrochem. Soc. 87, 365 (1945). 
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sence of the effect; and the simplicity of the set- 
ups was a persuasive factor in deciding against 
more complicated apparatus. 


EXPERIMENTAL 


The first type of experiment has already been 
described recently.' In a jacketed lamp device 
the phosphor was irradiated from within, ex- 
posed to both \2537 and excited mercury vapor. 
The light output fell 25 percent during about 600 
hours of irradiation, but the recovered phosphor, 
after proper demercurization, was found to be as 
bright as it was originally.* 


* External to a lamp, the brightness of phosphor powder 
was measured as follows: A layer of powder about 2 mm 
thick, with flat and well-defined surface area, was irradiated 
with a beam of A2537 under oblique incidence. The 
amount of visible light generated by the excited phosphor 
was then measured with corrected photo-cells in terms of 
the brightness of an arbitrary standard of equal area, 
measured similarly and immediately before or after the 
sample. The accuracy of measurement was 4} percent or 
better. A differential method might be more desirable. 

Unfortunately, the mercury resonance line \2537 is the 
only monochromatic radiation which can be conveniently 
and abundantly produced in the spectral region where 
phosphors for use in F lamps are most sensitive. This fact 
necessitates complete absence of mercury and its com- 
pounds, solid, liquid, or gaseous, before the phosphor can 
be excited for correct measurement. With an equally good 
source of, say, 42500 or 42600 radiation, disturbance by 
mercury vapor would be eliminated but a reduction in 
brightness of non-demercurized phosphor would still be 
found. It would be in proportion to the area and density 
shaded by condensed Hg and Hg compounds, provided 
the latter are opaque to that radiation, which is likely to 
be the case. Demercurization is easily effected by heating 
the phosphor in a good vacuum to low temperatures 
(300°C or even as low as 150°C). 

Removal of phosphor from a lamp involves complete 
random reorientation of the billions of individual crystals 
comprising the coating. Such bulk measurements, then, 
showed no significant difference in brightness between 
demercurized phosphors from unburned lamps and those 
from lamps several hundred or thousand hours old. 
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In the second type of experiment the same 
principle was used but the technique was re- 
versed. The phosphors were irradiated from the 
outside, through ultraviolet transmitting glass. 
Phosphor brightness was measured in terms of 
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Fic. 1. Depreciation of half-coated lamps. 


lamp brightness with current passing through 
the lamp, but without disturbing the orientation 
and the bulk of the phosphor crystals. Regular 
“white” phosphor was coated on 18” T8 germi- 
cidal bulbs of Corning No. 9741 glass, and 
regular lamps were made. The lamps were 
seasoned, photometered in a sphere for their 
lumen output, and then divided into several 
groups. One group was stored to be re-read at 
intervals to check accuracy. Another group was 
burned_normally and also photometered at cer- 
tain intervals. The last group was irradiated 
with (2537 in the following manner: Each lamp 
was surrounded by 6 regular germicidal lamps 
which were burned in order to provide a fairly 
uniform and reasonably high amount of A2537 
radiation. No attempt was made to measure the 
wattage and density of \2537 radiation actually 
entering the test lamp. It was estimated, how- 
ever, that it should be of the order of one-half 
to one-fourth of the amount generated in a 
normally burning lamp. The lamp to be irradi- 
ated was not burned except for the time when it 
was being measured. 
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In the first series the irradiated lamps were 
only half coated—i.e., half around their circum- 
ference but full length. In the second run, a 
uniform coating of less than normal thickness 
was applied all around the bulbs. A few controls 
of fully coated lamps in germicidal as well as 
lime glass were also made which were only 
burned. All lamps were photometered in a 
sphere at constant current of .300 amp. The 
data obtained are averaged and considerably 
condensed in Tables | and II. 

Figures 1 and 2 show plots of the data in 
Tables | and II. It is apparent that the irradiated 
lamps decline in their output and follow the same 
qualitative pattern as the burned lamps. Be- 
cause of the nature of the 9741 glass the deprecia- 
tion curves of burned lamps cannot be expected 
to follow the same uniform pattern as regular 
lead or lime glass lamps. The scale of the time 
ordinate for the irradiated lamps in groups C 
had to be reduced in order to make the similarity 
apparent; this was pointed out above. The re- 
duction factor in Fig. 1 was 3.6, while that for 
Fig. 2 was 2.6, as found by trial. To make the 
analogy more complete, the irradiated lamps 
were rejuvenated by application of external heat 
in a tube furnace at 300°C. The averaged lumen 
values were as follows: 


After 2 hours After 1501 hours Rejuvenated 
409 302 427 


The apparent slight improvement of the re- 
juvenated lamp over the 2-hr. reading is be- 
lieved to be caused,-at least in part, by a slight 
deterioration during the first 2 hours. 

Finally, a test was devised where depreciation 
and no depreciation could be produced in a single 
lamp. A 36’’.T8 lamp was made with a thin 
coating of white phosphor on a bulb of 9741 

TaBLe I. First run with half-coated lamps of normal 


coating thickness. (9741 glass). Averaged lumens at 0.300 
amp. 





Depreciation due to burning for 
2 44 204 


Lamps 421 hours 
A. Full coat 675 600 523 480 
B. Half-coat 452 367 359 323 


Depreciation due to irradiation for 

2 736 ©1501 hours 

C. Half-coat 409 361 311 302 
Hours of irradiation, 

divided by factor 3.6 46 204 
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glass. One-half of the lamp was surrounded ‘by 6 
regular germicidal lamps, while the other half 
was surrounded by 6 lime glass lamps which 
generated light and heat but no external ultra- 
violet; these lamps were also of the 18” T8 (15- 
watt) size. All 12 lamps burned in order to pro- 
duce nearly uniform radiation and temperature 
distribution within and around the test lamp in 
the center. Again, the test lamp did not burn 
except for the time when it was being measured 
after removal from the “light bath.” It proved 
impractical to take sphere readings of lumen 
output. The two ends of the lamp, therefore, 
were photometered against a known standard by 
taking readings with a cell box arrangement 
which caught the light from about two-thirds of 
each section. After each interval of irradiation 
the test lamp was burned for at least 15 minutes 
before output readings were taken. The lamp 
was photometered at constant current of 0.350 
amp. and under these conditions the meter de- 
flections should be proportional to the lumen 
output. Initial readings for the two ends of the 
lamp were not equal because a coating thickness 
differential existed. It was decided to irradiate 
the thicker coated, brighter section of the lamp 


TABLE II. Second run with full coated lamps of less 
than normal coating thickness. (9741 glass). Averaged 
lumens at 0.300 amp. 








360 928 hours 








A. Normal coating, 
lamps burned 
B. Thin coating, 
lamps burned 550 
C. Thin coating, 
lamps irradiated 539 


607 536 
420 


Hours of irradia- 
tion, divided by 
factor 2.6 8 117 357 


with \2537, because this might produce less 
favorable results (lower depreciation). 

In addition to several other effects, a notice- 
able change of color was observed in the section 
which was irradiated with ultraviolet; the color 
temperature went up by a few hundred degrees 
K (estimated). The decline in output, however, 
was several times greater than might be expected 
on that basis alone. This and the other effects 
observed require further study. 
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Averages of several measurements are given 
in Table III and plotted in Fig. 3. Again it is 
apparent that the light output of the ultraviolet- 
irradiated section declined as though the lamp 
had been burned with normal passage of current, 
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Fic. 2. Depreciation of thin coated lamps. 
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while the other half of the lamp remained sub- 
stantially unchanged. The slight departure (up 
and down) from a straight horizontal line re- 
flects the accuracy of measurement. 


DISCUSSION 


The experiments have given two results which 
are of fundamental interest. First, they have 
shown that fluorescent lamps can depreciate in 
the presence of excited mercury alone, without 
actual passage of current and in the absence of 
positive ions. Secondly, they have shown in a 
semiquantitative way that the depreciation due 
to excited mercury follows the same pattern as 
the depreciation of burning lamps: the decline is 
relatively greater during the initial periods than 
during the later stages. Excited mercury was 
produced by irradiation of its vapor with its 
resonance radiation \2537. This seems to prove 
that the phenomenon of F lamp depreciation is 
primarily of photochemical nature, involving ex- 
cited mercury and/or excited phosphor. It would 
be of great interest to separate the two possi- 
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bilities in repeat experiments with other than 
mercury resonance radiation. However, the latter 
should be near the peak of phosphor sensitivity, 
i.e., close to 42537. 

One of the chief limitations of the present ex- 
periments lies in their semiquantitative nature. 
In plotting the curves with a linear time abscissa 
it is assumed that the wattage input of 2537 
into the test lamp remains constant. Actually, 
the output of the germicidal lamps declines grad- 
ually (about 25 percent) during the course of the 
experiment, and the transmission to \2537 of the 
glass of the test lamp also declines through 
solarization. Therefore, the amount of \2537 
actually getting into the test lamp becomes 
progressively smaller. In a quantitative experi- 
ment this amount should be determined, at 
least from time to time, and the output curves 
should be plotted against watt-hours of \2537 
rather than against time alone. Such a plot will 
' show a change of absolute slope compared with 
the present data, but qualitatively the decline 
pattern will remain unchanged. The scope of the 
present paper did not attempt to go beyond 
these limitations. 

What are the forces which cause the deposition 
of mercury and its compounds on the phosphor, 
thus masking in part its otherwise unaffected 
capacity as a light transformer? One is forced to 
assume that a surface reaction (or, if one wishes, 
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“deterioration’”’) does take place between the 
phosphor or glass and excited mercury atoms. 
This reaction does not proceed into the interior 
of the crystals, but is confined to the surface; it 
requires a relatively low energy of activation, 
and is reversible. If it were to go deeper, it would 
necessarily result in a physical destruction of the 
phosphor. Recovery of initial phosphor bright- 
ness indicates that such disintegration does not 
take place. Exact measurements of the extinc- 
tion coefficient for 42537 in silicate, etc., phos- 
phors would be another means of following the 
stability of these crystals. To the author’s 
knowledge, such measurements have not been 
made with fine powders and \2537 radiation. 

In the absence of mercury, phosphors do not 
depreciate; in its presence, they do it to varying 
degrees. There are phosphors which give good 
maintenance in a lamp, and others of about the 
same chemical composition which show bad 
maintenance. The tangible and intangible fac- 
tors which play an additional role in the art of 
phosphor manufacture cannot be discussed here. 


TABLE III. Averaged meter deflections at 0.350 amp. 








2 162 647 1625 hours 





Ultraviolet-irradiated 
section 63 55 51 51 
Visible-irradiated section 45.5 46.2 46 45 
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Suffice it to say that chemical composition alone 
is not a criterion for phosphor performance in a 
lamp. If a phosphor gives good maintenance, less 
mercury reacts with it and deposits on it com- 
pared with phosphors of poor maintenance. One 
might even speak of ‘‘mercurophobic”’ and “‘mer- 
curophilic”’ phosphor surfaces although, strictly 
speaking, a lamp coated with phosphor of true 
mercurophobic surface should have no deprecia- 
tion whatever. 

It is easy to eliminate unexcited, atomic 
mercury as the agent responsible for deteriora- 
tion. Atomic mercury could reach the phosphor 
surface through condensation from the vapor 
phase, as a wetting or non-wetting continuous 
liquid or in form of innumerable discontinuous 
droplets. Either assumption can be shown to be 
absurd. Consider, for example, a regular 40 
watt, 48712 White lamp which has an inner 
(geometric) surface of about 1320 cm?, contains 
about 2.2 g of phosphor and about 25 mg mer- 
cury. If the mercury were spread out into a uni- 
form cylindrical film, it would form a layer of 
about 52 atoms or 156A thickness. It is incon- 
ceivable that the large surface tension forces 
would permit the free existence of a film stretched 
to that extent. Electron microscope pictures of 
metal films in thicknesses up to 100A or more 
show that they are completely broken up into 
separate particles.? Differences in thermal co- 
efficients of expansion of metal and substrate 
likewise would cause rupture and coagulation of 
the ‘‘film.”” The concept of a continuous film of 
mercury seems untenable. 

Simple condensation in form of innumerable 
tiny droplets also has its limits. If it is assumed 
that the normal 10 or 20 percent lumen loss 
of a fluorescent lamp be due merely to the shadow 
cast mechanically by subdivided mercury drop- 
lets, then the 25 mg of Hg in a 40-watt lamp 
should be expected to cover an area of 10 to 20 
percent of the lamp surface or 132 to 264 cm’. 
Assuming the droplets to be spheres, it can be 
calculated that this would require a subdivision 
into particles of the order of .1- to .05-micron 
diameter. Such small drops, however, would 
have a greatly increased vapor pressure compared 
with the larger drops from which they form. The 





2 R. G. Picard and O. S. Duffendack, J. App. Phys. 14, 
291 (1943). 
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well-known Thomson equation permits calcula- 
tion of these increases in vapor pressure. 


dp 20M 
p % rXsp.g.X RT 





(ong=450 dynes/cm; M=200; r=5X10-* cm; 
sp.g. = 13.55; 
R=8.3X10' erg/degree mole; T = 323). 


The equation holds so long as the surface ten- 
sion oug is independent of the radius; in other 
words, r must be large against atomic dimensions. 
Substituting the proper numerical values it fol- 
lows that a mercury particle of about .1-micron 
diameter (1000A or about 330 Hg atom diam- 
eters) would have a 10 percent higher vapor pres- 
sure at 50°C than a drop of 1-mm diameter. 
Such small particles would be thermodynamically 
unstable and would disappear by distillation. A 
few large particles would grow at the expense 
of most small particles. It would follow, then, 
that a lamp whose lumen output declined 10 or 
20 percent should regain its initial brightness if 
it were simply left lying for weeks, months, or 
years; this has not been found by experiment. 
Conversely, a new lamp has likewise not been 
found to depreciate if it was neither burned nor 
irradiated. 

Mercury droplets of thermodynamically stable 
size can account for not more than about 1 
percent lumen loss if the latter is caused by 
mechanical shading of the phosphor against 
\2537. The calculation left entirely untouched 
such energetic difficulties as are connected with 
seed and phase formation when vapor atoms con- 
dense to form stable sub-micron particles. The 
situation becomes even more difficult if the 
(geometric) phosphor surface is considered in- 
stead of the glass surface. A 40-watt (G.E.) 
lamp contains about 2.2 g of white phosphor of 
sp.g. 4.5. Particles of all sizes from 1 micron and 
lower, up to 4, 5 microns or more are represented. 
Their total number is of the order of 10", and 
their total (geometric, not true) surface is of the 
order of 10* or 105 cm*—many times the area of 
the glass on which they rest. To summarize, the 
depreciation cannot possibly be produced by 
chemical action or by condensation from the 
vapor phase of neutral, atomic mercury. 
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This leaves excited mercury atoms as responsi- 
ble for a reaction with the phosphor surface which 
can tie them up in a form in which their vapor 
pressure is lower than that of free Hg at the same 
temperature. Excited mercury could attack mer- 
curophilic phosphors or compounds by a reaction 
of the following types: 


ZnO + Hg*—HgO+Zn+ EE, 
SiO. +Hg*—> —SiO-—O—Hg+E2 
~Si0-—O-—Hg+Hg— 
—SiO —O—Hg—Hg-—-:--, etc. 


If we abandon the narrow chemical viewpoint 
‘of oxides and consider the phosphor or glass 
lattice as a whole, it is known that the latter 
consists of SiO, (or WQ,, etc.) tetrahedra with 
interspersed Zn, Mn, etc. ions or atoms. The 
crystal surface layer contains more ions of oxy- 
gen than of any other species, and unsaturated 
oxygen valencies may conceivably become tied 
up with Hg*. Existing stronger bonds may even 
be broken by the energy available from excited 
Hg (4.9 ev or 113 Cal.). Excess energy would be 
dissipated as heat as indicated by E in the equa- 
tions. Thus the oxygen atoms may become a 
bridge linking Hg atoms to the bulk of the phos- 
phor crystal by semi-chemical forces. Further 
growth involving —Hg—Hg bonds is quite 

plausible. After one or more additional “‘layers’’ 
of Hg atoms, if such should be present, these 
semi-chemical forces could be assumed to de- 
teriorate into physical adsorption, providing a 
condensed and opaque phase of sufficient density 
and stability to account for 10 to 20 percent ab- 
sorption of \2537. Thus, only a fraction of the 
25 mg of mercury in a 40-watt lamp would be 
required for the entire process of deterioration, 
and this is borne out by experience. Visual in- 
spection of lamps of many hundred or thousand 
-hours burning indicates that a good portion, if 
not the bulk, of the mercury is still present in 
form of free droplets. The condensed phase of 
mercury would also provide an effective barrier 
_Tor more excited atoms to react with the silicate 
surface underneath; the latter is then protected 
against further deterioration and the loss of 
brightness should occur at a greatly diminished 
rate, which is also observed. Complete clean-up 
of mercury by the phosphor occurs only when the 
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amount of mercury originally present was very 
small, or for other non-pertinent reasons. 

It should be noted that this proposed mech- 
anism of deterioration differs from somewhat 
similar suggestions which, however, go so far as 
to assume complete chemical reduction of the 
crystal surface to elemental zinc, silicon, ete. 
According to this theory, the oxygen is disposed 
of by assuming that it is liberated in gaseous 
form and is then free to diffuse and react else- 
where in the lamp. The loss of brightness is 
believed to be caused by the absorption of \2537 
by this opaque coating of elemental zinc, silicon, 
etc. The present tests indicate that this assump- 
tion is erroneous. The recovery of initial lamp 
brightness after simple application of external 
heat of depreciated lamps speaks against the 
presence or formation of material with a high 
vapor or decomposition pressure. Moreover, sili- 
con in its elemental form is a very stable sub- 
stance which reacts with oxygen only at rather 
high temperatures. It is hard to visualize a re- 
action which would make it innocuous when the 
phosphor is heated in a vacuum to 150°C. On 
the contrary, the tests on the full recovery of 
phosphor and lamp brightness indicate that a 
permanent deposit could not have formed on 
the phosphor. 

The change in color of the ultraviolet-irradi- 
ated white lamp, which was mentioned above, 
can be understood by assuming that the oxygen 
atoms or ions in tungstates are bound more 
tightly than in silicates. A relatively stronger 
attack by excited mercury atoms on the Zn-Be- 
silicate phosphor compared with the magnesium 
tungstate component would result in a shift of 
color temperature toward higher values. 

These speculations have been advanced chiefly 
in order to show what possibilities exist for chemi- 
cal reactions which would stop at the crystal 
surface, tie the mercury to it, yet make it re- 
coverable by such mild treatments as heating 
in a vacuum to low temperatures. With reference 
to the performance of fluorescent lamps, on the 
other hand, the term “depreciation” covers 
quite a complexity of processes. Each of these is 
contributory, though not individually responsible 
for the entire loss of lumen output. They include 
formation and deposition of a number of mercury 
compounds produced from excited mercury and 
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a variety of gaseous trace impurities left or 
formed in the lamp. This has been discussed in 
detail by Kenty and Cooper.’ 


SUMMARY 


It has been shown that the process of lamp 
depreciation is fundamentally of photochemical 
nature, involving excited mercury atoms and the 
surface ions of phosphor or glass. Fluorescent 


’C, Kenty and J. Cooper, Trans. Electrochem. Soc. 87, 
287 (1945). 


lamps containing mercury could be depreciated 
simply by irradiation with the mercury resonance 
radiation 2537, without passage of current. 
Semiquantitatively, this depreciation is of the 
same magnitude and of the same character as 
that of burned lamps; it is greater during initial 
periods than during later stages. Theories and 
suggestions for the tie-up of mercury are ad- 
vanced. Superficial oxygen ions of the phosphor 
crystals are assumed to function as a bridge, link- 
ing mercury atoms to the bulk of the crystal. 
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An elegant way of studying the geometrical meaning of 
a Harker synthesis is to transform it into a diagram here 
termed the implication diagram. This is done by trans- 
formation of polar coordinates, the characteristics of the 
transformation depending on the number of operations in 
the cyclical group of the axial symmetry for which the 
Harker synthesis is prepared. The implication diagram 
has the important property that it is a map of the location 
of atoms in the crystal structure plus additional locations 
(ambiguities) which would give rise to the same Harker 
synthesis, plus satellitic locations. Fortunately the satellitic 
locations can be identified as such. Satellitic peaks are 
caused by powers of rotation operations and by reflections. 
Satellitic peaks occur in specialized locations which are 
discussed for the several cases. On the other hand, it is 
not generally possible to decide between ambiguities. 
An implication diagram which is a map of the crystal 
structure, with ambiguities, but without satellitic loca- 
tions, exists for the following 18 of the 21 possible 
parallel-axial symmetries: P2, P2,, C2, H3, H3;, H32, R3, 
P4, I4, P4,, Po, P4s, T4,, H6,, H6;, H6., H6, sand H6s3. 
This is also true of certain axial symmetries combined with 
reflections, such as 2¢ and 3c. Only for the three axial 
symmetries P4, /4, and H6 do no implications exist without 


INTRODUCTION 


HE major difficulty besetting crystal struc- 

ture analysis is the loss of phase relations 

in the recording of the diffraction. If these phases 
could be determined, all crystal structures could 
be solved. Patterson? investigated the alterna- 
1A, L. Patterson, “A Fourier series method for the 
determination of interatomic distances in crystals,” Phys. 


Rev. 46, 372-376 (1934). 
2 A. L. Patterson, “‘A direct method for the determina- 
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satellitic peaks. The implication diagrams of the Harker 
syntheses P(xy}) for the space groups H6,, H6;, H6,2, 
and H652 and for Harker synthesis P(xy0) for R3 lack 
both ambiguities and satellitic peaks. For these space 
groups the implication diagrams are true maps of the 
crystal structure. Harker syntheses greatly exaggerate the 
electron densities of atoms in the structure, and for this 
reason their greatest usefulness lies in the location of the 
more compact atoms. Diffuse atoms, and_ therefore 
especially anions, provide Harker peaks of such low mag- 
nitude as to be often lost in the background. Background 
due to non-Harker interactions is discussed. An important 
characteristic of Harker syntheses, which is independent 
of any interpretation of their meaning with regard to the 
location of atoms in the crystal structure, is that they 
provide criteria for distinguishing space groups which 
cannot be distinguished in ordinary qualitative x-ray 
crystallography. This is because the Harker synthesis 
provides definite criteria for identifying reflection planes 
and rotation axes. It fails only to distinguish pairs of space 
groups which differ by a group of inversions alone. Harker 
synthesis thus realizes the ultimate possibilities of x-ray 
crystallography because it is quantitative. 


tive question, namely, what information can be 
ascertained about a crystal structure if the phases 
of the diffraction spectra are unknown. He dis- 
covered the important function, 


P(xys)= 5 SE | Ful? 


h=—oa k=—w |=—@ 
Xcos 2r(hx+kyt+lz), (1) 
tion of the components of interatomic distances in crys- 
tals,”’ Zeits. f. Krist. A90, 517-542 (1935). 
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Fic. 1. Transformation of position of two points related 
by operation A(a) from crystal to Harker diagram and 
return. a: Relations in crystal. b: Relations in Harker 
diagram. c: Relations in implication diagram. 


and proved that this function is equivalent to 
another giving the average value of the product 
of the electron density at the two ends of a vector, 
as the vector, fixed in direction and length, is al- 
lowed to range over the crystal cell. This function 
has a high value when two atoms in the crystal 
structure are separated by this vector. Conse- 
quently (1) has maxima at points whose locations, 
with respect to the origin, define vectors which 
correspond to interatomic vectors within the crys- 
tal. 

Harker*® pointed out that, for crystals having 
any symmetry elements other than inversion 
centers, certain special regions of Patterson’s 
three-dimensional function are the loci of the 
ends of vectors relating symmetry-equivalent 
_atoms. These are therefore very useful in locat- 
ing the positions of such atoms. For the purpose 
of locating atoms related by axes of symmetry, 
certain plane sections of the graphical repre- 
sentation of Patterson’s function are useful, 
‘ while certain lines are useful in aiding in the 
location of atoms related by symmetry planes. 

* David Harker, “The application of the three-dimen- 
sional Patterson method and the crystal structures of 


proustite, AgsAsS;, and pyrargyrite, Ags;SbS;,”" J. Chem. 
Phys. 6, 381-384 (1936). 
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In this paper, the discussion is concerned with 
atoms related by axial symmetry elements, with 
or without symmetry planes parallel with the 
axes. For axial symmetry, Harker recommends 
the particular Patterson functions P(xy0) and 
P(xy}). For these he gives the following form 
(axis assumed parallel to Z). 


Screw Component of Axis =0 


P(xy0)= > DY Ccos2r(hx+ky), (2) 


h=—2% k=-x 


where 


C= > x Frxi| ?. . 
l=—@ 


Screw Component of Axis=c/2 


P(xv3)= © DY Ccos2a(hx+ky) (3) 
h x k=—x 


where 


x 


C= >. (—1)! Fix , 
l=—m 


It will be evident from what follows that these 
two cases, P(xy0) and P(xy}), are not sufficient, 
and are not even recommended in some cases 
where they apply. A more general form for 
P[xy(1/w) ] is 


1 w ra 
p(s») => > Cicos2r(hx+ky) 
h=—a k=—@ 


Ww 


-> > C2 sin 2r(hx+ky), (4) 
h=—a k=—@ 
where 
. wo 2n 
Ci= Y cos—]| Fu\?, 
l=—a Ww 
and 
oa : Qn ; 
C= > sin —l| Fru ?. 
l=—a@ WwW 


THE IMPLICATION DIAGRAM 


The interpretation of Harker diagrams for 
certain axial symmetries is greatly simplified by 
a transformation of the data to another diagram 
which will be termed the implication diagram. 
The transformation is as follows: In Fig. 1a, 
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suppose that the two atoms, 1 and 2, are related 
by an operation A(a), of a group of rotations. 
(The notation here is that ordinarily used in dis- 
cussions of crystal symmetry, namely, the axis 
of the rotation is designated by A and the 
amount of the rotation by a.) The atoms 1 and 
2 determine two interatomic vectors, namely 1,2 
and 2,1. Consequently the Harker diagram, 
Fig. 1b, contains two peaks, 1,2 and 2,1, which 
correspond to these vectors. Thus, to each dis- 
tribution of atoms in the crystal there corre- 
sponds a unique distribution of peaks on the 
Harker diagram. 

Consider the reverse transformation. What 
inference can be drawn from the Harker diagram 
regarding the distribution of atoms in the crys- 
tal? To find the possible atomic positions which 
have given rise to the Harker peaks, all that is 
necessary is to find two points in the crystal 
equidistant from axis A which subtend an angle 
a at A, and which are separated by distance 
1,2=r. There are two solutions to this implica- 
tion, Fig. 1c, unless a=, in which case the two 
solutions coincide. 

The mechanism for transformation of any peak 
from the Harker diagram to the implication 
diagram is easily formulated. If the polar co- 
ordinates of the Harker point are r and ¢ (Fig. 
2a), the polar coordinates of the corresponding 
two implication points (Fig. 2b) are 





sS= = 
Tv 
2 sin — 
n 


a 
2 sin — 


; : Qn 
where 1 is defined as — 


a >(5) 


Tv Qa 
h=0+/ tains ), 

- a 

Tv a 
mois 

2 a ) 


An especially easy way of making use of this re- 
lation for certain axial symmetries will be dis- 
cussed after certain other properties of the impli- 
cation diagram are mentioned. 

From (5) it is evident that each pair of points 
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on the implication diagram corresponding to 
one point on the Harker diagram is separated 
by an angle of (4/2—a/2)+(x/2—a/2)=r—-—a 
at the axis A. The following discussion will show 
that if each of the m points of a symmetrical set 
on a Harker diagram is transformed to the 
implication diagram, the points of the implica- 
tion diagram with azimuth ye in (5) may be 
disregarded. 

Harker syntheses are cosine summations, so 
that for each point, 1, at coordinates x;, y; there 
is another point, p2, at coordinates #1, 9:. The 
point, p2, subtends an angle w at the axis A 
with point p. Because of the symmetry of the 
axis A, there is a third point, 3, which subtends 
an angle +a with p;. Now consider the implica- 
tion points corresponding with p,; and p3. The 
second implication point of p; has an azimuth of 
g—m/2+a/2; the first implication point of ps; 
has an azimuth of ¢+2+a+(2r/2—a/2)=¢ 
+32/2+a/2=9—2/2+a/2. Thus the yo of p; 
is identical with the y of p3. As a consequence, 
if all points of the axially symmetrical set of n 
Harker points are to be transformed into their 
respective implication points, it is sufficient to 
locate on the implication diagram merely one 
point with azimuth y, for each Harker point, 
since the second point with azimuth yz will ap- 
pear as the ¥ point corresponding with another 
point of the symmetrical Harker set. 

It follows from this discussion that to every 
Harker set of axially equivalent points there 
corresponds a symmetrical set of m points on the 
implication diagram. Furthermore, these ” im- 
plication points, except for the ambiguities dis- 
cussed in the following sections, are the posi- 
tions of atoms in the projection of the crystal 
structure on the plane normal to the axis. 

The several kinds of ambiguities involved in 
implication diagrams will now be considered. 


Implication 





Fic. 2. Relations between polar coordinates of points on 
Harker diagram and implication diagram. 
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AMBIGUITIES IN THE IMPLICATION 
Ambiguities Due to Centrosymmetry 


Since, for every interatomic vector there 
exists the reverse vector, all Harker diagrams as 


A’ 
78... 
a . 
_ New A a 
ax - = 
e ’ ™s 
>? . 
-— 
Fic. 3. The combination of an operation of rotation, 


A(a), with an operation of translation, ¢, and the conse- 
quent operation A’(a). 


~ 
~ 





well as all implication diagrams are centrosym- 
metrical. Thus the group of symmetry opera- 
tions of the axis A of the Harker diagram or 
implication diagram must contain the operation 
A(x) whether the symmetry of the crystal axis 
A contains this operation or not. The only 
crystallographic case where it does not is for 
n= 3 (not including the symmetry designated by 
3”). This requires sixfold axes to appear on 
the Harker diagrams and implication diagrams 
corresponding to »=3 in the crystal. Thus for 
n=3 only, a twofold ambiguity arises in the 
interpretation of implication diagrams due to 
this cause. 


Ambiguities Due to Translations 


The interpretation of Harker diagrams and 
implication diagrams is complicated by the trans- 
lation symmetry of crystals. This is funda- 
mentally caused by the well-known fact* (Fig. 
3) that if the operation of a rotation, A(a), is 
combined with a translation ¢ at right angles to 
the axis, the combination is equivalent to an- 
other rotation, A’(a), of identical magnitude 
about a parallel axis, A’. This causes the parallel 
symmetry axes (as well as the axes representing 
their subgroups of operations) of space groups to 
occur in several sets which are not necessarily 
translation-equivalent. These two axes are trans- 
lation equivalent, however, for n=6, where they 
are related by a translation at right angles to 


_ ‘Harold Hilton, Mathematical Crystallography and the 
Theory of Groups of Movements (Clarendon Press, Oxford, 
1903), p. 147. 
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the axis; they may also be equivalent by a trans- 
lation which is not perpendicular to the axis for 
other values of m when they occur in the centered 
and rhombohedral lattices. 

The consequences of this for the interpretation 
of Harker diagrams and implication diagrams 
may be illustrated for the simple case of a one- 
dimensional ‘‘space group” represented by a 
linear array of twofold axes, Fig. 4b. The twofold 
axes occur in sets A and A’, which are not transla- 
tion-equivalent. If an atom is placed at 1, then 
the twofold axis requires another at 2, and the 
translation reproduces these at lf, 2, etc. In 
addition to the interatomic vector 1,2 and its 
reverse, there arises vector 1/, 2, and its reverse, 
due to translation. The Harker diagram, Fig. 4c, 
then displays peaks at the ends of these two 
vectors, both of which radiate from the same 
origin. This kind of multiplication of peaks does 
not give rise to any ambiguity, for in transform- 
ing the Harker diagram to its implication, one 
need only consider any section of total length ¢ 
(most conveniently the section within the range 
+t/2 of the origin). Only vectors corresponding 
to points within this range are distinctive, and 


— 
x A A 





“ A = 


Fic. 4. Relations between crystal, Harker diagram, and 
implication diagram, and the cause of ambiguity. Locations 
of symmetry elements on implication diagram, which are 
not present in crystal, are indicated by outlined symbols. 
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one therefore disregards vectors corresponding 
to points beyond this range. 

There is, however, a second consequence of 
principle illustrated in Fig. 3. If the same dis- 
tribution of atoms with respect to axis A, shown 
in Fig. 4b, occurs about axis A’, as shown in 
Fig. 4a, the same Harker distribution as that 
shown in Fig. 4c results. It is therefore impossible 
to determine whether the Harker diagram im- 
plies Fig. 4a or Fig. 4b. This is a general feature 
of implication diagrams, and it arises because 
the space group, in consequence of the principle 
illustrated in Fig. 3, contains several sets of non- 
equivalent axes of equal n. In space groups, the 
number of sets of such parallel axes depends on 
the symmetry of the axis. For uncentered 
lattices there are always four sets of twofold 
axes, three sets of threefold axes, two sets of 
fourfold axes, and only one set of sixfold axes. 
This number is reduced in centered and rhombo- 
hedral lattices by the factors 2 and 3, respect- 
ively. Let this number of sets be represented by 
m. Then, in transferring a Harker synthesis to 
the implication diagram, an m-fold ambiguity 
is involved, since one cannot tell which of the 
m sets of axes in the crystal relates the atoms so 
as to give rise to the Harker points. 

In order to take care of this ambiguity, the 
implication diagram requires that the inferred 
locations of atoms shown in Fig. 1c be plotted 
about each of the m non-equivalent axes of the 
space group, as illustrated in Fig. 4d for this 


TABLE I. 


Ambiguity of 
location, m 
Centered 
or 





Harker Primi- rhom- 
Axial syn- tive bohedral Subgroups of 

n symmetry thesis lattice lattice same n 
2 2 P(xy0) 4 2 4?, 4.7, 6.3, 6,3 

2 P(xy}) 4 2 4,2, 43%, 613, 65° 
3 3,3(=3+3), P(xy0) 3 1 63? 

6( =3/m) 

31, 3e P(xy}) 3 2 6:’, 6;?, 627, 62 
4 4,4 P(xy0) 2 1 

4, 43 P(xy}) 2 2 

4o P(xy}) 2 2 
6 6 P(xy0) 1 — 

61, 65 P(xy$) 1 — 

62, 4 P(xy}) 1 — 

63 P(xy}) 1 ‘and 
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TABLE II 
n 

| 2 3 4 6 

| 1 1 v2 
Ss | mM ww 3° r 
y | @ g+30° e+45° ¢+60° 
Shrinkage factor | ; S v 1 
Rotation 0° 30° 45° 60° =0° 
simplified case. This has interesting conse- 


quences. Since identical implications are plotted 
about non-equivalent axes, the latter take on 
the characteristics of translation equivalence in 
the implication diagram. This diagram therefore 
always has shorter translation periods than the 
crystal (except in the case of m=1). 

Table I succinctly summarizes some important 
properties of the Harker diagrams and implica- 
tion diagrams for the several axial symmetries. 
Unfortunately, for n>3, one cannot deal with 
sets of parallel axes of equal m. This is because 
the powers of the operations of such axes form 
subgroups of lower n. In the last column of Table 
I a superscript, r, appended to the symbol of a 
symmetry axis, designates a subgroup composed 
of the rth, 2 rth, 3 rth, etc., powers of the funda- 
mental operation of the axis. All subgroups on 
the same line are identical in character with the 
group indicated by the symbol in column 2. 
The arrangement is intended to bring out the 
fact that for purposes of constructing and in- 
terpreting Harker and implication diagrams, 
atoms related by a 6; axis, for example, can be 
studied by the Harker synthesis P(xy$); or by 
the synthesis P(xy}), thus studying, in implica- 
tion, atoms related by the second power of the 
operation of the screw; or by the synthesis 
P(xy}), thus studying, in implication, atoms re- 
lated by the third power of the operation of the 
screw. Harker’s original paper recommends the 
last two cases. The study should always be made 
by P(xy%), however, for the ambiguity factors are 
very different, namely 1, 3, and 4 for the powers 1, 
2, and 3, respectively. When the implication of 
P(xy%) is prepared, there is no ambiguity due to 
translation, and the implication diagram is a map 
of the crystal structure projected on the plane 
normal to the sixfold axis. 
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Fic. 5. Relations between the symmetry elements of the crystal, Harker diagram, and implication diagram for the 
four possible values of n. The outlines of the crystal cells are indicated by broken lines. The figure shows the domains 
(outlines in full lines) of the Harker synthesis and their transfer to the implication diagram for each case. The locations 
of new symmetry in the implication diagram for the case of primitive cells are indicated by outlined symbols. 
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CONSTRUCTION OF IMPLICATION DIAGRAMS 


It is evident that an implication diagram can 
be constructed by taking the desired features, 
such as peaks, of the Harker diagram and trans- 
forming these, point by point, to an implication 
diagram. A better way of handling the manipula- 
tion, however, will be evident from the following 
discussion. : 

It will be observed from (5) that every point 
on the Harker diagram transforms into a point 
on the implication diagram by a shrinkage of its 
radial polar coordinate and a rotation of its 
azimuth. The shrinkage and rotation are con- 
stant for a given value of the axial symmetry 
parameter, m. The specific values of shrinkage 
and rotation for various values of m are given in 
Table II. 

Evidently an implication diagram can be 
plotted from the Harker synthesis by transforma- 
tion of the polar coordinates of the synthesis. A 
practical way of performing this is to plot the 
Harker synthesis to a scale reduced by the cor- 
rect radial factor of Table II, rotate the plot by 
the correct angle given in Table II, and lay 
down the resulting plot on the m locations of 
similar axes discussed under ‘‘Ambiguities Due 
to Translation.’”’ The implication diagram con- 
structed in this way shows all the detail of the 
Harker synthesis, instead of just the peaks. 

No interpretation need take place up to this 
stage. Any proposed crystal structure can then 
be plotted directly on the features of the implica- 
tion diagram. This procedure has the advantage 
of deferring any decisions about the meaning 
of doubtful features of the Harker diagram until 
they can be tested directly in comparison with 
the features of the proposed crystal structure. 

Figure 5 shows the relations between sym- 
metry elements in the crystal, the Harker dia- 
gram, and the implication diagram for the four 
possible values of m for pure axial symmetry. 
This figure also shows domains in the Harker 
synthesis and their transfer to the implication 
diagram for the several axial symmetries. Sym- 
metry elements in the implication diagram shown 
in outline represent locations of symmetry not 
present in the crystal. Cell outlines are shown in 
broken lines. 
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SATELLITIC PEAKS ON IMPLICATION DIAGRAMS 


Satellitic Peaks Due to Powers of 
Rotation Operations 


The groups of operations represented by axial 
symmetries (pure rotations, screws, and roto- 
reflections) are cyclical. Hence only powers of the 
fundamental operation are members of the group. 
This gives rise to important consequences. In 
the case of screw symmetry, each atom can only 
be related by any selected operation of the screw 
to two other atoms, namely, one atom above 
and one atom below. There are thus only two 
Harker vectors from the atom due to this oper- 
ation and its reverse. The two vectors have equal 
length and diverge by an angle of r —a=a—2n/n. 
For the case of m = 2, the two vectors coincide, but 
they are distinct for all other cases. Thus, unless 
further complications occur, implication diagrams 
for screws are comparatively simple, and indeed, 
each peak on the implication diagram is the 
location of an atom, except for ambiguities. 

No further complication occurs for symmetries 
with »=2 or 3, nor for the screws 4;, 42, 43, 61, 
and 65. Unfortunately, however, complications 
do arise for the axial symmetries 4, 6, 62, 63, and 
6,. This is due to the fact that for these five 
symmetries (Fig. 6) there occur in the same 
plane not only the atom related to the original 
atom by the first power of the symmetry opera- 
tion, but also one or more additional atoms re- 
lated to the original atom by higher powers of 
the same symmetry operation (possibly reduced 
by a translation parallel to the axis). This gives 
rise to Harker vectors in addition to the two 
mentioned in the last paragraph. The several 
cases are illustrated in Fig. 6. 

The additional vectors complicate the in- 
terpretation of Harker syntheses. The implica- 
tion diagrams not only contain peaks at location 
of atoms, plus their ambiguities, but also addi- 
tional peaks where there are no atoms. Of course, 
implication diagrams for such symmetries are 
easily interpreted if the crystal contains only one 
set of equivalent atoms, an extremely rare case 
indeed. In the more usual case of a crystal con- 
taining several sets of atoms, the diagram is 
more difficult to interpret. The interpretation is 
aided by taking account of the relations between 
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the multiple sets of points, both with regard to 
geometry and comparative intensities. 

The confusion due to the multiple sets of 
points is mitigated, for the symmetries 62, 6s, 
and 6,4, by the fact that Harker syntheses can 
be prepared for these symmetries for powers 
other than the first power. Thus, the symmetry 
6; can be treated as 6;° by using P(xy0). The 
implication diagram then has the characteristics 
of n=3. Similarly, the symmetries 62 and 6, can 
be treated as 6° and 6, by using P(xy0). The 
implication diagrams then have the character- 
istics of n=2. Each of these treatments involves 
the ambiguities characteristic of its nm. It is 
possible to make use of this change in character 
of the ambiguity so as to provide a criterion for 
discarding satellitic peaks. This can be done by 
preparing implication diagrams for the several 
different powers of the same operation. Each 
implication contains the true implication plus a 
false implication due to ambiguities plus satellitic 
peaks. The peaks which these diagrams have in 
common are the true locations of the atoms, the 
false ones being, in general, different on the 
implication diagrams prepared for two different 
powers of the operation. 

The characteristics of the possible implication 
diagrams can be formulated by considering the 
implication diagrams of the possible axial space 
groups instead of merely the individual axes. 
There are 21 space groups involving only opera- 
tions about sets of axes parallel to a single line. 
The characteristics of the implication diagrams 
of these groups are set forth in Table I 11. The first 
column of this table lists the Hermann-Mauguin 
symbol of each of the 21 parallel-axial space 
groups. The second column provies the possible 
Harker sections for each such space group. The 
third column gives the possible values of » for 
which the implication can be constructed. For 
each of these values the axial symmetries and 
their ambiguity factor are listed in the next two 
columns. In the last column are shown the sub- 
groups (if any) of the space group which give 
satellitic peaks in the implication diagram for 
that implication construction. 

The table brings out the point that for several 
parallel axial space groups there exist implica- 
tion diagrams which have neither translation 
ambiguity (i.e., 2=1) nor satellitic peaks (no 
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interfering subgroups). These are R3, H6,, and 
H6;. The perfection of R3 in this respect is 
marred by the fact that there is an ambiguity 
due to centrosymmetry. The group R3, however, 
is perfect in this respect. Evidently all axial 
space groups having R3, H6;, or H65 as sub- 
groups have implication diagrams whose peaks 
are located at the positions of the atoms in the 
crystal structure. Thus implication diagrams 
exist for the following space groups which are 
maps of the crystal structure: R3, H6;, H6,2, 
H6;, H6;2. The following axial space groups have 
trigonal implication diagrams which are maps of 
the crystal structure plus a centrosymmetrical 











TABLE III. 
Subgroups 
causing 
Axial Harker satellitic 
groups synthesis n axis m peaks 
P2 P(xy0) 2 2 4 
P2, P(xy}) 2 21 4 
C2 P(xy0) 2 2 2 
P(xy}) ; +t 2 
H3 P(xy0) 3 3 3 
H3, P(xy}) 3 31 3 
H32 P(xy}4) 3 3 3 
R3 P(xy0) 3 3 1 
_ P(xy}) 3 31 2 
Pa P(xy0) 2 42 =2 4 
13 P(xy0) 2 q2 2 
P(xy}) 2 21 2 
P4 P(xy0) . 4 2 P2 
\2 2 4 P4 
PA, P(xy}) 4 4, 2 
P(xy}$) 2 4? =2, 4 
P4; P(xy}) 4 4, 2 
P(xy}) 2 47=2, 4 
P4s, P(xy}) 4 2 2 
P(xy0) 2 4.2 =2 4 
14 P(xy0) 4 4 1 12 
2 4.2=2 2 I4 
P(xy}) {4 4, 1 I2 
\2 2 2 14 
T4, P(xy}) 4 4; 2 
P(xy}) 2 42=2) 2 
P(xy0) 2 2 2 
H6 P(xy0) (6 6 1 3. P2 
43 3 3 P2, H6 
(2 2 4 H6, H3 
H6, P(xyi ) 6 6; 1 
P(xy}) 3 6; =3, 3 
P(xy}) 2 6;3 = 21 4 
H6; P(xy?) 6 65 1 
P(xy}) 3 6;? = 32 3 
P(xy}) 2 65 = 2; 4 
H6. P(xy}) (6 62 1 H32 
\3 62? =3. 3 H62 
P(xy0) 2 623 =2 4 
H6, P(xy}) {$ 64 1 H3, 
3 63 = 3, 3 H6, 
P(xy0) 2 63 =2 4 
H6; P(xy}) (6 6s 1 P2, 
\2 2, 4 H6; 
3 3 


P(xy0) 
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Fic. 7. The appearance of satellitic peaks due to reflection symmetry. The relative 
heights of the peaks are indicated by small numerals on the implication diagram. 


ghost: R3, R32, and all the axial isometric space 
groups, namely P23, F23, 123, P2,3, 12:3, P43, 
P4.3, F43, F4,3, 143, P4;:3, P43, and 14,3. 
Further useful information of this character is 
assembled in Table IV. 

The conclusions of this section can be sum- 
marized as follows: For 18 of the 21 possible 
parallel-axial symmetries (i.e., all symmetries 
except P4, /4, and H6), there exists an implica- 
tion diagram which is a map of the crystal struc- 
ture repeated m(=4, 3, 2, or 1) times by transla- 
tion (and, in the case of trigonal crystals lacking 
a center of symmetry, repeated again by in- 
version). For each axial space group containing 
one of these parallel-axial space groups as sub- 
groups, therefore, there exists an implication 
diagram which is a map of the crystal structure 
plus m—1 ghosts (all repeated by inversion in 
the case of trigonal crystals lacking a center of 
symmetry). For the three axial symmetries P4, 
I4, and H6, all possible implication diagrams 
contain, in addition, one or more times this 
number of satellitic peaks due to subgroup 
operations. 


Satellitic Peaks Due to Reflection Symmetries 


The addition of reflection symmetry in planes 
normal to the symmetry axes does not give rise 
_ to any additional peaks. The addition of sym- 
metry planes parallel to the symmetry axes may 
or may not give rise to satellitic peaks on implica- 
tion diagrams. For purposes of discussing this, 
suppose that all symmetry axes are regarded as 
screws (the screw component of which may be 
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zero in the special case of a pure rotation axis) 
and that all symmetry planes are regarded as 
glide planes (the translation of which may be 
zero in the special case of a pure reflection plane). 
Suppose that the Harker synthesis, 


1 
P ( xy— ) , 
WwW 


is made for a translation component, 1/w, of the 
screw. Then, if the symmetry plane does not 
also have a translation component, 1/w, it does 
not provide additional atoms separated by this 
level difference. Hence it does not give rise to 
new interatomic vectors of this separation, and 


TABLE IV. 


Symmetries without satellitic peaks 


Translation 


Total Centro- ambiguity 
ambiguity symmetry factor, 

factor factor m Parallel-axial groups 
1 1 1 R3, H6,, H6s 
2 z 1 R3 _ 
2 1 2 C2, 14, P41, P43, P42, [4 
3 1 3 H3, H3;, H32, H63*t 
4 1 4 


P2, P2,, P4,t H62,t H6st 
Symmetries with satellitic peaks 


Number of Translation 


sets of ambiguity 
satellitic factor, Parallel-axial 
peaks m groups 
1 1 14, H6;* 
1 2 P4 
2 H6.,* H6,* 
3 1 H6 


* Subgroup is without satellitic peaks. 
Tt Subgroup only. 
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Fic. 8. The relations between peaks and satellitic peaks for the symmetries 2mm, 3m, 4mm, and 6mm. 
The relative heights of the several peaks, where these are different, are indicated by small numerals on 
the implication diagram. Points and symmetry not present in the crystal are indicated by outlined symbols. 
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consequently no new peaks appear on the Harker 


synthesis, 
1 
P(x | ). 
Ww 


In this case the implication diagram is a map of 
the crystal structure, with the exception of the 
ambiguities and satellitic peaks already dis- 
cussed. Important cases of this are two-, three-, 
four-, and sixfold axes combined with parallel 
glide planes whose glide component is c/2. 

If, on the other hand, the translation com- 
ponent of the screw axis and the axial component 
of the glide plane are the same, the reflection 
symmetry introduces additional atoms at the 
same level, and the Harker syntheses contains 
peaks representing vectors from each atom both 
to its rotation and its reflection equivalent. It 
is convenient to divide the discussion of such 
instances into two cases, according as the plane 
contains the axis or is parallel with the axis. 


Symmetry Plane Containing Symmetry Axis 


Figure 7 shows the simplest example of this 
case, namely a rotation and reflection combina- 
tion in which only one plane has the same transla- 
tion component as the axis. This example may 
serve as a basis for the development of more 
complicated examples. The solid circles and 
double circles each constitute sets which, in 
themselves, are related by the twofold screw 
only. The interatomic vectors for these sets 
within themselves are the usual ones, and are 
shown in the general positions of the correspond- 
ing Harker and implication diagrams. There is 
only one other kind of interatomic vector intro- 
duced by glide plane, namely a vector at right 
angles to the plane; this is shown in the lower 
part of the crystal diagram. Such a vector always 
gives a point in a specialized location on the 
‘Harker and implication diagrams, namely on the 
central line at right angles to the plane. Its 
coordinates are 2x,0 on the Harker diagram for 
the situation shown, but x,0 on the implication 
. diagram. Since the relation of scale on implica- 
tion to Harker varies with n, it is better to ob- 
serve the general geometry of the location of 
this sateliitic point with respect to the true impli- 
cation points than merely the coordinates of the 
satellitic points. 
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In Fig. 8 are shown the relations between the 
atoms in the crystal, the Harker diagram, and 
the implication diagram for the important cases 
of reflection planes containing the several n-fold 
rotation axes. The Harker diagram in each case 
is derived in a manner similar to the simple 
case just treated: The atoms are first regarded 
as divided into two axially symmetrical sets. 
The two sets are related to one another by the 
several reflection symmetries of the planes. The 
total number of interatomic vectors in the two 
sets consists of (a) the vectors within each set 
due to axial symmetry, plus, (6) the vectors 
from one set to the other due to each reflection 
plane. The first collection of vectors, (a), lo- 
cates, in the implication diagram, just the posi- 
tions of the atoms in the structure (plus attend- 
ant ambiguities characteristics of m, and satellitic 
peaks previously discussed). The second set of 
vectors all lie in special positions, namely on 
lines in the Harker diagram at right angles to 
the plane producing the vector in the crystal. 
In the implication diagram, this line is rotated 
out of normality with the plane producing the 
points by an amount characteristic of the number 
n (see Table 11). This rotation places the line 
normal to another plane of the centrosym- 
metrical group. The peaks at the ends of these 
vectors on the implication diagram do not repre- 
sent possible positions of atoms. They can be 
recognized by their geometrical relation to the 
peaks in the general position. In the diagrams of 
Fig. 8, the derivation of the satellitic peaks on 
the Harker diagram due to a single reflection 
plane of the symmetrical set is given. Vectors 
represented by solid lines are those due to re- 
flection. Vectors represented by broken lines are 
certain ones-due to axial symmetry which are 
included to show the relation of rotation vec- 
tors to reflection vectors in producing the ge- 
ometry of the pattern of true and _satellitic 
peaks. Note that for »>3, both axial and re- 
flection vectors occur for powers of operations. 
These are shown in bracketed sets in the left 
column of Fig. 8. 

The method of analysis of diagrams having 
axes with intersecting reflections is indicated in 
the right column of Fig. 8. To relate satellitic 
peaks to non-satellitic peaks caused by the same 
symmetrical set of atoms, first locate the non- 
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Fic. 9. The possible combinations of twofold screws with parallel reflections, and their implication 
diagrams. Points and axial symmetry in the implication diagram which have no.counterparts in the crystal 
are shown by outlined symbols. 
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satellitic peak. Unless the atom is in a special 
position on one of the reflection planes, this is 
the peak in the general position closest to the 
origin. From this peak, drop a pair of lines to 
each of the two nearest symmetry planes of the 
crystal (for »=3, ignore the spurious extra 
planes introduced by centrosymmetry; they can 
be identified by not containing satellitic peaks) 
and making angles with the plane of 90° for 
n=2, +60° for n=3, +45° for n=4, or +30° for 
n=6. (For n=2, these two lines coincide, causing 
a satellitic peak of double height.) At the points 
of intersection of the two lines with each sym- 
metry plane lie the two satellitic peaks due to 
reflection for n=2, 3, 4, or 6, respectively. 

This completes the process for »=2 or 3, but 
for n=4 or 6, there are further satellitic peaks 
due to powers of the operation and their at- 
tendant reflections. The following directions for 
analyzing these additional peaks rest on the 
fact that (left column of Fig. 8) for the fourfold 
case, a reflection vector is common to the two 
zones n=4 and n=2, while for the sixfold case, 
a reflection vector is common to the two zones 
n=6, n=3, and another is common to the two 
zones n= 3, n=2. Thus the outer reflection peaks 
just found are also related to the required axially 
satellitic peaks. For the fourfold case, draw a line 
at right angles to each of the two non-equivalent 
reflection planes at the positions of the outer 
satellitic peaks, found as described above. These 
two lines intersect each other in the axially 
satellitic peak. This completes the analysis for 
the fourfold case. For the sixfold case, the first 
subgroup has »=3. To find the axially satellitic 
peak corresponding ‘with this, drop lines from 
the two outer reflection-satellitic peaks (one on 
each plane) as found above, each line making 
an angle of 60° with its plane, the pair converg- 
ing outward. These intersect in an axially satel- 
litic peak due to the subgroup with n=3. Com- 
plete the »=3 diamonds, thus locating the »=2 
reflection-satellitic peaks at their outer points. 
At these draw perpendiculars from the planes. 
- Their intersections are the n=2 axially satellitic 
peaks. This completes the process for the sixfold 
case. 

In this manner, starting with a peak in a gen- 
eral position, all satellitic peaks caused by the 
same symmetrical set of atoms can be located. 
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Thus located, satellitic peaks can be discarded 
and only non-satellitic peaks need be considered 
in the further analysis of the implication diagram. 

By giving attention to the geometrical rela- 
tions between rotation and reflection peaks on 
implication diagrams, such diagrams can thus be 
easily interpreted for cases where the axis con- 
tains the plane. The only theoretical difficulty 
arises when the atoms responsible for the peaks 
are not in the general position but rather are 
located on the specialized position on the line 
through the axis and normal to the plane. When 
this occurs, the peak due to such an atom is 
correctly located on the implication diagram, 
but the location superficially suggests that the 
peak is satellitic. It can be distinguished from 
the reflection-satellitic peaks just discussed be- 
cause it is not geometrically related (in the man- 
ner discussed) either to (a) peaks in the general 
position, or (b) to other peaks on its own sym- 
metry line. 


Symmetry Plane Parallel with Axis 


Even when the symmetry plane does not con- 
tain the axis, the positions of reflection-satellitic 
peaks are confined to specialized locations so 
that these can be recognized. Several simple 
cases involving twofold screws and representa- 
tive kinds of reflections are shown in Fig. 9. 
Two different kinds of cases can be recognized. 
In general, the atoms of a set a1, d2, ds . . ., are 
related by the operations of rotation, and this 
set is related to another set, Gis, Gos, G3. . . ., DY 
the operation of reflection. If the Harker syn- 
thesis is carried out (symmetry permitting) so 
that vectors 1—1is, 2—2s, 3—3s ... result, 
then the reflection-satellitic peaks are confined 
to a line normal to the reflection plane (Fig. 9a 
and 9b), a situation similar to that resulting 
when the reflection plane contains the axis, as 
discussed in the last section. 

A slightly different situation occurs when the 
Harker synthesis is carried out so that vectors 
other than 1—1s, 2—2s ... result. Figs. 9c 
and 9d illustrate cases containing twofold axes 
where vectors 1—2s, 2—1s result. This general 
situation can be analyzed as follows: Let the 
operation of rotation relate a; to a;, of which the 
corresponding atoms by reflection in the plane 
are @;, and aj. Let the Harker synthesis be car- 
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ried out to that vector j—1s results. Now a; is 
related to a; by a combination of a rotation 
(to a:) followed by a reflection (to a:,). This 
combination is equivalent to a reflection in some 
other plane of the space group making an angle 
with the first of }3(a,/Aa;). Thus the vector 
j—1s is normal to this plane, and the reflection- 
satellitic peak again appears in the implication 
diagram on a line normal to some reflection plane 
in the space group. 

In both the above instances, if the glide of the 
reflection operation has a component normal to 
the axis of amount g (Fig. 9b and 9d), then the 
reflection-satellitic peaks are still on a line normal 
to the reflection plane, but are displaced by an 
amount g/2 from the axis. 


INTERPRETATION OF IMPLICATION DIAGRAMS 


The conclusions arrived at in the foregoing 
discussions may be summarized in more compact 
form as follows: Implication diagrams may be 
prepared for any space group involving axial 
symmetries with or without additional reflec- 
tion symmetries. The diagram is a map of the 
crystal structure, usually, but not always, com- 
plicated by the presence of peaks not correspond- 
ing with any atoms in the crystal structure- 
These unwanted peaks may be divided into two 
categories, namely ambiguities and _ satellites. 

The ambiguities are, in effect, ghost images of 
the crystal structure repeated by inversions (in 
the case of threefold axes only), and translations 
(not the translations of the space group). The 
second kind of ambiguity is a consequence of the 
translational symmetry of crystals. The number 
of ghosts and their placement, arising from this 
cause, is a function of the number of operations 
in the axial symmetry group. It is usually im- 
possible to resolve this ambiguity. 

Satellitic peaks differ from ambiguities in that 
they are not ghosts. They invariably occur in 
positions where they can be recognized and re- 
jected. They are of two kinds: rotation-satellitic 
peaks and reflection-satellitic peaks. The rota- 
tion-satellitic peaks occur in a special pattern, 
geometrically related to non-satellitic peaks. 
The reflection-satellitic peaks are concentrated 
in lines normal to reflection planes, and also 
occur in recognizable patterns. 

Satellitic peaks are missing in space groups 
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characterized by purely axial symmetries of cer- 
tain sorts. They are absent when the Harker 
synthesis can be made for the following axial 
symmetries: 2, 2;, 3, 31, 32, 3=3+2, 41, 42, 4s, 
4, 61, 65, and 6=3/m. They are also absent for 
certain Harker syntheses of 62, 63, and 64. They 
are also absent for certain axial symmetries 
combined. with glides having translation com- 
ponents parallel with the axes. For all such cases, 
the interpretation of implication diagrams is 
comparatively simple, for the implication dia- 
gram is the crystal structure plus m—1 ghosts. 
(The value of m is given in Table I.) 

The value of m is unity for sixfold axes. As a 
consequence of this, space groups containing the 
two kinds of sixfold screws listed in the last para- 
graph are without ghosts. For this reason, the 
implication diagram is a map of the crystal 
structure for the particular space groups H6,, 
H6;, H6,2, and H652 (C, C., Dé’, and Dé, 
respectively). This is also true of R3, except that 
the implication of this group has a Harker peak 
at the origin. This opens the possibility of solv- 
ing the structures of crystals belonging to such 
space groups, no matter how large the cell. 

In addition to the ambiguities and satellitic 
peaks, each Harker synthesis for P(xy0) has a 
satellitic peak at the origin. This is not true for 


syntheses 
1 
P(x) 
Ww 


where 1/w is different from zero. For this reason, 
other things being equal, the latter synthesis is 
to be desired. If it is desirable for other reasons 
to utilize P(xy0), then the additional synthesis 


1 
P(sy~) 


should also be carried out if possible merely to 
locate a possible atom at the origin. 


SPACE GROUP DETERMINATION BY 
HARKER SYNTHESIS 
It is generally recognized that ordinary x-ray 
diffraction studies do not provide criteria for 
distinguishing the 230 space groups from one 
another in all cases. Rather, one is able to dis- 
tinguish 121 diffraction groups’ from one an- 





‘M. J. Buerger, X-Ray Crystallography (John Wiley 
and Sons, New York, 1942), p. 510-511. 


593 








other. Where a diffraction group contains only 
one space group, then that space group is 
uniquely characterized, otherwise not. This is 
fundamentally because such methods depend on 
extinctions, which are able to detect symmetries 
with translation components only. Thus it is 
impossible to distinguish between 1, 2, and m, 
between 4 and 4, or between H3m1 and H31m. 
For example, within the diffraction group 
mmmP- — —, the same extinction effects and 
other qualitative diffraction effects are given by 
P222, Pmm2, and P2/m 2/m 2/m. 

Harker synthesis may be described as quanti- 
tative x-ray crystallography, since by its aid it 
is possible to distinguish all space groups except 
pairs differing only by an inversion. In other 
words, by means of Harker syntheses, it is 
possible to push x-ray crystallography to the 
limits set by the fact that diffraction effects are 
centrosymmetrical. This is because a Harker 
synthesis provides unique characterizations of 
the symmetry elements m, 2, and 4 vs. 4, as well 
as of other symmetry elements. This is an out- 
come of the fact that x-ray data is utilized in a 
quantitative way as against the qualitative cri- 
terion based on the mere presence or absence of 
reflections used in ordinary x-ray crystallography. 
Harker synthesis provides these space group cri- 
teria without any attempt being made to interpret 
the synthesis for the locations of the atoms. 

To give an example of the way the synthesis 
can be used in space group identification, sup- 
pose one has a crystal belonging to the diffrac- 
tion group mmmP —- — -. The possible space 
groups are P222, Pmm2, and P2/m 2/m 2/m. The 
distinction between these can be settled on the 
basis of reflection-satellitic peaks. Representing 
the presence of satellitic peaks by S and their 
absence by 0, the Harker syntheses for the sev- 
eral axes give the following combinations of 
- satellitic peaks: 


P222 000 
Pmm2 SS2 
P2/m 2/m 2/m SSS 


It is even easier to distinguish between such 
cases as H3m1 and H3im. The distinction be- 
tween these two depends on whether the sym- 
metry planes are parallel to the long diagonal or 
the short diagonal of the cell. In Fig. 8, the right 
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column shows that the threefold case has sym- 
metry planes parallel to both directions, but the 
true symmetry plane is marked on the implication 
diagram by containing the satellitic peaks. The 
false symmetry planes lack the satellitic peaks. 
(On the Harker diagram the situation is re- 
versed.) Since the long and short diagonals are 
easily distinguished on the implication diagram, 
this analysis quickly fixes the space group. 


NON-HARKER BACKGROUND 


Although Harker pointed out® that certain 
one-dimensional and two-dimensional Patterson 
syntheses could be used to find interatomic vec- 
tors between atoms not related by symmetry, it 
will nevertheless be convenient to adopt the 
following nomenclature: 


Harker peak: a peak corresponding to an interatomic 
vector between symmetry-equivalent atoms. 

Non-Harker peak: a maximum corresponding to a vector 
between atoms or parts of atoms which are not 
symmetry-equivalent. 


One of the confusing aspects of Harker and im- 
plication diagrams is that they contain a non- 
Harker background. This is always present in 
varying degree. Its characteristics may be such 
as not to seriously affect the interpretation of the 
diagram, or, on the other hand, it may furnish 
peaks of magnitude equal to, or greater than, the 
Harker peaks. Since one does not usually know 
at the outset of the analysis of a crystal structure 
which of these cases he is dealing with, it is 
dangerous to assume that every peak occuring 
on a Harker or implication diagram is a Harker 
peak. This point is emphasized here because 
there appears to be a prevailing opinion that any 
peak on a Harker diagram is a Harker peak; or, 
if non-Harker peaks are recognized at all, that 
at least all strong peaks are Harker peaks. 
Non-Harker background would be absent, ex- 
cept for special cases such as mentioned beyond, 
if atoms were composed of matter concentrated 
at discrete points. The background arises from 
the circumstance that the electronic matter of 
atoms is distributed over a comparatively large 
volume. For pairs of symmetrically equivalent 
atoms, every point in one atom is separated from 
a corresponding point in the other atom by the 
translation:component of the symmetry opera- 
tion. But there are usually parts of the volumes 


JOURNAL OF APPLIED PHYSICS 














of pairs of atoms not related by symmetry which 
are separated by the same translation. The in- 
tegration of the Patterson function extends over 
the entire height of symmetry-equivalent atoms, 
but only over as much of the height of non- 
equivalent atoms as overlap by the translation 
component of the symmetry element. This effect, 
of course, usually causes the peak due to the in- 
teraction between non-equivalent atoms to be 
comparatively weak compared with that due to 
equivalent atoms. The difference is augmented 
by the fact that the dense centers of equivalent 
atoms interact, while in the case of a pair of 
non-equivalent atoms either the diffuse periph- 
eries interact, or at best, the dense center of one 
interacts with the diffuse periphery of the other, 
unless the atoms are fortuitously separated by 
exactly the translation component of the sym- 
metry element. 

Yet non-equivalent atoms can give rise to 
strong peaks. For example, suppose that a very 
heavy atom is on the symmetry axis. Its Harker 
representation may be a peak at the origin which 
is not distinct from the origin peak. The Harker 
diagram selected for this axis will contain strong 
non-Harker peaks if any other atom has a co- 
ordinate differing from that of the heavy atom 
by the translation component, 1/w, of the syn- 
thesis, and the peak will also be strong if its 
coordinate has approximately this value. 

Furthermore, in simple structures, or in struc- 
tures derived by substituting different atoms in 
simple structures (substitution structures), pairs 
of atoms not symmetrically related may occupy 
levels differing by 0, 3, 3, 3, or 3, which are the 
level differences utilized for finding interatomic 
vectors between symmetry-equivalent atoms for 
various axes. In these instances, intense non- 
Harker peaks appear on the Harker diagrams. 
While it is not generally possible to distinguish 
between Harker and non-Harker peaks, the 
latter should be suspected in crystals of simple 
structure and in crystals expected to show sub- 
stitution relation to simpler structures, for ex- 
ample superstructure crystals. In syntheses in- 
volving satellitic peaks, the non-Harker peaks 
are not accompanied by satellitic peaks located 
in the same way as those to be expected for the 
particular implication. In these cases, non- 
Harker peaks can be recognized. 
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HEIGHTS OF HARKER PEAKS 


It will be recalled that the value of the three- 
dimensional Patterson function at a _ point, 
P(xyz), is the average value of the product of 
the electron densities at the ends of a vector 
defined by a line extending from the origin to 
the point xyz, as the vector, fixed in direction 
and length, is allowed to range over the volume 
of the cell. For a pair of atoms, A, related by 
symmetry, the two atoms at opposite ends of the 
vector, are the same, and the height, H, of the 
Harker peak, above any non-Harker background, 
is given by 

Ha = (pa? )cett wv (6) 


where (pa*)ccit w is the average squared density 
of the atom A, the average being taken over 
the cell. Now, the average over one volume 
differs from the average value over another 
volume merely by a scale constant. It is possible 
to determine (p4”), 4 the average squared density 
of the atom over an arbitrary volume, v, provided 
only that the arbitrary volume is greater than 
that of any atom which might come up for con- 
sideration. For any cell, the peak height, Ha, 
due to the Harker interaction between a pair of 
identical atoms can then be expressed 


H,sa=k(pa?)» Ave (7) 


The proportionality constant holds for all atoms 
in the cell, but changes when a new cell is con- 
sidered. The value (p4?)» » is characteristic of an 
atom, and the interpretation of Harker diagrams 
would be aided by a list of such values for the 
several chemical atoms. 

The function (7) has a high value for compact 
atoms and a small value for diffuse atoms. Con- 
sequently, Harker peaks tend to be high for the 
compact cations and low for the diffuse anions. 
Thus, for the silicates, the silicon and oxygen 
atoms, which are approximately isoelectronic, 
provide Harker peaks between which there is 
great contrast. Indeed the oxygen peaks are 
swamped by background. Peak contrasts in 
such cases are even more pronounced than the 
peak heights would lead one to believe, for the 
background contrast is enhanced by the fact 
that atoms with large volumes have peaks which 
decline gradually over larger areas than peaks 
due to small atoms. 
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Dissociation Energies of Surface Films of Various Oxides as Determined by 
Emission Measurements of Oxide Coated Cathodes* 


HAROLD JACOBS 
New York University, New York, New York** 
(Received March 11, 1946) 


A mechanism of falling emission in vacuum tubes is discussed. It is demonstrated that when 
an electron achieves a critical kinetic energy in moving from the cathode to the anode, and if 
the anode is an oxide, the electrons will cause a dissociation of the oxide. The liberated oxygen 
will return to the cathode and falling emission will result. The critical energy of the electron 
starting the dissociation is found to be equivalent to the heats of formation of the oxides bom- 
barded in the case of five different oxides. Using this equivalence principle, the heats of forma- 
tion of two compounds are found which have not been recorded in the literature, TasO, and 
ZrO. In addition, the contact potential between barium oxide cathodes and seven various 
oxides are determined and the work function of the seven oxides computed. 


PART I 
Introduction 


OR some years it has been well known that, 
when a metallic oxide, such as nickel oxide 

or copper oxide, is bombarded by electrons, this 
oxide will decompose releasing large quantities 
of gas. If such bombardment of an oxide takes 
place in a vacuum tube, this may result in 
“gassy’ tubes and ‘‘low emission’ tubes.! 
Usually, the emission in the tubes will noticeably 
decrease and the gas pressure may increase. This 
falling emission will be shown in a series of 
experiments later in this paper. Although it has 
been thought that oxides decompose upon elec- 
tron bombardment, there has not been very 
much reported as to the fundamental manner by 
which the phenomenon of the release of gases 
from oxides upon electron bombardment occurs. 
In this report, we shall propose an hypothesis, 
and then subject this hypothesis to experiment. 
An electron in motion will have a certain kinetic 
energy. When this electron bombards a surface 
film of oxide, it transfers this energy to the 
‘ bombarded material. As the kinetic energy of 
the electron reaches a critical value, the oxide 
will decompose. We propose that the critical 
energy of the electron required for decomposition 
of the oxide corresponds to the heat of formation 


* Summary of thesis submitted in partial fulfillment of 
the requirements for the degree of doctor of philosophy 
at New York University, July 10, 1945. 

** The author is now employed at Sylvania Electric Pro- 
ducts, Inc., Kew Gardens, New York. 

'L. B. Headrick and E. A. Lederer, Phys. Rev. 50, 
1094 (1936). 
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of one molecule of the oxide. We shall determine 
the energy of decomposition by bombarding the 
oxide at various applied voltages and observing 
emission from an oxide coated cathode. As the 
bombarded material emits gases, these gases 
“poison” the oxide coated cathode and falling 
emission is observed as a function of bombarding 
voltages. 

If a relationship can be shown between the 


‘heats of formation of molecules and the energy 


required for dissociation, we shall be able to 
determine the heats of formation of molecules 
whose values have not been determined by 
physical chemical methods. However, before 
discussing the hypothesis further, we shall discuss 
in some detail the nature of alkaline earth oxide 
cathodes, since this is the tool used throughout 
all of these experiments. 

A great deal of work has been done on the 
alkaline earth oxide cathodes during the nineteen- 
twenties and nineteen-thirties. It was found that 
because of an extremely low work function these 
materials produced a copious electron emission 
at relatively low temperatures.? These cathodes 
are usually composed of physical or chemical 
mixtures of barium oxide and strontium oxide 
(BaO+SrO). Sometimes, another mixture is used 
containing a third component in addition to the 
other two, calcium oxide (CaO). In the rest of 
this paper, we shall use barium in the written 
work, but it should be remembered that the 
properties of calcium and strontium can _ be 


2 A. Wehnelt, Ann. d. Physik 14 (1904). 
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described throughout in a somewhat similar 
fashion. We choose only to speak of barium 
since this is the component which has the lowest 
work function* and is most often referred to in 
the literature. 

To make an indirectly heated cathode, the 
alkaline earth carbonates are often sprayed on a 
nickel sleeve and mounted in a tube. During 
exhaust the cathodes go through a process known 
as “activation.”” The cathodes “are heated to 
about 1000°C and the carbonate decomposes 
according to the reaction 


BaCO;—BaO+ COs. (1) 


The carbon dioxide is pulled out through the 
exhaust tubing leaving a deposit of barium oxide 
on the nickel sleeve. Upon further heating and 
application of voltage to tube elements, free 
barium is formed, some of which diffuses to the 
surface of the cathode in a process which will 
be discussed later. It is the presence of this free 
barium which is responsible for the low work 
function and high electron emission from the 
cathode.‘ In the “‘activation”’ of a cathode, there 
are many processes at work, all of which combine 
to form a sort of equilibrium. 

Some of the factors at work are as follows. On 
the side of those factors producing free barium, 
we have 

(a) reduction by core metal,' 


(b) electrolysis of the cathode coating,* and 
(c) reduction of the cathode by gases.” 


Several factors have been indicated by experi- 
ments as hindering the production of free 
barium. 


(a) High voltage operation has brought about a loss of 
surface barium in the form of sputtering by ions 
formed from residual gases in the tubes. 


3 J. P. Blewett, Gen. Elec. Bull. No. 1001 (December, 
1931). 
oe A. Becker and R. W. Sears, Phys. Rev. 38, 2193-2213 
1931). 

5 E. F. Lowry, Phys. Rev. 35, 1367-1378 (1930). 

6 J. A. Becker, Phys. Rev. 34, 1323-1351 (1929). 

7J. H. DeBoer, Electron Emission and Absorption 
Phenomenon (Cambridge University Press, Cambridge, 
England, 1935), p. 365. ‘ 

8A. W. Hull, Gen. Elec. Rev. 32, 213 and 390 (1929). 
Massachusetts Institute of Technology Department of 
Electrical Engineering, Applied Electronics (J. Wiley & 
Sons, New York, 1943). 
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(b) Residual gases have been found chemically to com- 
bine with free barium. Furthermore, the barium 
oxide is quite permeable to these gases®; i.e., 


2Ba+0.—2Ba0. (2) 


(c) Some free barium is continually being lost by evapo- 
ration.!° 


Thus we have a balance of factors, the slightest 
disturbance of which can cause large variations 
in emission. This balance was the indicator used 
in detecting at what voltages dissociation of 
oxides started to occur. 


PART II 


Method of Determining Critical Kinetic Energy 
of Electrons Required for Falling Emission 
to Occur 


A. Vacuum Tubes 


Tubes were made as per sketches in Fig. 1. 
These tubes consisted of a Nonex glass envelope 
with tungsten leads in which a rotating nickel 
cylinder that could be turned by a magnet 
outside the tube was assembled in the mount 
structure. On this cylinder was welded a strip 
of the metal to be investigated. This strip was 
as long as the cylinder and somewhat wider than 
the cathode. Great pains were taken to insure 
purity and lack of contamination of the strip. 

The metal strips chosen for test purposes were 
taken from metals such as are used for spectro- 
scopic standards. All free metal parts were 
acetone washed, rinsed in distilled water, acid 
washed, distilled water washed, acetone-rinsed 
again, and hydrogen fired for ten minutes at 
800°C. During sealing, N» flushing gas and low 
fires were used to prevent oxidation. 


B. Exhaust 


The tubes were then sealed on an exhaust 
system and the metal was oxidized at exhaust. 


9 J. A. Becker, Phys. Rev. 34, 1323 (1939). Several tubes 
were made by the author containing an oxide coated 
cathode and a square anode (1 cmX1 cm) a distance of 
4 cm from the cathode surface. These tubes were placed 
on a vacuum system and emission measurements made. 
In the first tube with the pressure<one micron, the 
emission was 12ya at 5 v on the anode and 650 ya at 30 v. 
In the presence of 15 yu of oxygen, the emission dropped to 
0 wa at both 5 v and 30 v. In the second tube, the emission 
read 15 wa and 360 wa at a pressure<one micron. In the 
presence of three microns of oxygen, the emission read 
3 ya at 5v and 1.9 ya at 30 v. All measurements were 
made at 7.0 v, .66 amp. on a 6.3-v filament. At lower 
wattages on the cathode filament the effect of oxygen was 
found to increase. 

10. R. Koller, Phys. Rev. 25, 671-676 (1925). 
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Magnet to 
rotate nickel cylinder 





Getter - 


Nickel Cylinder 


Glass envelope 







Cathode 
(clone spacing to anode) 


Metal strip to be studied 


Support wire 


Fic. 1. Experimental tube for measurement of dissociation energies of surface films of oxides. 


This exhaust system contained a mercury cut-off 
and an AgO flask. A high frequency furnace 
was used to heat treat the anode for 30 seconds 
at 800°C brightness in 500 » of oxygen. 


C. Method of Measurement 


The filament in the cathode sleeve was usually 
adjusted to roughly about half the rated voltage, 
i.e., 3.5 v, .38 amp instead of 6.3 v, .65 amp. 
This would bring the temperature down to an 
estimated temperature region of 350° to 450°C. 

The metal oxide to be studied had been kept 
free of contaminations from the cathode activa- 
tion and aging by rotating it away from the 
cathode. After the cathode temperature was 
’ lowered, the sample was rotated by means of a 
magnet to a position directly in front of the 
cathode. The use of low cathode temperatures 
was essential because 

(a) metallic evaporation from the alkaline earth cathode 
could be kept to a minimum, and 


(b) the effect of oxygen in inhibiting emission was more 
marked with lower temperatures of the cathode. 


A voltage was then applied to the anode and 
increased in }-volt steps for one minute each. 
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At lower voltages, the current would rise with 
the voltage and not decay with time; however, 
at a critical voltage, the emission would slowly 
start to fall. Beyond this voltage, the decay in 
emission became more marked. 
voltage was noted. 

By varying filament currents, the cathode 
emission could be varied. It was found that 
falling emission is independent of the wattage 
but is only dependent on electron volts of energy 
applied. Furthermore, time was not important, 
for below the critical voltage, the emission could 
be maintained. 

Now the voltage applied is not the true 
measure of energy of the electron because of 
contact potential and thermal kinetic energy of 
the electrons. The true energy of electrons was 
approximated as follows. The measurement of 
falling emission was made. This was determined 
by the beginning of falling emission as recorded 
by a milliameter. The smallest amount of falling 
emission which could be determined by this 
instrument was in the order of one microampere. 
Therefore, the anode was biased down to that 
voltage (E,) on the anode such that emission 


This critical 
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from the cathode of one microampere was just 
reached. This latter voltage was arbitrarily called 
the retarding cut-off potential (E,). Then E,—E, 
was used as the real potential (Fig. 2). 
E,—E.=E real. (3) 

These measurements were repeated several 
times at different wattages and on _ several 
samples of similar types. In addition to this, a 
further check was made on sample to insure 
that the same surface was being bombarded. 
This was done by measuring the contact potential 
between each surface and the cathode. If the 
contact poteritials were the same from tube to 
tube of the same sample, the consistency of 
producing the same oxide surface, could be 
further assured. 

The method of measuring the contact potential 
was done as follows. The filament voltage was 
still further lowered and both accelerating and 
retarding potentials were applied, emissions 
measured and extrapolated. Against this, ac- 
celerating potentials were applied and extrapo- 
Jated (Fig. 3). 

The point at which these extrapolations meet 
was assumed to be the contact potential, for if 
there were no contact potential, the retarding 
potential would change to an accelerating po- 
tential at 0 volt (Fig. 4). 

Thus the shift in the theoretical curve is 
caused by a contact potential E contact between 
the work function of the barium cathode ¢ and 
the work function of the oxide being studied ¢’, 
or in other terms, the difference in the work 
functions of the two materials 








¢—¢' =E contact. (4) 
PART III 
Data 
The data in Table I were taken on tubes, 
ad 
Log Ig Pa 
| 
1 ua | 
Eo E, Potential 


Fic. 2. Measurement of electron energy. J,—emission, 
E.—applied potential, E.—cut-off potential. 
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Fic. 3. Shift in emission due to contact potential. 
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Fic. 4. Ideal emission curve with no contact potential. 


anodes of which were prepared as in the discus- 
sion above. It might be pointed out that the 
values “‘critical potential’’ for falling emission to 
occur were independent of time, distance, and 
current density of the bombarding current pro- 
vided that the anodes were not caused to increase 
in temperature to any great extent. 

Figure 5 illustrates a typical case of how falling 
emission occurs as the voltage is raised. In this 
case the voltage was raised in 0.5-volt steps and 
held at each level for three minutes. The partic- 
ular sample being studied under bombardment 
was FeO. 

In the study of falling emission while bom- 
barding oxides with electrons at various voltages, _ 
the idea came up that if the oxides were re- 
sponsible for the falling emission, a clean metal 
surface would not cause falling emission when 
bombarded. 

In order to test this hypothesis, a tube was 
made with a copper strip on the anode. Instead 
of oxidizing this strip on exhaust, however, 
hydrogen was introduced and all parts of the 
anode were heated at 750°C for three minutes 
in 500 microns pressure of hydrogen. The gas 
was removed and the process repeated. After 
this treatment, the copper appeared very bright 
and it was assumed oxides were reduced. 

Here it was seen that with pure Cu (hydrogen 
fired at exhaust) there is no falling emission. 
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TABLE I. 
ns Po- 
Potential tential 
applied - for 
tor Retarding falling 

lube Metallic Filament falling cut-off emis- 
No. oxide voltage emission - a — 
101 FeO 3.5 7.0 Vv 1.5v 5.5v 
101 FeO 3.1 7.0 1.7 5.3 
101 FeO 4.1 6.5 1.3 5.2 
103 FeO 4.0 7.0 ia 5.8 
103 FeO 3.7 7.0 1.4 5.6 
110 FeO 3.6 ran 1.95 5.55 
110 FeO 4.0 8.0 1.8 6.2 
104 CuO 34 4.5 aun 2.2 
104 CuO 4.0 4.5 2.4 2.3 
105 CuO 3.5 3.5 1.9 6 
105 CuO 4.1 3.5 2 1.7 
106 CuO 3.3 4.0 1.9 re 
106 CuO 3.5 4.0 1.9 2: 
114 MoO; 4.1 9.0 1.5 7.5 
114 MoO; 3.4 9.0 2.4 6.6 
115 MoO; 4.0 9.5 2.2 ( 
115 MoO; 4.0 9.5 yr 7.3 
115 MoO; 3.3 9.5 2.4 7.1 
116 NiO 3.6 7.0 4.2 2.8 
116 NiO 4.2 7.0 4.0 3.0 
117 NiO 3.4 4.5 1.8 2.7 
117 NiO 4.0 4.5 a 2.8 
117 NiO . & 5.0 2.2 2.8 
118 NiO p eS 6.0 3.8 2.5 
118 NiO 4.1 ao 3.3 2.2 
216 Wo: 3.3 7.0 2.1 4.9 
216 wo, 3.7 6.5 1.8 4.7 
218 wo: 4.1 8.0 3.1 4.9 
218 wo, 3.2 7.5 3.5 4.0 
218 wo. 4.5 7.5 2.9 4.6 
203 Ta20; 3.1 7.5 3.0 4.5 
203 Ta2O; 3.5 7.0 2.7 4.3 
204 Ta20; 3.4 7.0 1.8 5.2 
204 Ta2O; 3.1 7.5 2.5 5.0 
214 Ta2O; 3.3 7.5 2.8 4.7 
210 ZrO2 3.5 4.5 2.6 1.7 
210 ZrO2 4.1 a5 y Be 1.3 
211 ZrO2 4.1 3.5 1.2 1.7 
211 ZrO2 3.4 4.0 2.1 1.9 
212 ZrOz a6 . 4.0 aun 1.8 
212 ZrO2 3.3 5.0 29 2.0 
212 ZrO2 3.3 3.5 1.4 an 
212 ZrO2 3.5 a 2.0 1.5 
600 


If anything, there seems to be a slight tendency 
towards rising emission. This may be caused by 
hydrogen still remaining in the copper metal or 
it may be caused by some other effect. In any 
case, this shows that falling emission did not, 
occur when the oxide was not present. 


Summary of Data 


In Table II we have listed the experimental 
value of the energy required for falling emission 
to occur. Next to this value is the value of the 
heat of formation (#7) required for the probable 
chemical reaction to occur. It is assumed that 
23 kilogram calories per mole is equivalent to one 
electron volt per molecule. FE real=electron 
energy at first sign of falling emission and 
E contact=contact potential between surface 
studied and barium oxide cathode (Fig. 6). 


PART IV 
Conclusions and Discussion 


From the section containing the data we might 
first conclude that the phenomenon of falling 
emission in this study is related to the anode to 
which the current is being drawn. If falling 
emission were caused by something occurring in 
the cathode itself, i.e., electrolysis or caused by 
residual gases in the tubes, there would be no 
reason for the change in E, (energy necessary to 
show falling emission) from one anode material 
to another. However, this critical energy of the 
electron is different for every chemically different 
surface. Furthermore, in studying different tubes 
with the same chemically prepared anode surface, 
we find the same E,. It was also shown that when 
a “clean’’ metal was bombarded, no falling 
emission occurred. Falling emission did occur, 


TABLE II. 


Handbook H 








(Heat of Probably 
Com- E real formation E contact dissociated 
pound volts for reaction) volts to: 
2FeO 5.64.2 A.D. 5.5 volts 2.74.3 A.D. 2Fe+O.- 
CuO 2.0+.2 1.5 3.2+.1 Cu+ 40, 
MoO; 7.15+.03 7.4 3.14.2 Mo+30: 
NiO 2.74.2 °‘ 2.5 4.4+.2 Ni+40, 
WO, 4.74.15 5.5 3.8+.6 W+0:, 
TasO; 4.8+.3 unknown 3.5+.2 Ta204+ 402 
ZrO. 1.94.3 unknown 3.6+.4 ZrO+ 40,2 








A.D. =average deviation. 
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Fic. 5. Falling emission as voltage is increased 
(three-minute intervals). 


however, when this same type of metal (Cu) 
was oxidized and then bombarded at or above a 
critical voltage. We might then conclude that 
the falling emission was a function of some 
property of the oxide on the anode. 

It was seen that a certain amount of energy 
(electron volts) was needed to cause the effect. 
When the quantity of energy reaches a critical 
value, something occurs in the cathode which 
causes falling emission. We might assume that 
the bombarded molecule decomposes and that 
the free oxygen returns to the cathode, combines 
with the barium on the cathode, and lowers the 
emission. The statement that oxygen from the 
oxide is the agent of falling emission is further 
verified by a relationship in the cases of five 
metals between the heats of formation and the 
energy of the electrons at falling emission. These 
values are found to be equivalent. 

There were in general two considerations 
involved. First, an oxide was chosen which was 
very common and: very likely to be present. 
Second, the compound in the same series of 
oxides which had the next lower heat of forma- 
tion was considered. The difference in energy of 
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Fic. 6. Relationship between heat of formation of 
molecules and energy of electron required for falling 
emission. 


the two molecules was taken as the energy 
required to release oxygen in terms of the heats 
of formation. 

For iron the likelihood of the presence of FeO 
was high since this compound is common and 
easy to form. A molecule of 2FeO was assumed 
to exist. In terms of Heats of Formation approxi- 
mately 5.5 electron volts would be required to 
strip this of its oxygen. The value obtained by 
direct electron bombardment was found to be 
5.6+.2 A.D. electron volts. 

In the case of copper, it was assumed that 
the oxide was CuO. This was done since CuO is 
one common oxide and Cu,0 is the other. Now 
the H for CuO is slightly higher than the H for 
CuO (roughly 39 kg cal. as against 34). There- 
fore, if energy were supplied to a _ system, 
CusO—Cu+CuO and this reaction would not 
liberate oxygen and thus could not be perceived 
by measuring emission. However, as more energy 
was added in the form of higher voltages, the 
CuO would decompose to Cu+4Q, and it is this 
reaction which would affect the emission from 
the cathode. 

In the case of molybdenum, MoO; was com- 
mon and its heat of formation was so close to 
the experimental value that this compound was 
chosen. Why the compound seemed to decompose 
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TaBLeE III. Work functions of several oxides. 








Compound ¢ of oxide 

FeO 3.85 electron volts 

CuO 4.35 

MoO; 4.2 

NiO 5.5 
4.9 
4.6 
4.7 





wo, 
Ta2Os . 


> 
5 
5 
5 
ZrO, 5 


completely to molybdenum and oxygen instead 
of going down the transition steps, i.e., to MoO, 
is not known. 

NiO was chosen for nickel since it is the most 
prominent form and because of good agreement 
with the H value. 

WO, was chosen for tungsten since it is a 
common form of oxide and since the agreement 
was fairly good with the H listed in the Hand- 
book. It was also noticed in the making of the 
tube that the yellow WO; had not been reached. 
In one tube not recorded in the data, the oxi- 
dizing process was allowed to continue and it 
took from two to four times the time used to 
obtain the yellow color. 

In considering tantalum and zirconium the 
same line of reasoning was used. It was assumed 
that during the oxidation process a common 
oxide was formed. Since zirconium and tantalum 
both have a great affinity for oxygen, it was 
assumed higher oxides were formed, i.e., ZrOs 

and Ta.Osx. 

' The H (heat of formation) for both of these 
compounds is known: the H for their lower 
oxides, however, have not been reported. The 
experimental values for the higher forms did not 
coincide with their values of H; therefore, just 
as in the previous cases, it was assumed a low 
oxide form was reached in the reducing process, 
rather than free metal. 


For tantalum and zirconium the next common 


lower oxides were considered ZrO and Ta,Q,. 
The energy measured thus would be for the 
.following reactions 


Ta,O;-Ta,O.+ $02, 
ZrO.—Zr0+ 402. 


Using our experimental values of H, we arrived 
at the conclusion that the lower oxides had the 
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following values ot H: 


Ta.O; H—300 kg cal. by Handbook 
Ta,04+ 402 H—110 kg cal. by measurement 
TaO, H—190 kg cal. by calculation 
ZrOz H—178 kg cal. by Handbook 
ZrO+ 30, II— 43 kg cal. by measurement 
ZrO II—135 kg cal. by calculation 


PART V 


Further Considerations of Conclusions 


(1) Conclusion No. 1 concerns the mechanism 
by which falling emission occurs. It describes a 
manner in which oxides may decompose when 
bombarded in vacuum and explains why oxidized 
parts are in general harmful when used in 
vacuum tube construction when subjected to 
electron bombardment. It also suggests the 
possibility that other compounds may decompose 
in a similar manner upon bombardment, i.e., 
salts, silicates, etc. As far as the author knows, 
no tests of this type have been reported in the 
literature. 

(2) Utilizing the data concerning contact po- 
tentials the work functions of various oxide films 
may be determined. 

By ¢ oxide=E contact+1.15 volts (¢’ of 
Ba: BaQO)." 

Values found experimentally for the contact 
potentials were used to compile Table III. 

(3) We have found in the methods outlined a 
means of predicting the heats of formation of 
two compounds the H value of which has not 
been reported in the literature. 


H Ta,0,=190 kg cal. per mole, ~ 
H ZrO =135 kg cal. per mole. 


This method may possibly be used further in 
determining the heats of formation of other 
compounds whose values are undetermined. 

There is one more point worthy of mention. 
While the voltages found necessary to decompose 
the oxides correspond to the heats of formation, 
it may not be the actual heats of formation that 
are being studied, but rather the free energy of 
dissociation. In this case however, the heats of 

"J. H. DeBoer, reference 7, p. 357. This value of work 
function is not identical with all experiments. Thus there 


is an error probably of +.3 volt in this value and the 
values of oxides listed below. 
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formation and the free energies of dissociation 
are so close that the means by which we measure 
these quantities would not separate them. 
Actually then, we may be reading free energies 
of dissociation by our measurements rather than 
heats of formation. 
Lewis and Randall give the following relation- 
ship” 
F=H-—TS=(E+PV)-TS, 
where 
F—free energy per mole, 
H—heat of formation per mole, 
T—temperature, 
E— internal energy, 
P—pressure, 
V—volume, and 
S—entropy. 


It is known that there is a slight variation in 
F and H with P and T, but in our case P=0 by 


2G, N. Lewis and M. Randall, Thermodynamics 
(McGraw-Hill Book Company, Inc., New York, 1923). 








approximation, S is small for metallic oxides," 


- even though the temperature of the anode oxides 


is slightly higher than room temperature. Hence, 

F is nearly equal to H which is nearly equal to 

E, and we have a rough means of measuring 

these quantities in the method described above. 
(a) Free energy of dissociation per mole, F; 


(b) Heat of formation per mole, H; 
(c) Internal energy per mole, E. 


Acknowledgment should be made to Dr. E. 
A. Lederer for suggesting this problem and to 
Mr. D. H. Wamsley, Dr. L. B. Headrick, Mr. 
G. Widell, and Dr. L. Malter of the RCA 
Laboratories and Professors S. A. Korff, C. W. 
Barnes, and A. F. Myers of New York University 
for their suggestions pertaining to this work. 
~ 18 Bulletin 384, U. S. Department of the Interior, 
Bureau of Mines (Entropy of CaO is reported to be 9.5 
calories at 298.1° absolute, page 16. Entropy for CuO— 


10.4 calories at same temperature, page 19. Entropy for 
FeO—14.2 calories at same temperature, page 20). 





An Experimental Investigation of Forced Vibrations in a Mechanical System 
Having a Non-Linear Restoring Force 


; CarL_ A. LUDEKE 
Department of Mathematics and Mechanics, University of Cincinnati, Cincinnati, Ohio 
(Received April 4, 1946) 


This paper introduces a mechanical apparatus capable of generating and recording forced 
vibrations in a system having a non-linear restoring force. The experimental wave forms are 
then compared with the theoretical results given by three graphical methods. These are the 
methods of Martienssen, Hartog, and Rauscher. Provided a suitable first guess for the ampli- 
tude can be made, one application of Rauscher’s method gives better results than the other 
methods. A graph is made of the experimental amplitude as a function of the disturbing 
frequency. These results are compared with the theoretical results of the above three methods. 
Again one application of Rauscher’s method gives satisfactory results. In all experimental 
results it is noted that even though the restoring force is distinctly non-linear, the wave forms 
of the resulting motion are nearly sinusoidal as long as the frequency of the observed motion is 
the same as the frequency of the disturbing force. However, steady oscillations can be main- 
tained for which the observed frequency is a sub-multiple of the disturbing frequency. Two 
such subharmonics were recorded and the experimental wave forms are shown. 





ONSIDER an undamped system of one 

degree of freedom with a curved spring 
characteristic and under the influence of a har- 
monic disturbing force. The equation of motion 
is given by: 


m(d?x /dt?) +f(x) = Po cos wt, (1) 
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in which m is the mass of the system, f(x) 
represents the spring characteristic, Po is the 
amplitude of the disturbing force, and w is the 
angular frequency of the disturbing force. An 
exact solution of Eq. (1) exists only for some 
very special forms of f(x). There are however 
three approximate methods available, namely: 
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Fic. 2. Diagram of experimental arrangement. 


the methods of Martienssen,' Hartog,? and 
Rauscher.* In each of these methods f(x) need 
only be known graphically, and the methods 
themselves can then be completed by an entirely 
graphical procedure. The method of Martienssen 
gives a one-term solution, the method of Hartog 
gives a two-term solution, and the method of 
Rauscher gives a graphical solution capable of 
being repeated until the desired accuracy is 
obtained. In this paper only one application of 
the Rauscher method is used. It would be of 
interest therefore to compare the results of the 
three approximate methods with exact solutions 
if they were in general available. Since this is 
not the case, they are compared with experi- 
mental results. 

These results were obtained by means of the 
system shown in the diagram of Fig. 1. By 
changing springs a variety of load deflection 
curves could be obtained; all however of the 
“increasing stiffness’ type. Every effort was 
made to keep the damping forces to a minimum. 
The unbalance was rotated with a constant 
_angular velocity by means of an electric motor. 


1Q. Martienssen, Physik. Zeits. 11, 448-460 (1910). 

2 J. P. Den Hartog, “The amplitudes of non-harmonic 
vibrations,”’ J. Frank. Inst. 216, No. 4, 459-473 (October, 
1933). 

3M. Rauscher, “Steady oscillations of systems with 
non-linear and unsymmetrical elasticity,” J. App. Mech. 
5, A169-A177 (1938). 
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Fic. 3. Load deflection curve. 


The equation of motion for this system is then 
1 (d°0/dt*) +f(0) = MrLe’* cos wt=Pycoswt, (2) 
in which 


I is the moment of inertia of the vibrating parts =.31 
slug-ft.?, 

(0) is the spring torque, 

w is the constant angular velocity of the unbalance, 

M is the mass of the unbalance, 

r is the radius of the unbalance, 

L is the distance from the center of rotation of the 
unbalance to the center of rotation of the beam, 
and 

Po = MrLo*. 


A diagram of the experimental set-up is shown 
in Fig. 2. A typical load deflection curve is 


oe © @ «@ > - © © ¢ © © © © & © «2 
0 I 
owed 
O1 rad. 
x J 
7 . . ~~ * . . . . . - « 
Fic. 4. Recorded wave form when using the load 


deflection curve of Fig. 3. w=64.0 cycles/sec.; Po=.872 
ft.-lb. 
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shown in Fig. 3. This spring characteristic was 
used when recording the wave form shown in 
Fig. 4. The series of dots at the bottom of the 
picture represent time intervals of 1/30 sec. The 
series of dots at the top represent the position of 
the unbalance at 90° intervals. The heavier dot 
indicates that the unbalance is vertically upward, 
i.e., that wt=0. The simplest approximate solu- 
tion! of Eq. (1) is based on the assumption that 
the motion x=f(t) is sinusoidal and has the 
“forced” frequency. That is, the solution is 
assumed to be 

xX =X 9 COS wt. 


(3) 
Using this solution, the forces involved are the 
inertia force —mw’xy) cos wt, the external force 
Py cos wt, and the restoring force f(x). We note 
that these are all zero when x=0, i.e., Eq. (1) 
is satisfied by Eq. (3) when x=0. We shall now 
require the solution to satisfy Eq. (1) when 
x=X». Thus 


—muxot+f(xo) =Po 
f (x0) = morxot+Po. (4) 


A very convenient way of finding the value of 
xo which satisfies Eq. (4) is as follows: Plot f(x) 
load versus deflection, superimpose on this curve 
a straight line having a load intercept P» and a 
slope mw’. The point of intersection of the 
straight line and the curve will determine Xp. 
This value is then used in Eq. (3) to give the 
approximate solution. Obviously this method is 
only true for all values of ¢ when f(x) is linear, 
for any other form of f(x) it is true only at the 
middle and the end of the “stroke.” 

Figure 5 shows a load deflection curve with 
the construction lines for the one-term solutions. 
The experimentally recorded wave forms are 
shown in Fig. 6, together with the Martienssen 
solutions. The Hartog solution uses the same 
basic idea of the Martienssen solution, namely, 
to hold the frequency fixed and to choose an 
amplitude corresponding to it; but it extends 
the solution to two terms instead of one term. 
This allows us ‘to balance the forces at three 
points of the “stroke’’ instead of at only two 
points of the stroke. 

In this method, the solution of Eq. (1) is 
assumed to be of the form 


or 


x=, COs wt +a; cos 3, 


(5) 
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and f(x) is assumed to be of the form 
f(x) =A, cos wt+A3 cos 3wt. (6) 
Substituting Eqs. (5) and (6) in Eq. (1) we have 


(A, —mw*a;— Po) cos wt 

+(A3;—9mw*a3) cos 3wf=0. (7) 
Equation (7) must be satisfied for all values of ¢. 
Thus 


A;=mo'a,;+Po, As=9mw"a3. (8) 
Therefore Eq. (6) takes the form 
f(x) = (mw’a,+Po) cos wt +9mw*a; cos 3wt. (9) 


We must now choose the a’s in Eq. (5) so that 
when the A’s are calculated by means of Eq. (8), 
f(x) as represented by Eq. (9) turns out correctly. 

We notice that when wt=2/2, Eq. (7) is 
satisfied whether or not Eqs. (8) are satisfied. 
Since Eqs. (8) contain two parameters a; and ds, 
we can satisfy Eq. (7) at two other values of wt. 
Let us choose wf =0, and wt =7/6. 








At wt=0, x=ai1+4a;3; 
f(x) = mw*a,+ Pot+9mw'a3 
. } (10 
at w=x/6, x=Vv3a;/2; 
f(x) = (mw*a,+Po)v3/2. 
a 
o:658 
R =1704 
24, 
@e- 64 
R = .872 
20- Y 
3. : rf 
2 SP 4 = 57.8 
3 // ® re =e 
£ oe , 
zr ee 
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J al on © @:458 
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Fic. 5. Load deflection curve with construction lines for 
one-term solution. 
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Fic. 6. Wave form recorded when using the load deflection curve of Fig. 5. 


These relations can be written in the form 


a,v3 a,;v3 v3 
ni ) =ma*( )+Pr 
2 2 2 


f (ai +43) = 9mw*?(a,+a3) +(Po—8mw'a)). 








(11) 


The right sides of Eqs. (11) represent straight 
lines with slopes mw* and 9mw’, respectively, and 
having intercepts Pov3/2 and (Po—8mw*a;), 
respectively. The left sides of Eqs. (11) represent 
the function f(x) evaluated at x=a,v3/2 and 
x=a,+a;. Thus a, and as can be solved for 
graphically as shown in Fig. 7. The results for 
the load deflection curve of Fig. 5 are shown on 
Fig. 6. 

In both methods presented thus far the pro- 
cedure was to determine an amplitude corre- 
sponding to a given frequency. Perhaps the most 
exact method? available does just jthe opposite. 
It chooses an amplitude and then ‘solves for the 
corresponding frequency. It is a method which 
can be repeated until the desired accuracy is 
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Fic. 8. Load deflection curve. 











Fic. 9. Load deflection curve. 


In this method Eq. (1) is reduced to an equa- 
tion in which the right side is zero, i.e., 


m(d?x /dt*?) + R(x, t) =0, (12) 


VOLUME 17, JULY, 1946 











Fic. 10. Load deflection curve. 


in which R(x, t)=f(x)—Po» cos wt. In a system 
in which f(x) is linear, the solution is x =A cos wtf, 
hence R can be written as a function of x alone, 
i.€., 


R(x) =f(x)—Pox/A. (13) 


For systems which are slightly non-linear this 
equation can be assumed valid. Thus Eq. (1) 
takes the form 


m(d?x /dt?) + R(x) =0, 


in which R is given by Eq. (13). 
This equation of a free vibration is solved as 
follows: A first integral of Eq. (14) gives us: 


(14) 


2 4 dx 
oo— f R(x)dx; v=—. (15) 
md, dt 


The time corresponding to any displacement x 
is then found from the equation 


A dt A dx 
a f = 
« a 2: iv 


(16) 


For a purely graphical evaluation of this integral 
a change of variable is introduced by letting 
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Fic. 11. Amplitude versus frequency when using load 
deflection curve of Fig. 8. 


x=A sin ¢. Thus Eq. (16) becomes 


"/2 A cos¢ 
i= f ————d¢. (17) 
a v 


The procedure is therefore as follows: After 
choosing A, a table of values of x and v is made 
by graphically integrating Eq. (15). The variable 
x is changed to ¢ by means of the relation 
x=A sin ¢, and a new table of values of ¢ and ¢ 
is made by graphically integrating Eq. (17). 
From these two tables a plot is made of x and f¢, 
and the frequency corresponding to the choice of 
amplitude A is then determined. The entire 
process can then be repeated by using the 
relationship between x and ¢ to determine a more 
exact value of the function R(x) used in Eq. (15). 

Having the experimental results as shown in 
Fig. (6), the initial choice of A is greatly simpli- 
fied. Choosing A equal to the experimental values 
-of the amplitude, the results of one application 
of this method are shown on Fig. 6. 

In Figs. 8-10 are shown three load deflection 
curves of varying non-linearity. On these curves 
are superimposed the construction lines for one 
term solutions. From these results we can plot a 
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deflection curve of Fig. 9. 
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Fic. 13. Amplitude versus frequency when using load 
deflection curve of Fig. 10. 


relationship between amplitude and frequency. 
The experimental results together with the theo- 
retical results of all three methods are shown in 
Figs. 11-13. The lack of agreement between the 
experimental results and theoretical results for 
the small amplitudes observed at frequencies 
higher than the “jump” frequency is caused 
primarily by the errors involved in recording 
the small amplitudes and in determining accu- 
rately enough the load deflection curves for very 
small amplitudes. 

In the three theoretical methods illustrated, 
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14. Sub-harmonic resonance wave forms recorded 
when using the load deflection curve of Fig. 15. 


Fic. 


the initial assumption was made that the solution 
of Eq. (1) should have a frequency approximately 
equal to the frequency of the disturbing force. 
This, however, does not cover all possible solu- 
tions. In fact one of the most interesting phe- 
nomena observed in non-linear vibration systems 
is that of sub-harmonic resonance.‘ This is the 
phenomena of a steady forced vibration whose 
frequency is a sub-multiple of the disturbing 
frequency. The recorded wave forms of sub- 
harmonic resonance (3) and (2) are shown in 
Fig. 14. These were observed when using the 
load deflection curve shown in Fig. 15. The three 
theoretical methods as used in this paper would 
not apply to these wave forms. 

Before drawing final conclusions, let us empha- 


*C. A. Ludeke, ‘“Resonance,”’ J. App. Phys. 13, No. 7, 
418-423 (July 1942). 
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Fic. 15. Load deflection curve. 


size the fact that the theoretical methods were 
applied as if the system contained no damping. 
This is of course impossible to achieve experi- 
mentally. Also, there is considerable experi- 
mental error when observing frequencies higher 
than the jump frequency. Bearing this in mind, 
however, we can draw the following conclusions. 
First, even though the load deflection curve 
varies considerably from a straight line, as 
indicated in Fig. 3, and even though the ampli- 
tude of the maintained vibration is quite large, 
as indicated in Fig. 6, A=0.0171, nevertheless 
the resulting wave forms, Fig. 6, are very nearly 
sinusoidal. This means that the very simple 
approximation of Martienssen is nevertheless a 
very good one, at least for the cases herein 
considered. Second, though the two-term method 
at some particular frequency gives the exact 
amplitude, see Figs. 11-13, nevertheless the 
values depend so critically on the values of Po 
and w.that for the most part the results are not 
as good as the one-term solution, see Fig. 6. 
Third, provided the experimental amplitude is 
known, the short cut method of Rauscher can 
be applied once with excellent results. Fourth, 
as soon as the case of subharmonic resonance is 
considered the recorded wave forms are no longer 
sinusoidal. 

It is a pleasure to thank Miss M. Seuberling 
for her assistance in preparing the drawings. 
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A New Method for Measuring Dielectric Constant and Loss in the 
Range of Centimeter Waves* 


S. ROBERTS AND A. VON HIPPEL 
Laboratory for Insulation Research, Massachusetts Institute of Technology, Cambridge, Massachusetts 


(Received April 29, 1946) 


In 1940, dielectric measurements in the centimeter range were considered as difficult and not 
very accurate. The authors, therefore, developed a “‘hollow-pipe” method which overcame these 
objections and required only a weak oscillator and small amounts of the dielectric material. 
The theory and its practical applications, as perfected by March 1941, are presented in this 


paper. 


HE capacitance of a condenser is propor- 

tional to the dielectric constant or per- 
mittivity «* of the dielectric material which it 
contains. 


C= (e*Co) /€o, (1) 


where Cp is called the vacuum capacitance and 
€o is the permittivity of free space. When sub- 
jected to a sinusoidal voltage V, a condenser 
passes a current J which normally does not lead 
the voltage by 90° but only by (90°—4) (Fig. 1). 
If the current is defined in the general way 

* 


€ 
I = jwCV = jwCo—V, (2) 


€0 


the existence of the loss angle 6 can be repre- 


sented by introducing a complex permittivity 
e* =e’ —je’’. (3) 


The charging current is proportional to e’, the 
loss current to e’’, and the loss tangent 


tan 6=e''/¢’ (4) 


is given by the ratio of loss current to charging 
current. e’ and e”’ or e’ and tan 6 are the param- 
eters normally used for describing the frequency 
response of a dielectric material. The loss tangent 
is identical with the power factor cos @ for small 
angles 4. 





* This research, about which a preliminary report was 
given at the meeting of the American Physical Society in 
April, 1940, was ready for publication by Massh, 1941, but 
had to be classified and withheld from general circulation. 
The method, now one of common knowledge, has been 
refined further in our laboratory during the war and is 
being used with coaxial lines and hollow wave guides in 
the range from one meter to one centimeter. A survey of 
this newer development and of the extensive dielectric 
information obtained will be given elsewhere. 
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The “lumped” circuit aspect used above be- 
comes invalid in the ultra-high frequency range, 
but the concept of a complex dielectric constant 
is not affected. It describes the physical fact 
that a dielectric in a periodic electromagnetic 
field carries an electric current of the density 


J =cE+aD/at=(e' —je’)joE, (5) 


consisting of two components, a charging current 
of the density 


J.=¢'jwE (6a) 


storing electric energy by the reversible displace- 
ment of charge carriers, and a loss component 


Ji = éwE (6b) 


measuring the dissipation of energy in heat 
caused by traveling charge carriers as well as 
by the time lag between the electric field E and 
the dielectric flux density D. 

The dielectric constant and loss can be found 
indirectly in the microwave range by transmis- 


_sion line or optical methods. In the former e* is 


derived from the intrinsic wave impedance Zo or 
propagation function 7» of the dielectric medium, 
while in optical methods it is evaluated from the 
complex index of refraction n* equal to n(1—jx). 
Standing or traveling waves can serve for meas- 


‘ ry — I 
J—=WC,V 
I ia 
! 
* $ 1 
v om Co 
] ) : 
Vv 
$ Woo 
Fic. 1. Relation between current and voltage in a condenser 
with loss. 
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uring any one of these interrelated sets of 
parameters. 

The normal optical-procedure employs travel- 
ing waves. The intensities of the reflected and 
transmitted beam are measured and m and «x 
evaluated by Fresnel’s equations. Esau, with his 
co-workers Baz! and: Kebbel,? has carried this 
method over into the ultra-high frequency 
range, but here it becomes definitely inferior. 
Serious boundary problems arise because the 
wave-length is of the dimensions of the sample, 
and standing waves are set up between trans- 
mitter and receiver, producing intensity fluctua- 
tions by interference. 

The macroscopical wave-length, on the other 
hand, provesof great advantage if standing 
waves, that is interference optics, are applied 
because a detector may travel directly through 
the profile of the.wave pattern. Drude’s two 
classical methods* utilize this possibility and 
have since been employed in many variations.‘ 
However, a handicap of the Lecher system re- 
mained: an empirical calibration of the con- 
denser system was necessary, and extreme care 
had to be taken to avoid perturbation of the 
waves by the detector system and to shield the 
detector against stray fields. These limitations 
have been overcome by the development of a 
“‘hollow-pipe’’ method, to be described in the 
following paragraphs. 


PRINCIPLE OF THE METHOD 


By limiting the electromagnetic field to the 
enclosure of a hollow pipe or coaxial line, all 
boundary and stray effects disappear auto- 
matically and small amounts of any dielectric 
can be measured with precision. With this con- 
cept in mind we developed the apparatus shown 
schematically in Fig. 2. A transmitter radiates 
waves of a given frequency into one end of a 
closed wave guide; they are reflected by the 
metallic boundary at the other end. Standing 
waves are set up and can be measured by a probe 
detector traveling along a slot in the pipe parallel 





1G. Baz, Physik. Zeits. 40, 394 (1939). 

2 W. Kebbel, Hoch:tech. u. Elek:akus. 53, 81 (1939). 

3P. Drude, Wied. Ann. 55, 633 (1895); 61, 466 (1897). 

‘See, for instance, G: Potapenko, Zeits. f. Physik 20, 
26 (1923); M. Seeberger, Ann. d. Physik 16, 77 (1933); 
R. King, Rev. Sci. Inst. 8, 201 (1937); K. Slevogt, Ann. 
d. Physik 36, 141 (1939); H. Slatis, Ann. d. Physik 36, 
397 (1939). 
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to its axis. The dielectric is inserted in the closed 
end of the pipe opposite to the transmitter, filling 
the volume to a height d. Above it the wave 
pattern is measured in air. e’ and e”’ are calculated 
from the ratio of the field strength in node and 
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Fic. 2. Diagram of the apparatus showing principle of 
the method. 


antinode Emin/Emax and the distance xo of the 
first node from the surface of the dielectric. 

It is a known procedure to calculate the ter- 
minating impedance of a transmission line from 
the voltage ratio Vimin/ Vmax and the distance xo, 
or an acoustic impedance, by the location and 
relative magnitude of the sound-pressure maxima 
and minima.® But measurements of the terminal 
impedance of a short wave guide filled with a 
dielectric material have not been applied, so 
far as we know, to the evaluation of e’ and e’’. 
We are giving, therefore, in the following pages 
a mathematical theory of the method, the de- 
scription of our apparatus, and some results 
illustrating its performance. 
$H. O. Rosenstein, Zeits. f. Hochfrequenztechnik 36, 
81 (1930); N. N. Malov, Physik. Zeits. Sowjetunion 11, 


539 (1937). 
*W. M. Hall, J. Acous. Soc. Am. 11, 141 (1939). 
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MATHEMATICAL THEORY 


A standing wave in the pipe can be repre- 
sented as the sum of two traveling waves pro- 
gressing in opposite directions parallel to the 
axis of symmetry. The transverse field com- 
ponents E(x), and H(x), in air (medium 1) at a 
variable distance x from the dielectric surface 
along a line parallel to the axis are given below. 


E(x), =Aa exp (mx) +An exp (— yx) 
=Aalexp (x) +r exp (— yx) ], 
Aa A 


rl 
H (x), =—— exp (yix) —-—— exp (— y1%) 
Zi Zi 


(7) 


Aa . 
=—T([exp (yx) —7o exp (— yx) ]. 
Z\ 


Aj, and A,,; are the amplitudes of the incident 
and reflected waves at 
Their ratio 


the dielectric surface. 


An/A a=nm=e**, (8) 
where 


o=p+ jy, 


defines a reflection coefficient which character- 
izes the terminating pipe section containing the 
dielectric medium 2. The wave impedance Z(0) 
of this section, given by the ratio of the electric 
to the magnetic field intensities at the dielectric 
surface, is related to ro and 9. 


E(0) 1+ 7 1+e-** 
Z(0)= “ =Z, = Z) : -=Z, coth d. (9) 
H(0) 1—7r l—e* 


Normally the attenuation in the air-filled part 
of the pipe can be neglected, that is, the propaga- 
tion 1, is reduced to its phase factor 


¥1=jBi=j2n/d, (10) 


where A; is the wave-length in air inside the 
wave guide. Since y: is imaginary, that is, 
_lexp(+yix)| =1, it follows from (7) that the 
maximum and minimum field amplitudes of the 
standing wave are 
Emax = |Ai|(1+ |r0|)=|Ai|(1+e-*), 
(11) 
Emin = | Ai| (1— |r|) =|A;i|(1—e7-*), 
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and their ratio is 


Beste 1 —e* 
——o =tanh p. 
Eusx 1+6°% 








(12) 


The first minimum above the dielectric bound- 
ary occurs at a distance x» where the incident 
and reflected electric waves have a phase dif- 
ference of w radians, that is, 


2 WX 2 TXo 
- = —+2y—— —T, (13a) 
Ay ° M1 
or 
Y=2r(} — Xo, ys). (13b) 


The wave impedance Z(0) can now be expressed 
in terms of the directly measurable parameters 
xo, Ar, and Enin/Emax. Expanding coth ¢ in Eq. 
(9) one obtains 


tanh p—j cot y 
Z(0) = 2,——____—_—_- 
1—j tanh p cot ¥ 


Ente 
—j tan 
Enex Ai 
=Z; . (14) 
Ewin 24X0 
i-j tan —— 
Smex A 


2 TXo 














The impedance Z(0) is determined by the 
length d and the propagation function y2 of the 
short section of pipe containing the dielectric, 
or conversely, if Z(0) and d are known, 72 can 
be derived from them. 

The standing wave in the dielectric is 


E(x)2=A 2 exp (vox) +A 2 exp (— 2%), 


A i2 A r2 
IT(x)2=—— exp (y2x) —-—— exp (— 2%), 
Z Z2 


9 


(15) 


and at x = —d (Fig. 2) the pipe is terminated by 
a metal plate, introducing at this boundary a 
voltage node. 


0=A,» exp (—yed) +Anexp (yd). (16) 


Solving (16) for A,2/Aj2 and substituting this 
value in (15) gives 


E(0)s 
Z(0) =——— = Z, tanh yd. 
H( 


2 


(17) 
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Fic. 3. Chart of the function Ce =tanh Tei7/Te’’, with ¢ and r expressed in degrees. 


The characteristic impedance Z: and the pro- 
pagation function y2 of any TEM or TE wave 
in a dielectric are related in the following manner: 


(18a) 
where ye is the permeability of the dielectric 
medium. The permeability of dielectric materials 
is generally yo of free space; therefore 


¥2Z2=jwpe, 


¥2Z2 = Jwpo = Zi, (18b) 


and 


71 
Z(0) =Z,— tanh yed. 


010 4 s4 ec’ 
—=:25 
£5 
C 
= s 6.68 
ooe4 4 € 
Eum | *o ee ¥ \ 
Ewan | * Xe Yi \ 
x _ 
006 4 3 es « 
\ x 


002 4 


(19) 














Fic. 4. Emin/Emax and xo/Ai, versus the ratio of thickness 
d of the sample to wave-length A: in the dielectric. 
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Equating (19) and (14) yields 

















Ean 24Xo 
—jtan 
tanh vod —jM Ene Mi 
= . =Ce*%. (20) 
yd 2nd Ennio 24x 
i-j tan 
max Mi 


The function Ce can be found, as shown 
above, by measuring the thickness of the sample 
d, the wave-length \, in the air-filled pipe, the 
ratio of minimum to maximum field strength 
Emin/ Emax, and the distance xo of the first mini- 
mum from the dielectric surface. 

The next step is the determination of 7e** 
= 2d from charts or from a series approximation 
of the function 

tanh Te’" 


Te’ 


= Ces, 





(21) 


In Fig. 3 the argument r is the ordinate, the 
absolute value 7 the abscissa, while C and ¢ are 
parameters of intersecting curves. 7 and ¢ are 
expressed in degrees. The hyperbolic functions 
are multivalued so that measurements with a 
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Fic. 5. Constructional details of the apparatus. 


single thickness d of the dielectric may leave the 
value of y2 in doubt. But if two different thick- 
nesses are used, only one set of values T7/d and 
r, will satisfy both experimental results. 

For low loss materials it is desirable to use a 
thickness near \2/4 or higher odd multiples of 
d2/4. As Fig. 4 shows, Emin/Emax reaches maxima 
for these values of thickness. The apparatus used 
in 1940 could measure the loss if Emin/Emax 
were larger than 0.003. For d=)./4 in Fig. 4, 
the loss would be below this limit, while for a 
thickness d=3).2/4 or T=32/2, the loss could 
be determined easily. In handling dielectric 
materials with low loss one can obtain more ac- 
curate results than are obtainable from the 
chart, by using the following series approxima- 
tion, valid in the range C>1 when T3/2z. 


0.212 0.0096 
Ter j(4.71- ,: es — 





et), (22) 


The complex dielectric constant can be found 
‘from y2 by the general relation 


2r 2 

y= (“) —w*ye*. (23) 
Xe 

The cut-off wave-length \, is determined by the 
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geometry of the pipe and the mode of propaga- 
tion. The frequency of plane waves with a wave- 
length \, in the medium to which y refers is 
called the cut-off frequency. The values of \, for 
the lowest order waves in different wave guides 
are 


rectangular pipe’ A. =2X width 
round pipe*® A. = 1.71 Xdiameter 
coaxial line A, = 


The value of y: is given by 


> 


2r\? . 2r\* 
1’?= (~) — w* We = -(=) . (24) 
Ae A 


Bearing in mind that yu for dielectrics is generally 
uo Of free space, one obtains 


(1/A0)?— (y2/2m)? 
(1) re)? (1/1)? 





€o =€E (25) 


For the coaxial line this equation reduces to 
€2* = €9(y2/71)”. (26) 


The determination of ¢:* may now be sum- 
marized as follows. The ratio Enin/Emax, the 
distances x9 and d and the wave-length A; are 
measured. From these measured values C and 
¢ are calculated with the help of Eq. (20). From 
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Fic. 6. Details of magnetron. 


C and ¢ the values of T and 7 are found from the 
chart or by a series approximation. yz is found 


7L. J. Chu and W. L. Barrow, Proc. I.R.E. 26, 1520 
(1938). 

8G. L. Southworth, Bell Sys. Tech. J. 15, 284 (1936); 
W. L. Barrow, Proc. I.R.E. 24, 1298 (1936). 
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by dividing Te" by d. Finally e.* is determined 
from ye by Eq. (25) or (26). 


APPARATUS AND PROCEDURE 


The apparatus used in 1940 is shown in Fig. 5. 
It is now obsolete in many details and will 
therefore be described only briefly. 

A round wave guide P is mounted vertically 
and is fitted at the bottom with a detachable 
end piece in which either solid or liquid dielec- 
trics may be inserted. A magnetron oscillator T 
with its magnet M is located at the upper end 
of the pipe. The output of the magnetron is 
monitored by a bolometer detector B. The stand- 
ing wave above the dielectric is measured by 
means of the traveling crystal detector E which 
is driven by the lead screw W while its position 
is indicated to an accuracy of .001 inch by the 
revolution counter C. 

The magnetron oscillator is of a special design 
shown in Fig. 6. An open anode construction, 
which serves to radiate the energy into the pipe, 
avoids the tuned parallel transmission line which, 
in conventional constructions,’ feeds the energy 
through the glass envelope to a dipole antenna. 
The principal advantage of this design is its 
simplicity of construction. It has an efficiency of 
about one percent and requires a critical adjust- 
ment of the magnetic field intensity and the 
orientation in the field. 

The output of the magnetron is filtered and 
stabilized by means of a resonator R, a cylindrical 
tube fitting in the pipe and closed on both ends 
by diaphragms. The approximate dimensions of 
this filter for a wave-length of 6.0 cm are: length 
3 inches, diameter of tube 13 inches and diameter 
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Fic. 7. Bolometer circuit and method of construction. 


* E. G. Linder, Proc. I.R.E. 27, 732 (1939). 
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of hole in the diaphragm 3 inch. The resonator 
R is inserted for part of its length in the pipe 
section holding the magnetron. The main pipe, 
which slides over the free end of the resonator, 
is adjusted to resonance and then clamped in 
position. 





Fic. 8. Details of traveling crystal detector. 


The construction and circuit diagram of the 
bolometer detector B are shown in Fig. 7. This 
arrangement proved very stable and troublefree, 
but of course less sensitive than a crystal de- 
tector of the silicon-tungsten type. 

The slotted pipe maintains a wave which is 
polarized so that the electric field is most in- 
tense and the longitudinal magnetic field zero at 
the slot. The electric field is shown in Fig. 8 
together with a more detailed drawing of the 
detector construction. Good metallic contact be- 
tween detector housing and pipe is required and 
is provided by the slider Z held firmly with the 
phosphor-bronze spring A. 

In principle, one measurement with the di- 
electric material in the pipe would be sufficient. 
After tuning to resonance by adjusting the 
plunger D (Fig. 5), x0, Emin/Emax and d, are 
found by the rectified current and position of the 
traveling detector. In practice, a second meas- 
urement is carried through with the empty pipe 
for calibrating the detector. Here the field 
strength is known. It should go down to zero at 
the minimum and increase towards the maxi- 
mum proportionally to the sine of the phase 
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angle 6, where 


6 = 2x(x—xX0)/d. (27) 
Then 


| E(x): | = Emax|sin 6). (28) 


The calibration curve of the detector is found by 
noting the detector current for several positions 
on both sides of the minimum. The amplitude 
at the minimum is an indication of attenuation 
in the pipe and frequency instability of the os- 
cillator, which should be very low for proper 
operation. The wave-length is determined by 
measuring the positions of two minima. 

With a low loss dielectric in the pipe, the 
measurements are carried out only in the neigh- 
borhood of the minimum, because the maximum 
is flatter and the detector is likely to be over- 
loaded. In this case a simple formula can be de- 
rived for calculating the ratio of minimum to 
maximum field strength: The absolute value of 
electric field strength is found from (7). 


| E(x), | = \Aa| |1+e-2(nz+9) | 
=|Aa||1—e*%t|, (29) 


Then the ratio of this field strength to the mini- 
mum, given by (11), is 


|E(x):| |et#®—e-+i® | |sinh (p+ 76)! 


} e°’—e? | sinh p | 











sinh p cos 6+ cosh p sin 6} 





= . (30) 
| sinh p | 
The square of this ratio is 
| E(x), |? sin? 6 
—— = cos* 6-+—_—_. (31) 
Emin? tanh? p 


Solving for tanh p, one obtains the desired result. 


atin sin 6 
=tanh p= . (32) 


) | E(x), | 2 4 
( : — cos? ) 
Baia* 











~The use of this indirect method of determining 
Emin/Emax makes it possible to measure much 
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TABLE I. Dielectric constant and loss at =6.0 cm. 
(Temperature 25+2°C except where noted.) 














Material e /eo tan 6 
Solids 

Polystyrene 2.53 .0005 to .002 
Lucite 2.57 .005 to .008 
Plexiglas 2.60 -006 
Cellulose acetate 3.24 .033 
Polyethylene 2.25 .0002 to .0012 
Glass 

Soda-lime-silicate 6.7 .014 

Pyrex No. 774 4.8 .009 

Pyrex No. 707 3.9 .0012 
Fused quartz 3.80 .0001 
Paraffin wax 2.25 -0002 
Plywood (Bakelite bonded) - 1.9 0.060 
Poplar, grain || electric field 2.05 0.075 
Poplar, grain | electric field 1.75 0.065 

Liquids 

Distilled water (28.5°C) 71 0.22 
Glycerine U.S.P. (28°C) 5.6 0.60 
Benzene (thiophene free) 2.25 <.001 


smaller values of this ratio than could be done 
directly. 


RESULTS 


Table I gives the values of e’ and tan 6 at 
4=6.0 cm for a number of common materials 
including solids and liquids. The first measure- 
ments on most of these substances were made in 
1940 on the apparatus described above. Since 
that time extensive and more accurate data at 
A\=3 cm, \=10 cm and at longer wave-lengths 
have been obtained for a wide variety of di- 
electrics in this Laboratory and summarized in 
the two volumes Tables of Dielectric Materials 
prepared under O.S.R.D. Contract OEMsr-191. 
Table I has been edited to some extent on the 
basis of these more recent data. 
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Interpolation with the Aid of a Plot of First Differences 


Gorpon S. FULCHER 
Washington, D. C. 
(Received February 28, 1946) 


Given the values of y_1, yo, y1, and y2 corresponding to 
X_1, Xo, X1, and X2, we first plot Ay/Ax, the ratio of first dif- 
ference to abscissa interval, as a function of mean abscissa 
2, and then draw a smooth curve through the graph points 
which gives the values of a function S(2). To determine 
the value of yp for any value of x, between xo and x1, we 
substitute values of S(z’) and S(2’’), obtained from the 
curve, in the following equations: y’ = yo+fAx-S(2’); 
y” =n — (1—f)Ax-S(2"); yp =y' + (9 —y’)(1+f)/3; where 
2 =4(xotxp); 2 =H(xit+xp) and f=(xp—xo)/Ax. With 
a carefully drawn smooth curve, the error in the result is 
only about one-sixth of the fourth difference term in the 
Newton-Bessel interpolation formula. In a simpler varia- 
tion of the method, the function L(%), whose values are 
given by straight lines connecting the graph points, is 
substituted for S(#) in the equations; however, the ac- 
curacy of the result is only about one-sixth. It is shown 


HE author became interested in methods 

of interpolation while making ballistic 
calculations during the winter of 1944-45 
for the National Defense Research Committee. 
Accurate interpolation formulas have been avail- 
able since the time of Newton; but not being a 
mathematician, I have a preference for graphical 
methods; so I amused myself investigating the 
use of a plot of first differences in connection 
with interpolation. This led to the method de- 
scribed below which has proved useful when a 
number of interpolations are required in a single 
set of data, such as can be represented by a 
single curve, and also when the interpolation is 
between datum points that are unevenly spaced. 


TABLE I. Data used to illustrate the method. 





Datum Abscissa Ordinate Difference Mean x Graph 





point x y Ay Ay/Ax point 

a 0 0 
3698 739.6 a5 Py 

b 5 3698 
2646 529.2 7.5 P, 

c 10 6344 
2007 401.4 12.5 Ps; 

d 15 8351 
1636 327.2 17.5 P, 

e 20 9987 
1411 282.2 22.5 P; 

f 25 ~—- 11398 
1247 249.4 27.5 Ps 


g 30 
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that the accuracy obtainable in any case can be estimated 
from the sagittas o; and oo for £ equal to x; and Xo, respec- 
tively, where o,=L(%:)—S(#1) and o9=L(Z%o)—S(%o), 
since the error is proportional to (y’’—y’) times 
(o1—¢@09)/(01+¢0). This is true for unequal datum point 
intervals as well as for equal intervals. When the sagittas 
are equal or approximately so, the results may be as ac- 
curate if the rectilinear variation is used. Taking advantage 
of this fact, a diagram is given which enables the logarithm 
of any number to be determined to eight places by sub- 
stituting values of L(Z), obtained by linear interpolation 
between given figures, in certain equations. The method 
has been found useful when a number of interpolations are 
required in a single set of data. An important advantage 
is that it can be used as easily when the abscissa intervals 
between datum points are unequal as when they are equal. 


1. DESCRIPTION OF THE METHOD 


To explain the method we shall first apply it 
to an actual case with evenly spaced datum 
points, using the data given in Table I. To 
determine the values of y, for a series of values 
of x, between 10 and 20 we proceed as follows: 

Step one: plotting first difference as a function 
of mean abscissa. If the abscissa interval Ax 
between datum points is constant, first difference 
Ay may be plotted as a function of mean abscissa 
£; otherwise the ratio Ay/Ax must be plotted. 
Figure 1 shows a plot of Ay/Ax as a function of 
& for the graph points P; to Ps whose co-ordinates 
are given in columns 5 and 6 of Table I. The 
original plot is 10 by 12 inches. 

Step two: drawing smooth curve through the 
graph points. My procedure is first to draw a 
smooth curve through the points freehand, in 
pencil; next to draw through each point in ink 
a short segment of the curve with the help of 
French curves; and then to complete the curve 
between segments, again using French curves. 
Great care in drawing the curve is required only 
when maximum accuracy is desired. The smooth 
curve gives the values of a function of # which 
we shall call S(Z). 

These two steps complete the graphical part 
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of the method. The straight lines shown con- 
necting the graph points in Fig. 1 are used in 
connection with the rectilinear variation of the 
method which will be described later. 

Step three: making the calculations. Assuming, 
for example, that we wish the value of y for x 
equal to 14, we make the following three simple 
calculations. First, using y, as starting point, 
we compute what we shall call yi4s’ by adding 
(14—10) times the value of S() for Z equal to 
4(14+10), as read from the graph; then, using 
Ya as starting point, we compute what we shall 
call yiss’"" by subtracting (15—14) times the 
value of S(Z) for equal to $(15+14), as read 
from the graph; and finally we obtain a weighted 
average of the two approximate results by adding 
(1+f)/3 times their difference to yi4s’, where f 
is the ratio of (14—10) to Ax. We thus obtain 
_the following results: 


VYuas’ = 6344 +4 < 411.6 = 7990.4; 
Yus’ =8351—1X 367.6 =7983.4; 
Yus =7990.4+.6X (—7.0) = 7986.2. 


As we shall see, this result happens to be correct 
within .3 unit. The subscript S is used to dis- 
tinguish values obtained with the aid of the 
smooth curve from those obtained by use of the 
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straight lines connecting the graph points, as is 

done in the rectilinear variation of the method. 
In general, the equations used in calculating 

the value of y, corresponding to abscissa xp 

between x» and x; are: 

Ys =Yot (Xp—Xo)- S(Z’) = yot+ fAx- S(x’); (1) 


Yps’ =Yi— (x1 —x,)- S(Z”’) 


=yi—(1—f)Ax-S(x"); (2) 


Yos=Yos + (¥ps” —Yps')-(1+f)/3; (3) 
where 
E =} (Xpt+Xo); (4) 
£" =3(x1+<x,); (S) 
and 
f= (x»—Xo) (x, — Xp) a (X»p—Xo), ‘Ax. (6) 


To check the result obtained above for yi, we 
shall calculate the correct value by use of the 
Newton-Bessel formula! which may be written 


1The Newton-Bessel, Gregory-Newton, Lagrange, and 

other interpolation formulas are given in E.T. Whittaker and 
George Robinson's Calculus of Observations (Blackie and Son, 
London, 1940) and by other books on interpolation ; tables 
giving the numerical values of the coefficients for the Greg- 
ory-Newton and Lagrange formulas are included in Herbert 
B. Dwight’s Mathematical Tables (McGraw-Hill Book Com- 
pany, Inc., New York, 1941) and in H. T. Davis’ Tables of 
Higher Mathematical Functions (Principia Press, Inc., 
Bloomington, Indiana, 1935). 
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as follows: 


1 
Yo =2(Yoty1) +(F— 2) Ans — a f(l — f)(A*y0+A7y1) “— —f)(f—2)d*9 


1 1 
+g fl DUDA-N(atw+an) + f0-NF+DQ—NF- Darn 
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To apply this formula we need the differences 
given in Table II. Using these differences we 
obtain : 


Via = 7347.5+602.1+40.4 —2.144 
— 1.922-+.020 — .048 = 7985.91. 


The Lagrange seven-term formula gives 7985.92. 
Comparison of the value for yi4s obtained above 
with these correct values shows an error of +.3 
unit. Other values obtained with the aid of 
Fig. 1 are given in Table II] together with accu- 
rate values given by the formula. The figures in 
the last column show that the maximum error 
comes out less than .5 unit in this case, that is 
less than Ay/3000. 

To test the accuracy of the method further, 
four new plots of the data in Table I were made, 
and the smooth curves drawn through the graph 
points were used to repeat the calculations. For 
the 40 values thus obtained, including those in 
Table III, the average error came out less than 
4 unit and in only two cases was the error more 
than 1 unit. The time taken to plot the six graph 
points and to draw the smooth curve was less 
than ten minutes in each case and only two 
French curves, not particularly suitable, were 
used. 

An important advantage of the method is that 
it may be used when the datum points are not 


TABLE II. Extension of Table I, giving second to 











Aty Aby Aty z 
5 3698 — 1052 — 150 22 5 
413 5 7.5 
10 6344 — 639 — 145 18 10 
268 23 12.5 
15 8351 — 371 —122 14 15 
146 37 17.5 
20 9987 — 225 — 85 10 20 
61 47 240 
25 11398 


— 164 — 38 6 25 
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evenly placed as easily as when Ax is constant. 
Before giving an example of such use, however, 
we shall describe the rectilinear variation of the 
method which is somewhat simpler, since it does 
not require the drawing of a smooth curve, and 
yet is sufficiently accurate for many cases. 


TABLE III. Values obtained with the aid of Fig. 1. 











(Og Graphical Consent — 
ed Yes. Ys” —Y,8") io of ro Yes 
11 6789.2 6781.8 74 6786.2 6785.95 +.25 
12 7210.8 7199.6 11.2 7205.6 7205.4 +.2 
13 7610.9 7599.8 111 7605.0 7604.7 +.3 
14 7990.4 7983.4 7.0 7986.2 7985.9 +.3 
16 8703.8 8699.0 4.8 8701.9 8701.65 +.25 
17 9042.6 9036.0 6.6 9039.5 9039.35 . +.15 
18 9368.6 9363.0 5.6 9365.6 9365.4 +.2 
19 9683.8 9679.8 40 9681.4 9681.0 +.4 





2. RECTILINEAR VARIATION OF THE METHOD 


After plotting the graph points, we draw 
straight lines connecting them, as shown in 
Fig. 1. If we designate by L(Z) the function 
whose values are given by these lines, then the 
equations used in calculating ypz are the same 
as Eqs. (1) to (6) except for the substitution of 
L(z) for S(#). That is, instead of Eqs. (1) to (3) 
we use the following: 

Vou’ =Yot(xXp—Xo)-L(z’) =yot fax-L(x’); (8) 
Vp’ =y1— (x1 — Xp) L(2”’) 

=yi—(1—f)Ax- LZ"); (9) 

VoL =Ypot +(Ypt" — Yr’) A+f)/3, (10) 

where Z’, %’, and f have the same meanings as 

above. For example, to determine yi4 we make 
the following calculations: 


vist’ = 6344+4 X 414.0 = 8000.0; 
yu” =8351—1X 372.0 =7979.0; 
yiat = 8000.0+.6 X (— 21.0) = 7987.4. 
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The error of this result is +1.5 unit. Other 
results obtained with the aid of Fig. 1 are given 
in Table IV..The average error for these results 
is about six times the average error for the 
results given in Table III. As we shall see, this 


TABLE IV. Results obtained using values of L(z) 
from Fig. 1. 


Xp (¥, a YreL Error of ypz Est. Epz 
11 —22.4 6787.6 +1.6 42.3 — 
12 —330 7208.2 +2.8 +3.3 
13 — 33.0 7607.4 +2.7 +3.3 
14 —21.0 7987.4 +1.5 +2.3 
16 — 13.0 8702.8 +1.2 +1.6 
17 — 18.8 9041.5 +2.1 +2.3 
18 — 18.6 9367.7 +2.3 +2.3 
19 — 13.5 9682.9 +1.9 +1.6 
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is the ratio to be expected if the smooth curve is 
accurately drawn. 

This rectilinear variation does not require a 
plot since the values of L(Z’) and L(Z’’) can be 
calculated by use of the following equations: 


Ax- L(x’) =Ay,—3(1—f) -A*yo; (11) 
Ax- L(x’) =Ay,+3f-A?y:. (12) 


However, when a number of interpolations are 
to be made, a plot is useful. 


3. GRAPHICAL METHOD OF ESTIMATING 
THE ACCURACY OF RESULTS 


If the smooth curve is drawn and also the 
chords connecting the graph points, as in Fig. 1, 
the accuracy obtainable by this method of 
interpolation can be estimated as follows. The 
error of ypz is approximately given by the 
following equation : 


. 3 (71-90) ” ’ 
E»y.= ae ——(Vpz —VpL ), (13) 
8 (o1+90) 


‘where o; represents the vertical sagitta at %, 
that is the vertical distance between the chord 
P;P, and the smooth curve at € equal to x;, and 


go represents the corresponding vertical sagitta 
- at Z. Evidently, 


o,=L(#,;) — S(%:), (14) 
oo = L(¥o) — S(Zo). (15) 


and 


According to Fig. 1a, ¢20, 715, and o19 are equal to 
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about 2.3, 4.6, and 8 units, respectively. Using 
these values, we obtain for E,, the estimates 
given in the last column of Table IV. These 
estimates are seen to agree approximately with 
the actual errors given in the preceding column. 

If the smooth curve is accurate, the error of 
Yps will be about one-sixth of E,,. With an 
inaccurate curve, the ratio may be greater or 
less than six. In the case of Fig. 1, in spite of 
the fact that the values of S(#) given by the 
curve are from .3 to .8 unit too large, the errors 
in the values of yp s happen to come out somewhat 
less than if the smooth curve had been accurate, 
averaging +.27 instead of +.37 unit. (See Table 
X, below.) Thus, even with a hastily drawn 
smooth curve, the values of y,s obtained with 
the aid of it are likely to be several times as 
accurate as values of y,, obtained from the 
chords. 

An interesting case is when.the two sagittas 
are approximately equal. In this case, according 
to Eq. (13), the error of results obtained from 
values of L(#) will be very small, and there is no 
point in using the smooth curve except to 
ascertain that the sagittas are equal. Therefore 
when the sagittas are found to be nearly equal 
accurate values of L(#), calculated by use of 
Eqs. (11) and (12), should be used to determine 
Ypxt if maximum accuracy is desired. 


4. CASE OF UNEQUAL INTERVALS BETWEEN 
DATUM POINTS 


In the example used above to illustrate the 
method, the abscissa interval Ax between datum 
points is constant. We shall now apply the 
method to a case where the abscissa intervals 
are unequal, using the data given in Table V. 


TABLE V. Data used to illustrate the method as applied to 
unequally spaced datum points. (y=log x.) 


Mean x Graph 


x y Ay Ax Ay/Ax Ps point 
40 ~=:1.60206,0 

.04139,3 4 .01034,82 42 P, 
44 1.64345,3 

.05551,7 6 .00925,29 47 P2 
50 1.69897,0 

.06445,8 8 .00805,72 54 P3; 
58 1.76342,8 

.06908, 1 10 .00690,81 63 P, 
68 1.83250,9 

.07058, 1 12 .00588,17 74 P; 


80 =: 1.90309,0 
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Whereas when Ax is constant either Ay or 
Ay/Ax may be plotted as a function of mean x, 
when Ax is not constant it is necessary, of course, 
to plot Ay/Ax. Figure 2 shows such a plot of the 
data in Table V. The original is about 10 by 13 
inches. Using the smooth curve drawn through 
the graph points, values of yps were determined 
for integral values of x, between 50 and 58. The 
resulting values of (yps’’—yps’) and of Es are 
given in Table VI. The maximum error came 
out only 0.5X10- less than Ay/10,000. Another 
similar plot, made independently, gave the same 
maximum error. 

Much to my surprise, the results obtained by 
use of the chords between the datum points came 
out as accurate as those obtained with the 


TABLE VI. Results for log 51 to log 57 from data in 
Table V. (Unit is .00001.) 











(¥,5" —Yps) Es On" “Ten) Ep 
x (Fig. 2) (Fig. 2) (accurate) (accurate) 
51 —4.1 0 — 15.09 —.21 
52 —6.2 +.1 — 25.87 —.18 
53 —7.6 +.3 — 32.34 —.13 
54 —8.2 —.5 — 34.49 — .08 
55 —8.2 —.3 — 32.34 — .06 
56 —5.8 —.4 —25 87 — .06 
57 —3.1 0 — 15.09 —.10 
Average: —6.2 +.2 — 25.87 —.12 





58 62 66 
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smooth curve. Then, in place of the values of 
L(#) read from the chords, I tried values accu- 
rately calculated by use of Eqs. (11) and (12). 
The resulting values of (ypz’’—ypxr’) and of E, 
are given in Table VI. The maximum error 
came out only —0.21X10~, less than Ay/30,000. 
Then I calculated E, by use of Eq. (13), substi- 
tuting for o; and oo the values 4.3 and 4.1, 
respectively, as determined from Fig. 2. For x 
equal to 54 the estimated value for E, comes out 
(vpx"’—Ypr’)/112 or —0.31X10. This is as 
close to the actual value as is to be expected. 
As will be described later, investigation of 
other series of uneven intervals shows that with 
uneven datum point intervals as with constant 
intervals, Eq. (13) gives, at least approximately, 
the value of Ex, to be expected. From this it 
follows that if when interpolating in a certain 
interval, the adjacent intervals can be adjusted 
so that the two sagittas come out equal, then 
the interpolation results obtained by using accu- 
rate values of L(#) will be unusually accurate. 
This principle has been used to construct a 
diagram, see Fig. 5 below, which enables the 
logarithm of any number to be determined to 
six to eight places, from the data on the diagram, 
by interpolation, using the rectilinear method. 


621 




















1s T T L$ 
y z 
10 to 
Sr B 
Z 
0 ~ L 1 0 
Qo 5 10 = 1S 


x orx 


Fic. 3. Plot of a function of x and of S(Z), a function of 2. 


The same method can be applied to any function 
for which the values of Ay/Ax do not vary too 
widely, for example to sine @ or cosine @ or to 
tangent of 0° to 45°. 


5. INVESTIGATION WHICH LED TO THE 
INTERPOLATION METHOD 


We shall assume that we are given the co- 
ordinates of the datum points a, b, c, and d of 
the function of x plotted in Fig. 3 and that we 
wish to determine the value of y, at a point p 
with abscissa x». We plot the datum points, 
draw a smooth curve through them in order to 
locate » approximately, and draw the chords 
between datum points and finally the chords bp 
and pc. Now since 


Vp = Yo + (xp —X») + (slope of chord bp), 
and 


Yp =Ye—(Xe—Xp) + (slope of chord pe), (16) 


if we knew the slope of bp or of pe we could 
determine y,. It occurred to me to investigate 
whether these slopes can be determined from 
the known slopes of the chords between datum 
points. 

If the point p is moved along the curve from 
‘a to b, the slope of pb varies continuously from 
that of ab to that of the tangent to the curve 
at b; if p is then moved on to c, the slope of bp 
varies continuously from that of the tangent at 
b to that of the chord bc, while the slope of pc 
varies continuously from that of bc to that of 
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the tangent at c; and if p is then moved to d, 
the slope of cp varies from that of the tangent at 
c to that of cd. Meanwhile, the mean abscissa 
of the chord to or from p increases from }(x,.+xs) 
to 3(x.+x4). From this it is evident that the 
slope of the chord between p and one or other 
of the two nearest datum points, considered 
always as extending to the right, is a continuous 
function of the mean abscissa Z. If we designate 
this function by F(#), then Eqs. (16) become 


Yp=Yot+ (Xp»— Xp) - F(Z’) 


=F" (x. — Xp) ; F(z”), 


(17) 


where #’ is the mean of x, and x», and #” is the 
mean of x, and x,. 

In Fig. 3 the slopes of the chords ab, bc, and 
cd, that is the values of Ay/Ax for the three 
intervals, are plotted as a function of mean 
abscissa and a smooth curve is drawn through 
the graph points. If we designate by S(#) the 
function of # represented by this curve, then 
S(#) has the same values as F(#) for % equal to 
the mean abscissas of the chords ab, bc, and cd. 

To investigate the relation between S(2), 
whose values can be determined graphically 
with more or less accuracy from the data, and 
F(Z), whose values we wish for substitution in 
Eq. (17), the function y=log x was chosen for 
study since it is an unusually difficult one for 
accurate interpolation. Given the data included 
in Table VII, the first problem was to determine 
the relation between S(#) and F(#) for this 
function. : 

Using a seven-place table of logarithms, the 
values of F(#) were computed for all integral 
and half-integral values of between 45 and 70. 
The values so obtained were plotted in Fig. 4 


TaBLe VII. Data for the function y=log x. 











Graph 
point x y Ay/Ax 2 point 
a 30 1.47712,12 
.01249,39 35 P35 
b 40 1.60206,00 
.00969,10 45 Ps 
c 50 1.69897 ,00 
.00791,81 55 Ps 
d 60 1.77815,12 
.00669,47 65 Ps 
e 70  1,84509,80 
.00579,92 75 Pr 
f 80 — 1.90309,00 








JOURNAL OF APPLIED PHYSICS 




















3 


SY xt0* 


a 
c 
e 


z, 





























4 


Z V 
v4 
/ 
/ 


= 
~ 
xX 
~ 


“a 
vs 
/ 
| ‘ ry, Kk 4 & 









































and a dashed curve drawn through them. The 
original plot is 10 by 15 inches; and in order to 
increase the vertical scale, lines of constant 
ordinate were drawn in at 45° to the lines of 
constant abscissa. The curve for F(#) passes 
through the graph points, of course. 

A smooth curve through the graph points will 
obviously lie near the F(Z) curve. To investigate 
the relation between these functions, the smooth 
curve shown in Fig. 4 was drawn through points 
whose ordinates were calculated by use of the 
following equation : 


10-S(#) =log (+5) —log (—5). 


This was done to eliminate effects of inaccuracy 
in drawing the curve. Finally, the graph points 
were connected by straight lines which give the 
values of the function L(%). It then became 
evident that the S(Z%) curve lies everywhere 
between the functions F(%) and L(#) and the 
difference S(Z) minus F(#) is about one-third of 
the difference L(Z) minus S(Z) in this case. It 
was found that for between 50 and 60 the 
ratio of the differences is in fact close to 3 on 
the average, being equal to 3.12. This suggested 
that if we subtract from S(Z) one-third of 
[L(#) —S(%)] we shall obtain a function of # 
whose values are approximately equal to F(Z). 
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If we call this function C(#), then 
C(£) = S(@) —3LL(x) — S(x) J; 


and if we substitute C(#) for F(#) in Eq. (17), 
we shall obtain approximately correct values 
for Yp. 

Since the error in the result obtained is 
(X»—X») or (x-—Xp) times the error in C(Z’) or 
in C(%”’), we shall reduce the error on the 
average if we write two equations, as follows, 


(18) 


Yoo’ =Yo+ (xXp—xXv) C(#’) =o + fdx- C(x’); (19) 
Voce’ =Ve— (Xe — Xp) * C(x"’) 
=Ye—(1—f)Ax-C(x’’); (20) 


and use Eq. (19) for cases where f is between 0 
and 3, and Eq. (20) for cases where f is greater 
than 3. 

Results for log 51 to log 59, obtained by 
substituting values of C(%) from Fig. 4 in Eggs. 
(19) and (20), are given in Table VIII. The 
average error comes out about .5X10-, or 
Ay/16,000. This shows that the method is suffi- 
ciently accurate providing an accurate smooth 
curve can be drawn through the graph points. 

To determine whether the relation between 
the functions F(#), L(#), and S(#) suggested by 
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Taste VIII. Results for log 51 to log 59, obtained by 
substituting values of C(#) from Fig. 4 in Eqs. (19) and 


(20). (Unit is .00001.) 





xp z F(2) Cc) a! n° 
51 50.5 860.0 859.9 — I 
52 51 851.7 851.5 — 3 
53 51.5 843.5 843.3 —- 7 
54 52 835.6 835.4 - § 
55 52.5 827.8 827.6 —1.3 
55 57.5 755.8 755.9 —.7 
56 749.1 749.2 —.5 
57 58.5 742.5 742.6 —.2 
58 59 736.2 736.2 —.1 
59 59.5 729.9 729.9 0 


Fig. 4 is a general relation, equations for these 
functions in terms of the differences were derived 
from the Newton-Bessel formula, Eq. (7). From 
these equations, which are given in Appendix A, 
it follows that 


Ax -[C(#’) — F(2’) ] 
fA-f)At+f) 


i eeeemesienei (Atyo+A‘y,)+ete. (21) 
144 
Ax: [C(#”’) — F(#”’) ] 
Plt ial f) 
= +——_—_—__—-(Atyy+ Aty:) ete. (22) 
144 


In absolute value these differences each amount 
to less than one fourth of the fourth difference 
term in the Newton-Bessel formula, Eq. 7. In 
our example above, since (A*yo+A*y:) equals 
— .00702, the value of these differences for f=4 
should be about +1.8X10- divided by Ax, that 
is +.18X10~. The differences between columns 
3 and 4 of Table VIII agree with this result. 
When an inaccurate smooth curve is used in 
obtaining the values of C(#) for substitution in 
Eq. (19) or (20), however, any error in S(Z) is 
multiplied by (4/3)(x»—x») or (4/3)(x-—xp) so 
‘that in the case of Fig. 1, for example, the results 
for ypc came out with errors which average four 
times the average error for ypc when. accurate 
values of S(#) were used. (See Table X, below.) 
It occurred to me that if values of S(#) are 
substituted in Eqs. (19) and (20) for correspond- 
ing values of C(%), since S(Z’) is greater than 
C(z’) if S(#’) is greater than C(#’’), the errors 
in the two results obtained, y,s5’ and yps’’, will 
be opposite in sign; therefore, a weighted mean 
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of the two values might be fairly accurate. To 
test this, accurate values of S(%) from Fig. 4 
were substituted in Eqs. (1) and (2), and values 
for yps’ and y,s” were obtained for integral 
values of x, between 50 and 60. From these, the 
figures in column 2 of Table |X were obtained. 
Figures for the difference between y,s5’ and the 
correct values are given in column 3, and the 
ratios of these differences to the differences in 
the second column are given in column 4. A plot 
of this ratio as a function of f showed that the 
values run parallel to the,values of (1+/)/3 and 
differ from them by only about .067. It follows 
that, approximately, 


Yp=Yps +(Yps”’ —Yps )(A+f)/3. (23) 


If we call this result y,s5, then Eq. (23) is the same 
as Eq. (3); hence we have arrived at the method 
described at the beginning of the article. 

Equation (23) is equivalent to the following 
equation : 


1 (f-3) 1 (f-3) 
e| —_ |-5’+| + |p” (24) 
2 3 2 3 


The value yps is thus a weighted average of 
Yps’ and y,ys”. It was found that the error of 
this average is less, with an inaccurate curve for 
S(#), than the error of the corresponding result 
Ypc obtained from the function C(Z), on the 
average. This conclusion is illustrated by the 
figures given in Table X for the errors in the 
results obtained, by use of various equations, 
from the data in Table I and from Fig. 1. 
Although with an accurate smooth curve Ec¢ 
averages only half Es, with the actual curve for 
S(#), which gives values which are up to .8 unit 
too large, Ec averages three times Es. Of course, 


TABLE IX. Results calculated from S(#) in Fig. 4. 
(Unit is ease.) 








Xp (Ys. —Yps5) (y Se ’) Ratio of 
(1) (2) (3) (3) to (2) 4f+1) 
51 — 6.45 —2.80 434 367 
52 —11.25 —5.27 468 .400 
53 — 14.49 —7.28 .502 433 
54 — 16.26 — 8.70 .535 467 
55 — 16.62 —9.44 .568 .500 
56 — 15.67 —9.40 .600 .533 
57 — 13.47 —8.51 .632 .567 
58 — 10.08 — 6.68 .663 .600 
693 633 


59 — 5.57 
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this is only one case. However, my experience 
shows that Eqs. (1) to (3) will give more accurate 
values than Eqs. (19) and (20), as a rule. 
Nevertheless, since only one of the two Eqs. (19) 
and (20) need be used in any particular case, 
it may be convenient to use the C(#) function 
rather than the S(%) function when this gives 
sufficient accuracy. 

Table X also gives the values of E,, the errors 
of results obtained by substituting values of 
L(#) in Eqs. (8) to (10). When accurate values 
of L(#) are used, 


_fA-NG+N2-f) 





~L= 
48 
X (Atyo+A*y:) —etc.; (25) 
and since 
fai-f) 
(ypx"’ —Ypt') = -———_A'*yj; (26) 
and 
(o,1—@0) 1 (A*yo+Aty;) 
- —etc., (27) 





(o,+00) 4 A*y, 

where o; and og» are the vertical sagittas defined 
above by Eqs. (14) and (15) 
mately, 


, therefore, approxi- 


3 (o,;—a0) 
a a ew 


8 (a, +09) 





You’ —Ypt’). (28) 


This is the same as Eq. (13) above. In the case 
of Fig. 4, o1 is equal to 51.33 and oa» to 89.55; 
hence E,, should be equal to about .1017 X (ypi” 
—Ypr') for values of xp between 50 and 6.0 


TABLE X. Errors of results from data in Table I, (1) 
using values of the functions from Fig. 1, and (2) using 
accurate values. 


(1) Values from as iii 1 (2) Accurate values 


Xp Ex Ec E, ‘Ss Ep Ec Es 
11 +1. 6 + 8 +. 3 +2.00 —.15 +.34 
12 +2.8 +1.2 +.2 +307 -—.3 +.50 
13 +2.7 —2.2 +.3 +3.03 —.25 +47 
14 +15 —-10 +.3 +196 —.1 +.33 
16 +1.2 + .7 +.25 +155 —.1 +.27 
17 +2.1 + 8 +.15 +240 —.25 +.42 
18 $2.3 — 4 +.2 +2.30 —.15 +.40 
19 +19 — 5 +.4 +147 -—.05 +.22 
a +20 +9 +.3 +2.22 —.17 +.37 
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TABLE XI. Results for log 54 interpolated in several series 
of datum points with various intervals. 








EuL 








Abscissas of Series of Eur 3 (1-00) Eus 
datum points intervals X10 Yi” tyr) Seite) X10 
42, 45, 50,58, 70,88  3-5-8-12-18 -—2.08  +.061 +.080 —3 
40, 44, 50, 58,68,80 4-6-8-10-12 -0.08  +.002 +015 +41 
( Figure 2) 
37, 43, 50, 58, 67,77  6-7-8-9-10 +146 —.041 —.030 +3 
36, 44, 50,58, 66,74 8-88-88 +3.10  —.085 —.085 +.5 
40, 50, 58,64,58 12-10-8-64 4640 —.165 —.158 +.8 








Actually the average ratio of Ex, to (ypr’’ —Yyp1’) 
comes out .101. 

With an accurate smooth curve, Es is about 
one sixth of E,. According to the figures in 
Table X, the ratio is one to six for accurate 
values and about one to six for values from Fig. 1. 

In the case of unequal intervals between datum 
points, after the results given in Table VI had 
been obtained, various other series of unequal 
intervals were investigated. Accurate values of 
L(&) were used. The results for E, for f equal to 
3, that is for log 54, are given in Table XI; also 
values of the ratio of Esa, to (¥sax’ —Ysar’) and 
of 3(¢:—«¢0)/(o:1t+00). The agreement between 
the figures in columns 4 and 5 of Table XI 
indicates that Eq. (28) gives the value of Ex 
whether the intervals between datum points are 
equal or unequal as accurately as the ratio of 
(¢:—o9) to (o;+09) can be determined graphi- 
cally. 

The values of Es given in the last column of 
Table XI were calculated by use of as accurate a 
smooth curve as could be obtained but are only 
approximately correct. They indicate that the 
ratio of E, to Es is practically equal to six 
except when the errors approach zero. The 
general conclusion is that the greater the value of 
the ratio of (a:—o9) to (¢:+¢09) in absolute 
value, the greater the advantage in accuracy of 
using the smooth curve rather than the chords 
between graph points in interpolation by this 
method. 


6. DIAGRAM WHICH ENABLES ACCURATE VALUES 
TO BE OBTAINED BY LINEAR 
INTERPOLATION 


For a certain interval, say log 50 to log 60, 
accurate values can be obtained by linear inter- 
polation between graph points, as we have seen, 
if adjacent datum intervals are correctly chosen. 
Now the effect of changing these intervals is to 
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Fic. 5. Values of 10L(2’) and of 10L(2’’) and other data for 
use in calculating log x. 


change the slopes of the lines connecting the 
graph point at the center of the interval with 
the adjacent graph points. These slopes are equal 
to [L(#),—L(#),]/}4x and [L(#)o—L(#),; ]/4Ax, 
where L(Z),, L(#),, and L(#)o represent the 
values of L(%) for & equal to x, x, and Xo, 
respectively. Given the value of L(%),;, we can 
adjust the values of L(Z), and of L(%)o so that (1) 
the error of yj, is zero, and (2) so that the error 
for f equal to some value, say .2, is the same as 
the error for f equal to say .8. The equations 
for L(#),; and L(#)»o are then as follows: 


Ax[L(#):+L(#)o] 

(y.s—Y.2) —(.8—.2)-Ayy 
2X .8X.2X(.8—.5)/3 

Ax[L(#)1—L(#)o]=4(y1+y0— 29); 


from which 2L(#); can be obtained by adding 
Eqs. (29) and (30), and 2L(%)o by subtracting 
Eq. (30) from (29). Substituting log 50, log 52, 
log 55, log 58, and log 60 for yo, y.2, i, ¥.s, and 
_¥1, respectively, we obtain: 


= 2Ay,— 





(29) 


(30) 


10[L (2): +L(#)o] =.15836,25006 — (.04742,46500 
— .04750,87476) /.032; 
= .15836,25006+.00262,8050 
; = .16099,05506. 
10LL(&)1—L(@)o]=4(3.47712,12546 
— 3.48072,53790) 
= —.01441,64976. 


Hence: 


L(#), = .00732,87026; and L (#)9=.00877,03524. 
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If we substitute these values in the equations: 
L(x’) = L(&),+(1—f)[L(4)o— L()4]; 
L(x") = L(x),— fLL(@),— L(@)1); 


(31) 
(32) 


and substitute the values of L(z%’) and L(#’’) so 
obtained in Eqs. (8), (9), and (10), we shall 
obtain values of log xp, for values of x, between 
50 and 60, with a maximum error of .00000,18; 
and if we substitute for Eq. (10) the following: 


Vp= Vp 4 + (yp - — YpL) 


1 (f—}) 
+ a(f—ayef, (33) 
2 3 


where in this case A has the value 0.0456, the 
maximum error will be reduced to about 
0.00000,005. 

The constants in Table XII enable the loga- 
rithm of any number to be calculated. The 
result will be accurate to from five to seven 
places, depending on which interval is used, if 
Eq. (10) is used, and to two places more if 
Eq. (33) is used. For example, in computing 
log 48 we get yas’ = 1.68169,9564; yas’ = 1.68092,- 
9260; (yas'’ — Vas’) = — .00078,0304; Va8s= 1.68123,- 
7382, if Eq. (10) is used, and yyg= 1.68124,1243 
if Eq. (33) is used. The errors of the two results 
are —.3856X10-> and +.0006X10-, respec- 
tively. If log 48 is obtained by computing log 96 
and subtracting log 2 the corresponding errors 
come out —3.36X10-* and —0.16X10-8, re- 
spectively. Thus, by using the interval 80 to 90 
or 90 to 100, logarithms accurate to eight places 
of decimals can be obtained from the data in 
Table XII. 


TasLe XII. Constants for calculating log x by use of 
function L(z).* 








yo and y1 


10L(Z)o 


xo 10L(%)4 10L(Z)1 A 
40 1.60205,9991 -11010,7277 
.09691,0013 0557 
50 1.69897,0004 .08770,3524 .08852,7149 
.07918,1246 0456 
60 1.77815,1250 .07289,6175 .07328,7026 
.06694,6790 0385 
70 =: 1.84509,8040 .06237,7539  .06258,6485 
.05799,1947 .0335 
80 1.90308,9987 .05451,7820 .05463,9540 
.05115,2522 0295 
90 1.95424,2509 .04842,0736  .04849,6400 
.04575,7491 .0263 
100 2.00000,0000 .04360,1931 











® The equations to be used in calculating log x are: Eqs. (31), (32), 
(8), (9), and either (10) or (33). 
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Although a graph is not necessary, such a diagram as that in Fig. 5 is useful as a guide to the 
algebraic signs to use in computing L(%’) and L(#’’); moreover, it gives all the data included in 
Table XII. 

APPENDIX. EQUATIONS FOR FUNCTIONS OF MEAN X. 


The datum points are equally spaced in all cases; and the functions are all accurate. The 
equations are derived from the Newton-Bessel formula, Eq. (7) above. 
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4 4 
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(1—f) (i—fP 
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The Use of Color and Fluorescence Indicators for Determining the 
Structure of Glasses 


W. A. WEYL 
Glass Science, The Pennsylvania State College, State College, Pennsylvania 
(Received April 25, 1946) 


The fact that light absorption and emission depend on the environment of the ions and 
molecules which absorb or emit light can be used to draw conclusions about the structure of 
liquids, glasses, and crystals. This method supplements, to a certain extent, information gained 
from other sources, such as x-ray diffraction data. The method is illustrated by several 
examples, like the influence of thermal history on the structure of glass and the constitution 


of the yellow fluorescing zinc orthosilicate. 


I. INTRODUCTION 
JR a long time chemists avoided studying 
macromolecular substances. The classical 
period of organic chemistry with its brilliant 
synthesis of natural and artificial dyes of alkaloids 
‘and pharmaceutical compounds left the field of 
high polymers practically untouched. The chem- 
ist of this period neglected all those reaction 
products which had the slightest resemblance to 
- a resin. Reactions leading to an amorphous 
product which could not be purified by crystal- 

lization were considered failures. 
This hostile attitude to everything with a 
high molecular weight changed only when the 
physicist began to cooperate with the chemist 
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and when physical measurements became routine 
methods in the chemical laboratory. With the 
increase of the molecular weight certain physical 
properties become more characteristic, whereas 
the strictly chemical features of the macro- 
molecules lose some of their significance. It 
becomes difficult to distinguish similar macro- 
molecules of different degrees of polymerization 
by purely chemical reactions. The distinguishing 
features are their diffusion speed, viscosity, 
cohesive forces, and their relaxation times. This 
is true for both organic and inorganic com- 
pounds. The work which shall be discussed here 
was chiefly designed to study inorganic high 
polymers, especially silicate glasses. 
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The vitreous state of matter has certain char- 
acteristics which are brought about by the high 
degree of polymerization and by the fact that 
the relaxation time of glasses is too high to 
allow the establishment of an internal thermo- 
dynamical equilibrium. Textbooks usually em- 
phasize that glasses are supercooled melts 
which are instable with respect to their crystal- 
line phases. This is no doubt true for many cases 
and it explains why certain glasses tend to 
devitrify, that is, to form crystals. For the 
understanding of the typical glass properties, 
however, it is much more important to realize 
that glasses possess an atomic configuration 
which had been stable at a certain higher 
temperature and which on its way to room 
temperature has been ‘‘frozen in’’ because of the 
high relaxation time of the glass forming units. 

Water can be cooled below its freezing point. 
The water is then a supercooled liquid but not a 
glass, it is instable with respect to its crystalline 
form, the ice, but it is in internal equilibrium. 
No matter how the supercooled water was 
prepared, its properties are well defined if the 
variables, temperature and pressure, are fixed. 
For glasses, however, this is not true. Depending 
on its thermal history, the properties of a glass 
of given composition may vary widely. The 
density, refractive index, chemical resistivity, 
or the electric conductivity of a glass is no 
longer defined by its composition, temperature, 
and pressure. 

The thermal expansion of a supercooled liquid 
is uniquely defined by a curve which gives the 





° 200 4006 «600 0 20 400 600 
Temperature in °C 


(a) (b) 


Fic. 1. Schematic curve of the thermal expansion of 
(a) an annealed glass and.(b) a rapidly chilled glass. 
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volume as a function of the temperature. If we 
allow a borosilicate glass to cool very slowly to 
room temperature and then determine its ther- 
mal expansion, we find a curve like that in Fig. 
la. If the same glass was chilled rapidly, so that 
the polymerization is incomplete, it shows an 
entirely different thermal expansion. As indi- 
cated in Fig. 1b the expansion at low temperature 
changes into a contraction around 350°C. This 
contraction is caused by a polymerization process 


LENGTH 


LENGTH 














200 400 
TEMPERATURE °C 





Fic. 2. Space diagram, showing the relationship between 
volume, temperature, and previous heat treatment. 


or by a change in the geometrical arrangement 
of the glass forming units which could not take 
place on cooling because of the high relaxation 
time and the speed of passing through this 
temperature region. 

In order to describe fully how a glass of given 
composition will change its volume on heating 
at a certain rate, we need a space diagram which 
represents the volume as a function of both 
temperature and previous heat treatment (Fig. 
2). 

Based on economic and technological consider- 
ations such as the source of raw materials, the 
methods of manufacture, and the processes used 
for forming the material into its final shape, one 
has repeatedly tried to draw a sharp demarcation 
line between inorganic glasses and amorphous 
organic materials. Modern developments such 
as glasses free of silica (phosphate glasses) or 
even free of oxides (beryllium fluoride) and 
especially the organic derivatives of silica make 
a sharp distinction very difficult. From a scien- 
tific point of view there was never a justification 
for a distinction between organic and inorganic 
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glasses. The vitreous state with its characteristics 
such as optical isotropy, typical thermal expan- 
sion and its dependence of all properties on the 
previous thermal treatment is found among 
elements, oxides, sulfides, salts, and a great 
variety of organic compounds. G. Tammann! 
to whom we owe much of the early work on 
glasses and their properties studied such com- 
pounds as brucin, salicin, betol, piperin as well 
as borates, phosphates, and elemental selenium. 
These glasses can be studied within a more 
convenient temperature range and their principal 
features are not different from those of the high 
melting silicate glasses. 

Just to give you another example of the close 
analogy of inorganic and organic high polymers, 
I want to mention that recently we synthesized 
a purely inorganic material with rubber elastic 
properties, a hydrated alkali-manganese-phos- 
phate. This was accomplished by cross-linking 
the long chain-like sodium-potassium-polymeta- 
phosphate molecules by means of a two-valent 
ion, in other words, by means of a unit having 
two functional groups. Sodium-potassium meta- 
phosphates can be obtained as high polymers 
with molecular weights of 100,000. Their aqueous 
solution forms a viscous liquid. If this solution 
is allowed to react with divalent ions, in partic- 
ular with manganese, a complex alkali-manganese 
phosphate is precipated. The divalent ions link 
together the long chains thus depriving them of 
their free mobility. The resulting rubber-like 
product shows birefringence when stretched. 

The polymerization of P,O,; can still be studied 
by similar methods as that of organic compounds. 
Despite its tendency to form high polymers, 
P.O; is soluble to a certain extent in such solvents 
as chloroform. SiO, glass as compared with PO; 
glass, however, is insoluble in all solvents. The 
high viscosity of fused silica indicates that even 
‘at 1800°C this compound is still highly poly- 
merized. j 

The methods for studying the structure of 
vitreous silica and silicates are greatly limited 
- by the insolubility of these high polymers and 
by the high temperatures required to decrease 
their relaxation times sufficiently to bring them 
within experimental range. So far the Fourier 


*G. Tammann, (a) Kristallisieren und Schmelzen (Leip- 
zig, 1903); (b) Der Glassustand (Leipzig, 1933). 
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analysis of the x-ray diffraction patterns has 
proved to be the most powerful tool for studying 
the structure of simple glasses. If the number of 
the atomic species does not exceed 3 or 4, one 
may obtain fairly reliable data concerning the 
average interatomic distances and the coordina- 
tion number of the cations. 

B. E. Warren? and his students have perfected 
this method and were able to derive the atomic 
structures of several simple glasses such as SiO», 
BeF.2, BO 3, and alkali borates and silicates. It 
must be considered as one of the major results 
of their work that they found the laws of crystal 
chemistry to hold true for the vitreous silicates 
as well as for the crystalline ones. 

Before it became possible to interpret the 
x-ray diffraction patterns of amorphous sub- 
stances, the author® and his collaborators used 
color and fluorescence indicators to obtain insight 
into their atomic structures. This method is 
based on the introduction of certain ions of the 
transition elements into the glass. Their light 
absorption and fluorescence is affected by the 
surrounding ions. Additions of new ions, substi- 
tutions, and changes in the structure due to 
temperature and heat treatment reflect on the 
indicator ions which in turn signal this change of 
their environment by changing their color or 
fluorescence. 


Il. THE USE OF COLOR AND FLUORESCENCE 
INDICATORS 


Changes of color have been widely used to 
indicate changes in chemical and physical condi- 
tions., The acid-base indicators such as methyl 
orange and phenolphthalein are well known. 
Cobaltous chloride has been used to indicate the 
humidity of the air. Since the use of silver 
mercuric iodide to indicate temperatures above 
35°C, by its transition from the red to a yellow 
modification, many color indicators have been 
developed which permit us to estimate tempera- 
tures of moving machine parts or furnace walls 
over a wide range. 





2B. E. Warren, J. App. Phys. 8, 645 (1937). 

3W. A. Weyl, Verdffentl. Kaiser Wilhem-Institut fiir 
Silikatforschung 7, 167-203 (1935). A rather complete 
bibliography on the subject of color and fluorescence 
indicators can be found in the author’s monograph on 
Colored Glass prepared for the British Society of Glass 
Technology. Part I. The Constitution of Colored Glass, 
J. Soc. Glass Tech. 27, 133-206 (1943). 
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WAVELENGTH IN MILLIMICRONS 


Fic. 3. Light transmission of four alkali silicate glasses 
containing neodymium as a colorant. 


The polymerization of anions to iso- or hetero- 
poly acids is associated with a change in light 
absorption. Studies in the ultraviolet and in the 
visible region revealed that the formation of 
polymer tungstates, vanadates, or chromates 
causes the absorption band to shift to longer 
wave-lengths. For example, the addition of 
sulphuric acid to a neutral chromate causes its 
yellow color of the CrO,?- to change into the 
orange of the Cr,O;*- and finally into the deep 
red of the Cr,Q,;?-. 

An analogous shift of the absorption edge 
with increasing molecular weight has been ob- 
served with alkali polysulfides. Both cases, 
chromates and polysulphides, have been studied 
in aqueous solutions as well as in silicate glasses. 

F. Weidert was the first to introduce color 
indicators into the glass melt. When studying 
the light absorption of glasses containing the 
rare earth element neodymium, he found that 
its spectrum is affected by the amount and the 
nature of the alkali present in the glass. Figure 3 
gives the light transmission between 550 and 
650 my for four alkali silicate glasses containing 
neodymium. The fine structure of the neo- 
dymium band is most pronounced in the rubid- 
ium glass. Replacing the large rubidium ion by 
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other monovalent ions of decreasing size causes 
the fine structure to disappear. The small 
lithium ion causes the strongest perturbation. 
The disturbance caused by the electric fields of 
the alkali ions expresses itself in many different 
ways. The stronger the electric field of the alkali, 
the more fluid is the melt and the greater is the 
tendency of the glass to devitrify. Some potas- 
sium silicate glasses require many months of 
heat treatment to produce crystals. Lithium 
silicate glasses are so unstable that only small 
quantities can be obtained without crystalliza- 
tion. Sodium silicates assume an intermediate 
position. Weidert showed that for a large number 
of glass compositions the fine structure of the 
neodymium band can be used as a means to 
determine the stability of a new glass composi- 
tion. Blurring of the fine structure was found to 
be an indication of pronounced tendency to 
devitrify. 

The light absorption of neodymium is caused 
by electron transitions taking place in rather 
well-protected inner orbits. The electrons which 
are responsible for the color are not identical 
with those which are involved in the formation 
of chemical compounds. In contrast to the 
colored salts of the elements of the transition 
group, all compounds of neodymium have nearly 
the same pink color. The oxide, the hydroxide, 
the sulphate, and the ammoniates have nearly 
identical absorption bands. Only minor shifts in 
the maxima and in the fine structure indicate 
that the deformation of the orbits of the outer 
valence electrons affects to a certain extent the 
incomplete inner shells. 

The author and his student used preferably 
those color indicators whose light absorption is 
caused by the outer or valence electrons. 

The light absorption of the salts of the ele- 
ments of the transition group such as Mn, Fe, 
Ni, Co, Cu, Cr, or V is caused by electron 
transitions prohibited in the free ions but made 
possible by a deformation caused by surrounding 
ions (Molecular Stark effect). The maximum 
molar extinction coefficient of these ions is there- 
fore usually low and of the order of 0.1 to 1.0. 
Under certain conditions these ions form com- 
plexes with much stronger light absorption which 
can be attributed to electron transfer bands. 

There is no fundamental difference between 
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Fic. 4. Light absorption of quinizarin in (a.) hexane, 
(b.) hexane-nitrobenzene (7 percent), and (c.) hexane- 
nitrobenzene (24 percent). 


the light absorption of these ions, in aqueous 
solutions and in glasses. The rigidity of glasses, 
however, prevents their frequent collisions with 
their neighbors and as a result enables many of 
them to fluoresce. None of the simple ions such 
as Mn?**, Cut, V5* fluoresces when dissolved in 
water, but they do so when introduced into 
glasses. The fluorescence—like the absorption 
spectrum—is also influenced by the environment 
as will be discussed later for the case of the 
divalent manganese ion. 


Ill. THE SPECTRUM OF INDICATORS AS 
AFFECTED BY SOLVATION 

The strength of the forces which different 
solvents exert upon the solute can vary very 
widely. They are weak if the solvent consists of 
molecules which have no permanent dipole and 
low polarizability. Such is the case for hexane 
and other “inert” solvents. Molecules with ‘a 
greater dipole moment and greater polarizability 
exert stronger forces. These forces can be com- 
pared with external electric fields which are 
asymmetrical and which constantly change their 
direction with the thermal motion of the mole- 
cules. The splitting of spectral lines in external 
electric fields is well known. In case of asym- 


metrical fields a blurring of the fine structure 


results rather than a sharp splitting of spectral 
lines into their components. Change in direction 
of the field, as caused by the thermal motion of 
the surrounding molecules, acts in a similar way. 

Many organic compounds which are soluble 
in non-polar hexane have an absorption spectrum 
which exhibits some fine structure. Figure 4a 
gives the light absorption of quinizarin in hexane. 
The two maxima are most pronounced in non- 
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polar solvents. Addition of 7 percent of the 
strongly polar nitrobenzene (Fig. 4b) blurs the 
sharp maxima and in a mixture of hexane and 
24 percent nitrobenzene they have completely 
disappeared (Fig. 4c). 

The light absorption of quinizarin has here 
been used as an indicator for the polarity of the 
nitrobenzene. The electric fields of the nitro- 
benzene molecules exert an influence which can 
be compared with that of the lithium in neo- 
dymium glasses (Fig. 3). 

The absorption spectra of many crystalline 
compounds such as chromium salts resemble 
those of glasses containing the same color centers, 
but they differ in their behavior on cooling. The 
broad absorption bands of the green chromium 
glasses do not change much in their appearance, 
if the glass is cooled to the temperature of 
liquid air. The crystalline chromium salts in 
this temperature range possess fine structure 
resembling line spectra. Why do the spectra of 
crystals and of solutions or glasses respond to 
temperature changes in a different way? In both 
cases the transitions of electrons in the chromium 
ion are affected by the external electric fields of 
the surrounding ions. The lack of symmetry of 
these fields in crystals is caused by thermal 
vibrations. As the temperature approaches abso- 
lute zero the intensity of the vibrations decreases 
and the color centers find themselves in electric 
fields of identical symmetry which now do not 
change any more with time. In glasses the alter- 
nating electric fields produced by the surrounding 
ions in thermal motion become also more con- 
stant on cooling, but the glass structure causes 
each color center to be in a slightly different 
environment. As the light absorption integrates 
over a multitude of color centers, a blurred 
spectrum results even if the thermal motion has 
ceased. 

This behavior excludes all those pictures of 
the glassy state which are based on a heapwork 
of submicroscopic crystals. The state of disorder 
must be in molecular dimensions. 

Besides its influence upon the fine structure, 
the solvent determines the exact position of the 
light absorption. Increased solvation often shifts 
absorption to shorter waves, indicating an addi- 
tional energy requirement for the electronic 
jump. Perhaps the best known example is the 
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color of elemental iodine in various solvents. 
Hexane as an inert solvent allows the iodine 
molecules to display optical properties similar 
to those in the vapor state. Solutions of this 
type are purple, like iodine vapor, and are 
characterized by a strong absorption band in the 
green region of the spectrum. Addition of mole- 
cules of higher polarity or greater polarizability 
to the hexane solution causes this band to 
broaden slightly and to shift to shorter wave- 
length (Fig. 5). The electronic transition which 
is responsible for the light absorption in the 
green region requires now a greater energy. For 
a theoretical interpretation of this shift we might 
use the same picture which D. M. Bose and S. 
Datta developed for explaining the different 
position of absorption bands and the variation of 
the magnetic susceptibility of an ion in different 
lattices. They assumed that a certain interpene- 
tration of the outer electrons of the paramagnetic 
ion and its environment has to be disentangled 
before the electron transition can take place. 
The energy required to loosen the ‘‘1-coupling”’ 
reflects on the position of the absorption band. 

Changing the solvent from hexane to benzene 
shifts the maximum absorption of iodine from 
525 mp to 500 mu. Using toluene instead of 
benzene causes a similar but greater shift of the 
absorption band to 485 my. The iodine can be 
used as an indicator for the solvation forces 
exerted by benzene, toluene, or similar com- 
pounds. The further the band is shifted to shorter 
waves, the stronger is the interaction between 
solvent and iodine, probably because of greater 
mutual deformation. 

If the solute can be excited to fluoresce, we 
find a similar influence of the solvent on its 
light emission, both on the fine structure and on 
the position of the bands. 

Addition of nitrobenzene to a solution of a 
fluorescing organic compound in hexane quenches 
the fluorescence. Even the strong fluorescence of 
mineral oil can be decreased to a small fraction 
of its original intensity if a few drops of nitro- 
benzene are added. Under standardized condi- 
tions, using a solution of nitrobenzene in hexane, 
one can use this effect as the basis for working out 
quantitative methods for fluorometric analysis. 

In glasses elemental selenium produces a pink 
color; its absorption band broadens if potassium 
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is replaced by sodium or lithium. The same 
change of composition causes the red fluorescence 
of the selenium to become fainter and to dis- 
appear.* 

More spectacular than the quenching of fluo- 
rescence in glasses and solutions is the shift of 
the emission spectrum. 

In the case of iodine we have noticed that 
increased solvation calls for larger light quanta, 
thus shifting the absorption to shorter waves. 
As far as the fluorescence is concerned the 
increased solvation expresses itself in the emis- 
sion of smaller light quanta and consequently in 
a shift of the fluorescence color to longer waves. 
Using uranium glasses as an example, one finds 
that the uranyl group can produce a group of 
five relatively sharp emission bands in the blue- 
green part of the spectrum (potash glasses, 
certain phosphates, and borates) and as the 
other extreme a broad emission band with little 
or no fine structure ranging from green to orange 
(lithium or high sodium silicates). 

The organic chemistry offers numerous ex- 
amples for the variation of the fluorescence color 
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Fic. 5. Light absorption of iodine in hexane-benzol. 


eS suk et Pe 
Curve 1: Hexane pure. 


Curve 2: Hexane 14 percent benzol. 
Curve 3: Hexane 28 percent benzol. 
Curve 4: Hexane 54 percent benzol. 
Curve 5: Benzol pure. 


* The fluorescence of cations in silicate or phosphate 
glasses can be improved if the O-~ ions are partly replaced 
by F~ ions. Introducing the more deformable Br~ and I- 
ions—even in small amounts—on the expense of the less 
polarizable O~— quenches the fluorescence. 
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TABLE I. 


Diaminoterephthalic acid ester 
Solvent 


para meta 
Hexane green blue 
Benzene green-yellow green-blue 
Ether green-yellow green-blue 
Chloroform green-yellow green 
Acetone vellow green 
Pyridine orange-yellow green 


Ethyl alcohol orange vellow-green 


with the nature of the solvent. The most striking 
shift has been reported for the dimethylnaphth- 
eurhodine. In hexane this compound produces a 
strong bluish green fluorescence. With increasing 
polarizability and dipole moment of the solvent, 
the fluorescence is shifted to green, yellow, orange 
and finally becomes red in those solvents which 
contain strong polar groups, such as the alcohols. 
One is easily tempted to give a more specific 
interpretation of this phenomenon and correlate 
the maximum of the emission bands with the 
dielectric constant or the dipole moment of the 
solvent. One should keep in mind, however, that 
the deformation or interpenetration of the elec- 
trons responsible for the change of light absorp- 
tion and emission depends upon a number of 
factors among which the permanent and induced 
dipole moments, the polarizability, size, shape, 
and the position of the responsible atomic group 
within the molecule are the most important. 

The interaction between solvent and solute 
can be intensified by introducing polarizable 
groups into the solvent, as can be seen from the 
influence of benzene and toluene on the light 
absorption of iodine. It can also be changed by 
introducing or varying the polar groups of the 
solute. 

Table | influence which the 
solvent has on the fluorescence of the esters of 
para- and meta-diaminoterephthalic acids. Table 
-II describes the change of fluorescence of the 
para diaminoterephthalic acid ester dissolved in 
ether and in alcohol if the amino groups are 
made inert by acetyl groups and partly or fully 
_ replaced by hydroxyl groups. 

Increasing the forces acting between a solute 
and its solvent shifts the light absorption to 
shorter and the light emission to longer wave- 
lengths, no matter how this increase in the 
mutual deformation and attraction 


describes the 


is accom- 
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plished, by modifying the solvent or by modi- 
fying the solute. 


2. The Absorption Spectrum as Affected by 
Chemical Changes 


We have interpreted the loss of fine structure 
and the gradual shift of the absorption bands to 
shorter wave-length as the result of increased 
solvation. Solvation must be pictured dynami- 
cally. The solute molecules attract and orient 
surrounding solvent molecules, but thermal mo- 
tion does not permit the formation of more 
permanent configurations. At low temperature 
especially if the forces acting between solvent 
and solute are strong, one may expect associa- 
tions of longer lifetime to form between solvent 
and solute. In this case one speaks of compound 
formation. lodine for example forms three groups 
of solutions. In hexane or carbon tetrachloride 
it produces the same purple as the vapor and 
one assumes little interference between the I. 
molecules and the solvent. For benzene and 
toluene we have found some interaction which 
causes the absorption band to be slightly shifted: 
these solutions are red. In alcohol, ether, or in 
an aqueous solution of potassium iodide, iodine 
dissolves with a deep brown color, and it is 
generally assumed that a more or less defined 
chemical compound has formed. In the potassium 
iodide solution it is the anion (I,)I~. Such an 
association of longer lifetime than the previously 
pictured dynamic solvation can give rise to 
different electron transitions.* 

Figure 6 indicates the change of the light 
absorption of iodine in hexane upon addition of 
ether. The formation of a compound between 
ether and iodine is indicated by the appearance 


TABLE II, 

Derivative Ether solution Alcoholic solution 
2,5-diamino ester yellow orange-yellow 
2-amino-5-oxyester green yellow-green 
2/5-dioxyester blue blue 
Monoacetyldiamino ester —_ blue-green green 


Diacetyldiamino ester violet-blue violet-blue 


* Whether a purple, red, or brown solution is formed 
depends not only on the nature of the solvent but also on 
the temperature. Red solutions become purple on heating 
because of decreased solvation. Brown solutions turn red 
because of the thermal dissociation of the labile compounds 


JOURNAL OF APPLIED PHYSICS 




















0.5+ 


VV 


0.254 











4000 5000 6000 A (A) 


Fic. 6. Light absorption of iodine in hexane-ether. 


urve 1: Hexane pure. 

urve 2: Hexane 2.28 percent ether. 
urve 3: Hexane 4.56 percent ether. 
urve 4: Hexane 7.15 percent ether. 
urve 5: Hexane 8.60 percent ether. 
urve 6: Hexane 28.6 percent ether. 
urve 7: Ether pure. 
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of a new absorption band. The position of the 
new band depends on the nature of the reactant. 
It is 450 my for ether and 440 my for alcohol. 
In contrast to the addition of benzene to the 
iodine hexane solution, we observe here the 
formation of a new color center on the expense 
of the original one which we attributed to free, 
but solvated, I. molecules. 

This type of change, appearance of a new 
spectrum and simultaneous disappearance of the 
original one, is characteristic for all chemical 
changes. It is widely used in colorimetry and in 
the field of glasses it has given us much valuable 
information about their chemical reactions. As 
an example for an oxidation-reduction process in 
glass, we give Fig. 7 demonstrating the reduction 
of CrO; to Cr2O; by means of As2O3. If chromium 
compounds are added to the glass batch, a green 
glass is formed. The color is produced by both 
Cr+ and Cr**+ which are in equilibrium. The 
position of the equilibrium depends upon the 
melting temperature, furnace atmosphere, com- 
position of the base glass and additions such as 
arsenic. The Cr** has its main absorption in the 
ultraviolet and in the short wave region of the 
visible. It does not absorb in the orange, where 
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Fic. 7. Light absorption of chromium containing glass, 
melted under oxidizing and under reducing conditions. 


the Cr*+ has an absorption maximum. The height 
of this band can, therefore, be used to evaluate 
the amount of Cr** and with it changes of the 
oxidation-reduction equilibrium. 

These chemical reactions are of great impor- 
tance to the glass technologist, especially to the 
one who has to make reproducible color filters. 
From the viewpoint of the glass structure another 
type of compound formation is far more signifi- 
cant, a change which cannot be expressed by the 
conventional chemical equations. 


3. The Spectrum as an Indicator of Changes 
in Coordination 


The different colors produced by iron, man- 
ganese, chromium, or vanadium result chiefly 
from these elements being capable of forming a 
number of different ions in glass. Iron occurs as 
Fe*+ and Fe*+; vanadium as V*+, V‘*+, and V*+; 
manganese as Mn** and Mn?*. Each ion is an 
individual and its optical properties are just as 
different as its chemical ones. 

Other elements such as cobalt and _ nickel 
when introduced into glasses can form only one 
type of ion, namely Cot* and Nit**. Nevertheless 
there are cobalt glasses which are pink and 
others which are deep blue. Nickel-containing 
glasses range from yellow over grey shades to 
deep purple. In these cases the different absorp- 
tion spectra are caused by different coordination 
complexes in which the Co** or Ni** participate 
as the central cation. This is known from the 
behavior of these ions in aqueous solutions and 
in various host lattices. 
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Fic. 8. Light absorption of the Co** ion in aqueous solutions 
containing increasing amounts of calcium nitrate. 


Cobalt ions diluted in aqueous solutions are 
pink. They have an absorption band in the green 
part of the spectrum. Color are the 
completely hydrated Cot*t* ions. Dehydration 
caused by increasing the temperature, addition 


centers 


of dehydrating compounds such as calcium 
nitrate or glycerol produces a slight change in the 
light absorption. The maximum is shifted to 
longer wave-length and its height 
Dehydration produces the opposite effect of 


increases. 


increased solvation. The electronic jump becomes 
more probable and is caused by light quanta of 
smaller size. Figure 8 represents the deepening 
of the pink of a diluted cobaltous nitrate solution 
upon addition of calcium nitrate. 

The addition of halides such as LiCl or CaCl. 
in its first stages produces the same effect. If the 
halide concentration exceeds a certain limit, one 
observes the formation of a deep blue complex 
ion such as (CoCl,)?~. This process is pictured 
in Fig. 9. 

Whenever an anionic complex CoX, is formed, 


one obtains an intensive light absorption in the’ 


long wave part of the visible spectrum. If X is 
Cl-, CNS~, or O*-, the color is deep blue. Using 
more polarizable ions such as I~ causes the bands 
to shift to longer wave-length and green com- 
plexes result. This has been observed in solutions 
as well as in borate glasses containing iodine. A 
blue pigment is obtained by introducing cobalt 
into a crystal lattice where it is surrounded by 
four oxygens. That is the case in the cobalt 
aluminate Al,CoO, (Thenards blue). Mixed 
crystals between MgO and CoO on the other 
hand are pink because they have NaCl structure 
and each Mg**+ and Co** is surrounded by six 
oxygens. This simple rule that CoO, complexes 
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are pink and CoQ, complexes are blue holds true 
only for those crystals which consist of rare gas 
ions. lons with incomplete outer electrons exert 
a strongly deforming influence upon the oxygen 
ions which in turn reflects upon the Co—O bond. 
Cobalt zincate for example is green. 

The same two types of color centers can be 
formed in glasses. CoO, groups can be considered 
complex anions comparable to the SiO, groups. 
It seems reasonable to assume that in the random 
three-dimensional network of SiO, tetrahedra an 
occasional SiO, group can be replaced by a CoO, 
group. In this case the cobalt ion plays the same 
role as the silicon. According to B. E. Warren 
ions with these functions are called ‘network 
forming” ions. They usually have the coordina- 
tion number four; only boron in the BOs glass 
has the coordination three. The polymerization 
of these randomly arranged XO, groups leads 
to a huge anionic three-dimensional network in 
whose interstices the ‘network modifying" ca- 
tions are distributed. 

The modern concept of N.W.F. and N.W.M. 
cations corresponds roughly to the older picture 
of acid and basic glass forming oxides. Only the 
acids such as SiQ., P2O;, BO; can form glasses 
by themselves. The glass chemist also used the 
term “amphoteric” for oxides which in some 
glasses assume the role of an acid and in others 
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Fic. 9. Light absorption of the Co** ion in aqueous solutions 
containing increasing amounts of lithium chloride. 


JOURNAL OF APPLIED PHYSICS 








that of a base. In modern parlance, Al**, for 
example, can act as an N.W.F. cation by entering 
the glass structure as AlO, tetrahedra or remain 
in interstices where it is surrounded by a larger 
number of oxygens in greater distance. 

The distinction between N.W.F. and N.W.M. 
cations is based on the different strengths of 
their electric fields. Small ions of high charge are 
the first to surround themselves with oxygen 
ions. Tetrahedral groups are particularly favor- 
able to glass formation. If the chemical compo- 
sition of the melt does not provide sufficient 
oxygen ions to permit the formation of inde- 
pendent SiO, groups, the tetrahedra have to 
share oxygens.* This polymerization increases 
the viscosity of the melt very rapidly and finally 
leads to a rigid three dimensional network, the 
glass. 

If the chemical composition of the melt pro- 
vides sufficient oxygens to form independent 
SiO, tetrahedra as in the case of sodium ortho- 
silicate, Na,SiOQ,, there is no more necessity for 
sharing the oxygens and with it for polymeriza- 
tion. Such melts are fluid and can crystallize 
easily on cooling. 

Whether or not a melt forms a glass depends, 
therefore, on the charge and size of the strongest 
cation and on their number as compared with 
that of the oxygen. V. M. Goldschmidt has de- 
rived certain rules for glass formation appli- 
cable to oxide and fluoride glasses, which were 
based on the probable coordination, because of 
size and charge of the cations and the cation to 
oxygen ratio. 

The validity of these rules is limited to the 
anions of the rare gas type. The strongly de- 
forming influence which ions with incomplete 
outer electronic shells exert upon the large anions 
leads to bonds which cannot be predicted any- 
more by simple rules, because the ionic radii of 
the deformed ions have lost their significance. 

The orthosilicates of noble gas ions do: not 
form glasses. However, lead and zinc orthosili- 
cates are known in vitreous condition. We have 
here a group of compounds which bridge the gap 
between the strictly ionic and the organic glasses. 


* The sharing of oxygens between two SiO, tetrahedra 
can be considered one principle of glass formation. It is 
the chief reason why silicate melts polymerize. 
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Whereas size and charge of the cobalt and 
nickel ions should classify them as N.W.M. 
cations or basic glass constituents, their incom- 
plete outer electronic shells enable them to form 
CoO, and NiO, complexes in equilibrium with 
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Fic. 10. Light absorption of the Ni** ion in a potassium 
silicate glass (purple type). 


those Co** and Ni*+ which have their places in 
interstices or N.W.M. positions. The relative 
number of XO, groups formed by their cations 
might be very small but they possess very 
intensive light absorption. The covalent bonds 
between Co?+ and O?- in the strongly bonded 
CoO, complexes probably makes possible hori- 
zontal electron transitions, producing electron 
transfer spectra instead of the prohibited vertical 
transitions of Co?* ion. 

Figure 10 shows the absorption spectrum of 
a purple Ni’ containing silicate glass. We at- 
tribute this characteristic absorption to NiO, 
groups. They form preferably in those glass 
compositions where the Ni?+ can successfully 
compete with the other cations and form 
independent anionic groups (rubidium silicate). 
Glasses built up from ions with higher electric 
potential (lithium silicate or phosphate glasses) 
prevent the formation of NiO, groups and shift 
the Ni?+ into N.W.M. or interstitial positions. 
In this case the Ni®+ has a higher coordination, 
maybe 6, is farther away from the anions, and its 
light absorption is different (Fig. 11). Glasses of 
this type are yellow. The Ni’ ion is thus an 
indicator for the availability of anions in the 
glass. Its light absorption is a function of the 
extent to which it succeeds in forming its own 
anionic or covalent group. A few results obtained 
with Ni** as color indicator shall bé mentioned. 
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Fic. 11. Light absorption of the Ni#* ion in a lithium- 
sodium silicate glass (yellow type). 


(1.) Sodium silicates containing Ni** are grey. 
Their absorption spectra show the presence of 
both NiO, and NiOg groups. Replacing Na* by 
the smaller and more effective Lit 
causes the color to change to yellow. Substitution 
of K*, and even more so of Rb*, for Na* produces 
purple glasses. The weaker the electric field of 
the other N.W.M. cations, the greater are the 
chances for the Ni** ions to compete for their 
oxygens and form NiO, groups. 

(2.) If the gray Ni** containing sodium silicate 
glasses are heated into the softening range, their 
light absorption indicates an increase in the 
number of NiO, groups (Fig. 12). On cooling to 
room temperature some of the nickel ions change 
back from N.W.F. to N.W.M. positions. The 
indicator ion signals a redistribution of the anions 
with temperature, not merely a change of the 
interionic distances. In the crystalline state such 
a change of geometrical arrangement of the ions 
produces a new modification. It has to take 
place abruptly because of the symmetry require- 
ments of the lattice and it is accomplished by 
cooperative maneuvers. Glasses show changes in 
their structure which are comparable with modi- 
fication changes of crystalline silicates but the 
adjustment takes place gradually; the single 
units do not have to cooperate to accomplish 
the change. 

(3.) The light absorption of a chilled Ni** 
containing sodium silicate glass indicates the 
presence of more NiO, groups than in the an- 
nealed condition (Fig. 13). This indicates that 
glasses of different heat treatment possess differ- 
ent atomic arrangement. 

The fact that glasses of different heat treat- 


therefore 
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Fic. 12. Light absorption of the Ni?* ion in a sodium 
silicate glass of different temperatures. 


ment possess very different physical and chemical 
properties can be explained on the same basis as 
the difference of properties of the high and low 
temperature form of a certain crystalline com- 
pound. Now we can understand the meaning of 
the volume contraction which the quenched 
glass exhibits on heating. It corresponds to a 
change of the high temperature modification 
into the stable low temperature form. The rapid 
chilling of the glass had frozen in a high temper- 
ature form just as the cooling of molten sulphur 
produces the rhombic rather than the monoclinic 
form. 

The change in coordination or the role which 
a cation plays in the glass structure effects the 
emission spectrum in a very similar way as its 
absorption. The divalent manganese is one of 
the best examples. MnO, groups produce green 
and Mn** in N.W.M. positions produces red 
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Fic. 13. Light absorption of a nickel-sodium silicate glass 
in annealed and chilled conditions. 
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fluorescence. The principles involved in using 
Mn? ions as fluorescence indicators are mostly 
the same as described for Nit+* as color indicator. 
N. J. Kreidl* has recently published a paper on 
this subject where he describes the use of the 
Mn** for exploring new fields of optical glasses. 

It might be interesting, however, to mention 
that this fluorescence indicator had been used to 
draw conclusions concerning the structure of the 
yellow fluorescing zinc silicate. This compound 
represents an unstable modification of the zinc- 


‘N. J. Kreidl, J. Opt. Soc. Am. 35, 249-257 (1945). 


orthosilicate activated by Mn which for a long 
time could not be obtained in well-developed 
crystals. It forms when ZnO and SiO; are allowed 
to react at relatively low temperature and by 
some authors it was considered to be amorphous. 
Only recently we could prepare this modification 
in better crystals as a devitrification product of 
some multicomponent zinc silicate glasses. These 
glasses precipitate zinc orthosilicate on heating 
and form the stable green modification at high 
temperature and according to ‘‘Ostwald’s Step 
Low” the yellow unstable modification at lower 
temperature. 
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30-October 3 American Society of Mechanical Engineers, Boston, 


Massachusetts 


VOLUME 17, JULY, 1946 


October 


!3- 5 National Electronics Conference, Inc., Chicago, Illinois 


9-11 American Society of Civil Engineers, Kansas City, Missouri 
25 American Mathematical Society, New York, New York 
November 
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and Industrial Research, National Research Council, Balti- 
more, Maryland 


29-30 American Physical Society, Minneapolis, Minnesota 


December 

2-6 American Society of Mechanical Engineers, New York, New 
York 

17 Institute of Aeronautical Sciences, New York, New York 

26-28 Geological Society of America, Chicago, Illinois 
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Here and There 








New Appointments 


N. B. Nichols has been appointed director of the re- 
search division of the Taylor Instrument Companies. 
Carlton W. Miller and William A. Mersman, both formerly 
of the Radiation Laboratory, have joined the research 
department, and Walter G. Countryman has returned to 
research work at Taylor after four years of military service. 

Radio Corporation of America recently announced two 
appointments. John L. Reinartz has returned to the com- 
pany, to be in charge of the amateur radio program in the 
Commercial Engineering and Power Tube Sections. John F. 
Rider, radio engineer and writer, has been retained by the 
Victor Division as a consultant on test equipment. 


Fellowship at Michigan State 


Michigan State College has announced a $1000 fellow- 
ship available for graduate study and research in physics 
for one year from September 1946. This fellowship has been 
established from a Frederick Gardner Cottrell special 
grant-in-aid, contributed by The Research Corporation of 
New York, for a study of some of the properties of ions and 
electrons in gases. Applications will be considered from any 
able individuals who hold bachelors’ or masters’ degrees 
from approved institutions and whose work shows promise 
of originality. The Fellow will not be required to teach. 
Inquiries should be addressed to Dr. R. C. Huston, Dean 
of the Graduate School, Michigan State College, East 
Lansing, Michigan. 


Laboratory Clinics on Polymer Research 


The Institute of Polymer Research at the Polytechnic 
Institute of Brooklyn is sponsoring a series of summer 
laboratory clinics for 1946. Five sessions will occur, June 
through September. A fee is charged for each clinic, and 
attendance for each one will be limited. 


Journal of Polymer Science 


Journal of Polymer Science is now being published bi- 
monthly by Interscience Publishers, Inc., and Elsevier 
Publishing Company, Inc. It is devoted to the advance- 
ment of our fundamental knowledge of the physics and 
chemistry of polymers. The new journal is an outgrowth 
of, and a successor to, the Polymer Bulletin, which has 
ceased to appear. It reflects the consolidation of independ- 
- ent plans of the Elsevier Publishing Company, Inc., of 
Amsterdam, and Interscience Publishers, Inc., of New 
York, to publish separate journals on the same topic. The 
Journal will be international in scope. All papers will be 
published in the English language. 


Duke University Fellowships 
Duke University is offering ten special research fellow- 
ships in microwave electronics and molecular physics 
with stipends from $1600 to $2400, Graduates of recog- 
nized colleges who have majored in physics are eligible. 
Some of the appointments are open to persons with only 


640 


an undergraduate minor in physics, provided they have a 
major in chemistry or an equivalent training in electronics. 
Individuals interested in these appointments may write 
for particulars to the Chairman, Department of Physics, 
Duke University, Durham, North Carolina. 


Conference at Texas A. & M. 


Research personnel in industry and college teachers were 
invited to attend a graduate conference at College Station, 
Texas, devoted to the behavior of mass spectrometers, 
electron microscopies, and other electronic devices whose 
operation depends on the action of electron or ion beams. 
The conference was conducted by the Electrical Engineer- 
ing and Physics Department of the Agricultural and Me- 
chanical College of Texas Juné 24 through July 6. The lec- 
turers were Dr. Ladislaw Marton of the Division of 
Electron Optics of Stanford University and Dr. John A. 
Hipple of the Westinghouse Research Laboratories. 


Film on Compressed Air 


“Our Industrial Air Power”’ is the title of a new educa- 
tional film on the many industrial uses of compressed air 
which is available for showing before students, engineering 
societies, and other groups. Among the sequences in the 25- 
minute, 16-mm sound-color film are: ‘‘Fundamentals of 
Air Compression,’ ‘Properties of Compressed Air,”’ “‘How 
Different Types of Compressors Work,” and “How Com- 
pressed Air Is Used.”’ The film is available without charge 
on a loan basis. Write Quincy Compressor Company, 
Quincy, Illinois. 


Army Research and Developments Program 


Early in May the creation of a Research and Develop- 
ments Division as a top-level General Staff organization to 
coordinate Army research with the activities of industry 
and educational institutions was announced by Secretary 
of War Robert P. Patterson. At that time General Dwight 
Eisenhower, Chief of Staff, had completed plans for the 
operation of the Division on a level with the Personnel 
Intelligence, Training, Supply, and Operations section of 
the General Staff. 

The Division has primary interest in the application of 
national scientific resources to the solution of military 
problems. The director will act as liaison between Army 
planning chiefs on one hand and industry and the uni- 
versities and research the other. Great 
civilian advances are foreseen as a result of much of the 
Army-sponsored research. 


laboratories on 


Cosmic-Ray Investigation 


In search of fuller information in regard to cosmic-ray 
intensities as they vary with latitude and altitude, the 
National Geographic Society, the United States Army Air 
Forces, and the Bartol Research Foundation jointly con- 
ducted in May a series of four round-trip flights in a spe- 
cially equipped B-29 bomber. The principal apparatus 
carried in the plane consisted of multiple banks of Geiger 
counters, designed by Dr. W. F. G. Swann, Director of the 
Bartol Research Foundation. 
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New Books 


Physics of the Twentieth Century 


By PAscuaL JoRDAN. Pp. 185+xii. Philosophical Li- 
brary, New York, 1944. Price $4.00. 








This is a simple little book, addressed to the non- 
physicist who is curious about the results of modern 
physics, and intended to instruct the philosopher. Its 


merit lies in its wealth of interesting examples designed © 


to clarify the reputedly paradoxical teachings of relativity 
theory and quantum mechanics. The difference in outlook 
between classical mechanics and modern theory is very 
impressively sketched. Unfortunately, crudities of transla- 
tion and limitation of scope render it inferior, on the whole, 
to the elegant expositions of Eddington. 

The interpretation of quantum mechanics, to which a 
large part of the book is devoted, involves an unreserved 
avowal of the early view which saw in the new doctrine a 
mathematical reflection of the uncontrollable effects of 
measurement on physical objects, an establishment of 
fundamental discontinuity in nature, and a thorough re- 
jection of causality. Very little attention is paid to the 
formulation, implicit in the work of Dirac and Schrédinger, 
according to which the concept of “‘physical states’’ has 
undergone modification. As a result, comments on crucial 
matters are, in the reviewer's opinion, sometimes super- 
ficial. Thus when the charge distribution of the electron 
in a hydrogen atom is asserted to owe its diffuseness to the 
unpredictable bombardment by gamma-ray quanta coming 
from the famous microscope, no explanation is given as to 
why the distribution should not depend on the frequency 
of the gamma-rays. 

As a philosopher, Jordan seems to urge two things: a 
recognition that materialism has been liquidated, and a 
bright-eyed positivism willing to accept all empirical facts 
even if they spell religion. His deliberations dealing with 
the theory of knowledge do not, it seems, rise to the level 
of technical competence that would impress philosophers. 

In spite of its faults, the book makes pleasant reading. 
Its misfortune is that it tries to enter a market—that of 
English scientific literature—in which the standard of 
excellence for works of this sort is extremely high, 

HENRY MARGENAU 
Yale University 





New Booklets 








School of Mineral Industries, The Pennsylvania State 
College, recently published its Circular 20, Opportunities 


for Productive Work through Mineral Industries Research, 


a 36-page pamphlet generously illustrated with photo- 
graphs and drawings. Free on request. Address the School 
at State College, Pennsylvania. 


American Optical Company, Scientific Instrument 
Division, Buffalo 11, New York, recently issued a four- 
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page folder describing the Spencer Dark Field Quebec 
Colony Counter which has been placed on the market. 


Televiso Products Company, 7466 West Irving Park 
Road, Chicago 34, Illinois, has published Bulletin No. 31, 
a four-page description of the Series 200 Vacuum Tube 
Voltmeter. 


Scientific Development Corporation, 614 West 49 
Street, New York 19, New York, began in May a new 
publication service entitled Public Domain. Each weekly 
issue will be a 250-page reference book containing over 
1000 patents due to expire four weeks after date of issue, 
which will then be in the public domain, and open to 
exploitation by whoever cares to use them. Each patent 
shown will include a good reproduction of a draftsman’s 
drawing together with a digest of typical claims and salient 
features. Charter subscriptions are offered for one year 
at $45; for six months at $25; and for 10 weeks at $10. 


Acme Scientific Company, 200 North Laflin Street, 
Chicago 7, Illinois, has issued its Bulletin 453, which gives 
an extensive listing of stock sizes of optical flats, both 
single and double surfaced. 


Nickel Steel Topics, published at 67 Wall Street, New 
York 5, New York, discusses in its March, 1946 issue the 
use of nickel alloys in tractor coil springs, aircraft engines, 
oil well drilling equipment, cutting saws, and testing 
machinery. 


Precision Equipment Company, 32 North State Street, 
Chicago 2, Illinois, has published a new catalog for all 
factory superintendents, maintenance engineers, purchas- 
ing agents, and other plant executives interested in in- 
dustrial equipment and supplies. Any items of industrial 
equipment or supplies not listed in this catalog will be 
procured by Precision’s Industrial Service Department. 


The Tin Research Institute, Fraser Road, Greenford, 
Middlesex, England, has issued a 32-page report of its 
work covering the years 1942-1944. Free on request. 
Technical representation for the Institute in the United 
States is furnished by the Battelle Memorial Institute, 505 
King Avenue, Columbus, Ohio. 


Castle Films, Inc., 30 Rockefeller Plaza, New York, 
New York, is offering free to users of training films a new 
catalog listing for the first time all U. S. Government 
16 mm sound films and 35 mm film strips, totaling 1158. 
They include training and other types of educational films 
produced by the Navy, War Department, Department of 
Agriculture, and U. S. Public Health Service. 


The March and April Ohmite News, published by Ohmite 
Manufacturing Company, 4835-41 Flournoy Street, Chi- 
cago 44, Illinois, featured a biographical sketch of Sir 
Joseph John Thomson, discoverer of the electron, and a 
description of the new Ohmite Laboratory for Precision 
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Illinois Institute of Tech- 
nology. Two pages. Free on request. 


Electrical Measurements at 


The Frontier is published quarterly by Armour Research 
Foundation of Illinois Institute of Technology, Chicago 
16, Illinois. The March 1946 issue features articles on 
“Limitations of Hardenability Indices,”’ “Recent De- 
velopments in the Use of Mold Bran,” “Electron Micros- 
copy,” and “Drinking Water Made Palatable by Chlorine 
Dioxide.” Twenty pages, illustrated. 


The General Radio Experimenter for April 1946, pub- 
lished by General Radio Company, 275 Massachusetts 
Avenue, Boston 39, 
tronic sequence timer which permits multiple “exposures 
on a single film to be made with the Microflash. Eight 
pages, illustrated. 


Massachusetts, describes an elec- 


The Main Library of General Electric Company, 
Schenectady, New York, issues monthly an eight-page 
review of recent periodicals, giving the title, source, and 
sometimes a very brief abstract, of current technical arti- 
cles. The leaflet is called Library Service, College Edition. 


A twenty-page pamphlet is available entitled Industrial 
Research Progress at Armour Research Foundation of IIli- 
nois Institute of Technology, 1944-45, describing its facilities 
and activities in physics research, as well as those in the 
other major fields of chemistry, metallurgy, engineering, 
etc. The Foundation is a non-profit institution established 
in 1936 to render a research and experimental engineering 
service to industry. Copies may be obtained by addressing 
Armour Research Foundation, Chicago 16, Illinois. 


Leeds and Northrup Company, 4934 Stenton Avenue, 
Philadelphia 44, Pennsylvania, has issued Catalog ED, 
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Galvanometers and Dynamometers. This is a new, completely 
rewritten edition of earlier Leeds and Northrup galvanom- 
eter catalogs. 37 pages. Free on request. 


The Gaertner Scientific Corporation, 1201 Wrightwood 
Avenue, Chicago 14, Illinois, has issued Bulletin 154-64, 
which describes Gaertner spectroscopes of the Bunsen- 
Kirchhoff type. Four pages. Free on request. The same 
company has also published Bulletin 153-64 on automatic 
linear dividing machines. Twenty pages. Free if requested 
on company letterhead. 


North American Philips Company, Inc., 100 East 42d 
Street, New York 17, New York, has released a new four- 
page booklet (R1022) entitled “‘Geiger-Counter Used in 
Powder Metallurgy.” Free on request. 


Electro-Voice, Inc., 1239 South Bend Avenue, South 
Bend 24, Indiana, has issued an illustrated four-page 
bulletin on its new Model 950 Cardax, a cardioid uni- 
directional crystal microphone. Free on request. 


Interchemical Review, winter issue of 1945, published by 
Interchemical Corporation, 432 West 45th Street, New 
York 19, New York, has the following table of contents: 

Survey of Sorption Chemistry 


Our Part in the War Effort 
A Tour of Bombed Nazi War Plants 


EDMUND N. HARVEY, JR. 
Joun O'CONNOR 
ERNEST W. PITTMAN 


Bausch and Lomb Optical Company, Rochester 2, New 
York, in its spring 1946 issue of The Educational Focus 
presents the following articles: 

Optics in Scientific Crime Detection J. EpGarR Hoover 

Color Photomicrography J. V. BUTTERFIELD 


Uses of Optical Equipment in the Aluminum Industry 
STANLEY V. MALculIt 
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